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FOREWORD 


This  :.eport,  which  presents  the  experimental  and  theoretical  results 
of  an  investigation  of  aerotheraoelastio  instabilities  of  inflaJ  >ble 
vehicles,  was  prepared  by  the  Aeroelastic  and  Strvei.ures  Rea  ear  ci.  Labora¬ 
tory,  Massachusetts  Institute  of  Technology,  Cambridge,  Massachu;  atta  for 
the  AF  Flight  Dynamics  Laboratory,  Research  and  Technology  Division,  Vright- 
Patterson  Air  Foroe  Base,  Ohio.  The  major  portion  of  the  work  was  performed 
at  MIT  under  the  supervision  of  Professor  H.  Ashley,  with  Mr.  J.  R. 
M&rtuccelli  and  Mr.  0.  Zartarlan  as  co-project  leaders.  The  flutter  tests 
were  carried  out  in  Ihnnel  A,  Oas  Dynaaioa  Facility,  ARO  Ino.,  AEDC, 
Tullahcma,  Tenn.,  the  transonio  tunnel  at  the  Cornell  Aeronautical  Labora¬ 
tories,  Inc.,  Buffalo,  New  Tori',  and  the  Wright  Brothers  Wind  Tunnel, 
Massachusetts  Institute  of  Technology,  Cambridge,  Massachusetts. 

The  research  was  accomplished  under  Air  Foroe  Contract  No.  AF33(657) - 
8L68,  Project  No.  1370,  "Dynamic  Pro  hi  sms  in  Flight  Yehialee",  and  Task 
No.  137003,  "Prediction  and  Prevention  of  Aaro therm oelastlo  Instabilities", 
with  Mr.  S.  J.  Pollock  of  the  Flight  Dynamic*  Laboratory  as  task  engineer. 
This  research  was  started  in  June  1562  and  was  completed  In  August  1962*. 

Authors  of  the  main  report  are  J.  R.  Kartu ocelli,  0.  Zartarlan,  and 
R.  B.  MeCallum.  Appendices  H  and  HI  are  mainly  the  work  of  Professors 
M.  L  an  da  hi  and  J.  W.  Mar,  respectively.  The  authors  with  to  express  their 
appreciation  to  Professor  J.  Dugundji  for  his  guidance  In  the  development 
of  the  vibration  theory  for  inflatable  vehicles,  and  to  the  personnel  of 
the  wind  tunnel  facilities  used  In  the  oourae  of  this  program  for  their 
help  in  conducting  the  flutter  teste. 

Portions  of  this  document  are  classified  CONFIDENTIAL  since  the  data 
revealed  can  be  employed  to  establish  design  orlterla  for  the  prevention 
of  aeroelaatlc  instabilities  of  future  Inflatable  vehicles.  This  report 
has  been  given  the  Aeroelaetio  and  Structures  Research  Laboratory  number 
ASKL  TO  lli*-l. 
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ABSTRACT 


This  report  covers  a  program  of  research  on  analytical 
studies  and  experimental  verification  of  the  vibration  and 
flutter  characteristics  of  inflatable  structures,  as  a  first 
step  toward  establishing  design  criteria  which  will  insure 
that  lifting  surfaces  for  future  inflatable  vehicles  will 
be  free  from  dynamic  instabilities.  The  research  effort 
was  mainly  confined  to  the  AIRMAT-type  structure,  where  the 
spacing  of  the  surfaces  is  achieved  by  chords  held  in  ten¬ 
sion  by  internal  pressure. 

A  method  is  presented  for  constructing  A1RMAT  models 
which  are  suitable  for  use  as  wind  tunnel  flutter  models. 

In  addition  to  the  standard  vibration  data,  the  vibration 
mode  shapes  of  a  delta  wing,  as  determined  by  use  of  mirrors 
mounted  on  the  model,  are  presented.  Shear-  and  bending- 
type  deformations  are  shown  to  be  equally  important  for  the 
models  used  in  the  programs  by  comparing  the  results  of 
static  and  dynamic  tests  with  theoretical  calculations. 

Some  experimental  results  on  the  properties  of  the  model 
surface  material  are  also  included. 

Sixty-five  flutter  points  were  obtained  for  cantilevered, 
half-span,  delta-wing  models  over  the  Mach  number  range  from 
0.116  to  6.0.  In  addition  to  determining  an  experimental 
boundary  for  a  basic  delta  model,  the  effects  of  leading- 
edge  sweep,  angle  of  attack,  and  the  pressure  differential 
between  the  model  and  ambient  air  pressures  were  investigated 
at  several  supersonic  Mach  numbers. 

A  series  of  theoretical  calculations  were  carried  out, 
employing  various  combinations  of  mode  shapes  (based  on 
eeasured  vibration  node  lines).  The  predicted  flutter 
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UNCLASSIFIED 


boundaries  are  to  fair  agreement  with 
The  lack  of  better  agreement  can  be  act.. a-  e 
plifying  assumptions  made  in  the  analysis,  particular  y 
those  dealing  with  the  assumed  code  shapes. 

A  theoretical  investigation  on  membrane-type  flutter 
of  an  Inflatable  panel  and  the  development  of  the 
governing  the  structural  behavior  of  an  Inflated  MSMAT  etc 
tion  are  presented  in  separate  appendices. 
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Appendix  I 
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Appendix  1 
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Pressure  coefficient,  defined  as 
/■■l  P  Pea 


Stagnation  point  pressure  coefficient 

Cylinder  model  diameter 

Axial  and  hoop  extensional  moduli, 
respectively 

Extensional  modulus 

Cover  "extensional  stiffness" 

Shear  modulus 

Cover  "shear  stiffness" 

Mach  number 

Mass  of  air  contained  in  the  seraichord 
formed  by  rotating  the  delta  model  about 
its  root  chord 

Generalized  masses  associated  with  modes  1 
and  j 

Generalized  mass  for  the  n-th  orthogonal 
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i-th  generalized  external  force 
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modes  i  and  according  to  slender-body 
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Modified  generalized  aerodynamic  forces; 
see  Section  B,  Appendix  1 
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edges 


xix 


SYMBOLS  (Coat'd) 

Temperature  in  degrees  Fahrenheit,  unless 
specified  otherwise;  also  torque  applied 
to  cylinder  models  (see  Section  2);  also 
kinetic  energy,  Eq-  (1-19) 

Strain  energy,  Eq-  (1-18) 

Free-stream  velocity 

Weight  added  to  cylinder  model  to  equalize 
hoop  and  axial  stress  levels 
Elementary  modal  functions  used  in  the 
transonic  flutter  analyses;  see  Section  B, 
Appendix  1 

Shear  angles  defined  in  Fig.  1*1 
Angle  of  attack;  angle  of  attack  at  root 
section  in  Section  3 

Initial  steady  constant  angle  of  attack 
Steady  and  unsteady  parts  of  oc  ;  see 
Section  C,  Appendix  I 

Adiabatic  exponent  (“  1**  for  air^;  also 
angle  related  to  twist  of  cylinder  models 

(see  Section  2) ,  • 

Entry  angle;  see  Fig*  1*3 

Axial  and  hoop  strains  in  the  cylinder 

models,  respectively 
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defined  in  Section  A,  Appendix  I 
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bration  frequencies 
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Hass  ratio  parameter,  defined  as^*-“~~-7p 

Also  Poisson’s  ratio,  when  latter  ia  dis¬ 
cussed  generally,  or  for  isotropic  materi 
als;  accompanied  by  subscript  when  referred 
to  specifically  for  nonisentropic  materials 

Poisson's  ratio  relating  hoop  strain  to 
axial  stress 

Poisson's  ratio  relating  axial  strain  to 
hoop  stress 

Poisson's  ratio  relating  warp  strain  to 
fill  stress 


Poisson's  ratio  relating  fill  strain  to 
warp  stress 

Poisson's  ratio  relating  strain  along  a 
direction  at  45  degrees  to  the  warp 
(or  fill)  to  a  stress  alc;\g  a  perpendicular 
to  this  direction 


n-th  generalized  coordinate 


Free- stream  density 


Free-stream  density 

Tangent  of  the  semiapex  angle;  also  the 
cotangent  of  the  leading-edge  sweep;  also 
semispan  made  dimensionless  with  respect 
to  root  chord,  see  Fig.  1.2 

Axial  and  hoop  stresses  in  the  cylinder 
models,  respectively 


Stress  level  in  the  model  cover,  see 
Section  2 
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SYMBOLS  (Cont'd) 

Complex  amplitude  of  the  velocity  potential 
associated  with  elementary  modal  function  Z; 
superscript  (1)  denotes  slender-body  so¬ 
lution.  When  with  subscript  x"  ,  it  denotes 
differentiation  of  with  respect  to  X 

n-th  orthogonal  mode  shape 


Frequency  of  oscillation  ir.  rad/sec 

Frequency  in  rad/sec  associated  with  free 
vibration  mode  n. 


Leading  edge  sweep  angle 


Superscripts,  Subscripts,  etc. 
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n 
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Denote  the  particular  set  of  mode  shapes 
to  be  considered 

Dot  indicates  differentiation  with  respect 
to  time 

Refers  to  fill  direction;  also  flutter 
condition  .. 

Refers  to  warp  direction 

On  ci  ,  /3  ,  ,  they  indicate  differen¬ 

tiation  with  respect  to  x,  y  (or  x,  y) 
in  Section  A,  Appendix  I 

Refers  to  static  divergence 

Refers  to  a  direction  at  45  degrees  to 
the  warp  (or  fill) 
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SECTION  1 
INTRODUCTION 


Among  the  many  configurations  proposed  for  re-entry  ve¬ 
hicles,  the  concept  of  an  inflatable  structure  (Refs.  1,  2,  3) 
offers  some  important  advantages  which  bear  directly  on  the  prob¬ 
lems  encountered  during  re-entry.  *ne  configurations  represents^ 
tive  of  the  type  likely  to  be  used  for  recoverable  inflatable 
vehicles  are  shown  in  Fig.  1.1.  The  basic  features  of  these 
configurations  are  a  hemispherical  nose,  a  swept-back  leading 
edge  (usually  a  large  diameter  tube),  and  a  sail,  AIRMAT,#  or 
several  tubes  connecting  the  right  and  left  leading  edges 
(see  Fig.  1.2). 

Two  trajectories,  which  are  typical  for  the  types  of  re¬ 
entry  vehicles  shown  in  Fig.  1.1,  are  presented  in  Fig.  1.3, 
in  terms  of  dynamic  pressure  vs  Mach  number.  The  additional 
parameters  T(temperature) ,  h(altitude),  a(acceleration)  and 
t(time  from  initiation  of  entry)  are  also  given  since  they  in¬ 
fluence  the  final  design  of  an  actual  vehicle.  The  trajectories 
were  obtained  from  Ref.  1  and  some  additional  calculations  based 
on  the  1959  USAF  model  atmosphere.  The  angle  of  attack,  ot  ,  and 
the  lift  drag  ratio,  l/D  -  0.577,  are  assumed  constant  during 
entry.  The  entry  height  is  400,000  feet,  the  initial  flight  path 
angle  is  -  1  degree,  and  the  leading  edge  radius  is  taken  as  one 
foot. 

From  the  foregoing  and  other  published  studies  of  inflat¬ 
able  re-entry  vehicles  (Refs.  2,  4)  one  can  draw  the  following 
general  observations  which  are  applicable  to  most  situations: 

'Trademark  of  Goodyear  Aerospace  Corporation 

Manuscript  released  by  the  authors  August  1964  for  publication 
as  an  FDL  Technical  Documentary  Report. 

S4/3-220: 
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(1)  With  the  light  wing  loadings  required  (W/S  of  the 
order  of  unity),  the  maximum  dynamic  pressures  will 
occur  at  relatively  high  Mach  numbers,  say  2  <  M  <  10, 
at  altitudes  in  the  neighborhood  of  200-250,000  feet. 

(2)  The  variation  in  dynamic  pressure,  for  the  range  of 
Hath  number  in  Fig.  1.3,  is  less  than  20  per  cent. 

(3)  While  the  descent  will  be  at  very  shallow  angles, 
the  lifting  surfaces  will  be  oriented  at  relatively 
large  angle  of  attack  (between  20  and  60  degrees), 
because  of  the  low  dynamic  pressures  which  will  be 
encountered  at  these  higher  altitudes.  The  minimum 
C.  will  be  dependent  on  the  maximum  tolerable 
leading-edge  temperature. 

(4)  Maximum  allowable  heating  conditions  dictate  blunt 
leading  edges  for  such  surfaces. 

(5)  The  lower  wing  loadings  give  lower  values  of  the 
radiation  equilibrium  stagnation-point  temperature 
because  lower  speeds  are  encountered  at  any  given  al¬ 
titude. 

(6)  The  very  compact  packaging  of  the  vehicle  for  launch¬ 
ing  into  orbit  and  assembly  in  space  removes  the  ne¬ 
cessity  of  having  to  place  a  large  lifting  surface 
vehicle  on  top  of  a  large  rocket  vehicle,  and  thereby 
puts  the  task  of  placing  such  a  vehicle  in  orbit  withii 
the  capabilities  of  present  day  boosters.  A  packaged 
volume  to  inflated  volume  ratio  of  0.03  to  0.05  is 
possible  for  an  AIRMAT  wing,  while  for  a  complete  re¬ 
entry  vehicle,  ratios  of  0.15  are  possible  (Ref.  2). 

The  hish  degree  of  flexibility  exhibited  by  these  inflat¬ 
able  structures  indicates  that  the  prevention  of  instabilities 
such  as  flutter  may  be  a  significant  problem.  This  report  covers 
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a  program  of  research  on  analytical  studies  and  experimental 
verification  of  the  vibration  and  flutter  characteristics  of  in¬ 
flatable  structures,  as  &  first  step  toward  establishing  design 
criteria  which  will  insure  that  lifting  surfaces  for  future  in¬ 
flatable  vehicles  will  be  free  from  dynamic  instabilities. 

The  research  effort  was  mainly  confined  to  the  AIRMAT- 
type  structure,  where  the  spacing  of  the  surfaces  is  achieved  by 
chords  held  in  tension  by  internal  pressure.  With  regard  to  the 
type  of  models  that  should  be  employed,  the  first  question  which 
naturally  arises  is  the  amount  of  effort  which  should  be  expended 
in  attempting  to  achieve  dynamic  similarity  between  the  wind 
tunnel  model  and  a  full-scale  vehicle.  A  strict  adherence  to 
the  requirements  for  dynamic  similitude,  as  found  from  dimension¬ 
less  analysis,  is  most  often  precluded  by  lack  of  suitable  con¬ 
struction  materials  and  the  properties  of  the  wind  tunnels  used 
in  the  test  program.  It  was  decided  that  the  goals  of  the  pro¬ 
gram  could  best  be  served  by  adopting  the  viewpoint  that  the 
models  should  provide  a  means  for  checking  theoretical  calcula¬ 
tions  based  on  model  parameters,  as  opposed  to  the  use  of  models 
to  check  or  sioulate  a  full-scale  vehicle.  At  the  least,  the 
models  should  have  the  same  planfora  and  the  same  type  of  in¬ 
ternal  construction  os  the  full-scale  vehicle. 

Section  2  contains  the  model  information,  which  includes 
s  description  of  the  model  construction  and  a  complete  summary 
of  all  the  preliminary  static  and  vibration  tests  carried  out  on 
the  models  prior  to  the  wind  tunnel  tests.  The  relative  iopor- 
tanr.e  of  shear-  and  bending -type  deformations  are  determined  by 
comparing  experimental  results  with  some  simplified  theoretical 
calculations.  The  initial  attempts  at  measuring  the  structural 
properties  of  the  surface  material  are  also  fully  described, 
along  with  a  discussion  of  the  problems  encountered  and  a  pre¬ 
sentation  of  results  obtained  using  some  new  testing  techniques. 
A  preliminary  account  of  some  of  the  material  presented  in 
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Section  2  can  also  be  found  in  Ref.  5. 

Section  3  contains  a  description  ox  the  wind-tunnel-test 
procedures  and  the  results  and  discussion  of  all  the  wind-tunnel- 
flutter  tests.  During  the  program,  a  total  of  65  flutter  points 
was  obtained  from  the  21  flutter  models  which  were  constructed. 

The  conclusions  and  recommendations  drawn  from  the  in¬ 
vestigation  comprise  Section  4. 

The  analytical  procedures  and  summaries  of  the  theoretical 
calculations  are  presented  in  Appendix  I.  Also  included,  is  a 
discussion  on  the  effects  of  leading-  and  trailing-edge  bluntness 
and  large  angles  of  attack  on  the  unsteady  aerodynamic  forces. 

Appendix  II  presents  some  theoretical  work  on  membrane- 
type  flutter  of  an  inflatable  panel,  which  is  principally  the 
effort  of  Professor  M.  T.  Landahl  of  the  Department  of  Aero¬ 
nautics  and  Astronautics,  MIT. 

Appendix  III  contains  the  development  of  the  equations 
governing  the  structural  behavior  of  an  inflated  AIRMAT  section. 
This  is  the  work  of  Professor  J.  W.  Kar  of  the  Department  of 
Aeronautics  and  Astronautics,  MIT. 

Appendices  II  and  III  are  presented  as  separate  entities, 
and  each  contains  its  own  list  of  symbols  and  references. 

An  extension  of  this  work  up  to  Mach  numbers  8  and  10  i* 
being  sponsored  by  the  Flight  Dynamics  Laboratory,  Research  and 
Technology  Division,  Wright-Patterson  Air  Force  Base,  Ohio,  under 
U.S.  Air  Force  Contract  No.  AF  33(615)-il55. 
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Tubular  -  Cylinder  -  Hemisphere  Construction 


AIRMAT  -  Hemisphere  Construction 

Figure  1.1  Sketches  of  Typical  Blunted  Delta  Inflatable 
Vehicles,  Showing  Type  of  Construction. 
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Fleur*  1.3  Typical  Re-Entry  Trajectory  for  an  Inflatable  Vehicle  y,  -  —  1  deg 
(Entry  Angle),  L/D  “  0.577,  “  400  K  ft. 
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SECTION  2 
MODEL  INFORMATION 


A.  MODEL  CONSTRUCTION  AND  INSTRUMENTATION 


The  AIRMAT  type  of  lifting  surface  construction,  where 
the  spacing  of  the  surfaces  is  achieved  by  chords  held  in  tension 
by  internal  pressure,  was  conceived  and  is  under  development  by 
Goodyear  Aerospace  Corporation,  Akron,  Ohio.  Special  looms  have 
been  developed  which  weave  the  two  surfaces  and  the  drop  chords 
in  one  integral  operation.  The  elastomer,  which  makes  the  ma¬ 
terial  surface  impermeable  to  air,  is  applied  after  the  weaving 
operation.  Variable  depth  mats  can  also  be  formed  on  these 
looms .* 


Because  of  the  relatively  small  number  of  models  required 
for  the  program,  highly  productive  methods  of  model  construction, 
with  their  correspondingly  long  setup  times  and  scheduling  dif¬ 
ficulties,  were  not  desirable.  Even  mote  importantly,  previous 
experience  indicates  that  eventually  some  recourse  will  have  to 
be  made  to  the  use  of  models  which  will  simulate  the  static  and 


dynamic  response  characteristics  of  the  full-sire  vehicle,  to 
assist  in  the  development  of  theoretical  techniques  and  to  evalu¬ 
ate  the  accuracy  and  reliability  of  these  methods.**  The  full- 


size  vehicles  are  expected  to  be  constructed  of  a  stainless 
steel  or  Rene'  41  mesh  sealed  with  an  elastomer  capable  of  with¬ 
standing  the  1500  °F  temperatures  which  will  ba  experienced 


Development  work  on  these  looms  was  sponsored  by  the  Manufac¬ 
turing  Technology  Laboratory,  Wright -Patterson  Mr  Force  Base, 
Dayton,  Ohio  under  Contract  No.  AF  33(600) -43036. 

s  program  no  effort  was  made  to  fully  simulate  the  static 
amic  characteristics  of  full-size  vehicles.  The  planform 
and  the  relative  outs' de  dimensions  were  chosen  as  repre- 
ntative  of  the  type  currently  being  proposed  for  use  as  re¬ 
entry  vehiclet. 
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during  re-entry.  To  achieve  dynamic  similarity,  the  materials 
used  in  the  model  will  necessarily  have  to  be  different  from 
those  used  for  the  full-scale  vehicle.  In  such  cases,  it  is  con¬ 
sidered  likely  that  simple  rubber  or  impermeable  cloth  material 
will  suffice.  On  the  basis  of  these  considerations,  it  was  de¬ 
cided  to  explore  alternate  methods  of  model  construction  which 
would  require  very  little  equipment  and  permit  a  large  degree  of 
flexibility  insofar  as  model  properties  were  concerned.  The 
following  method  of  model  construction  was  finally  adopted. 

A  styrofoam  mandril,  shaped  to  the  dimensions  of  the  de¬ 
sired  model,  is  attached  to  a  wooden  root  fitting.  A  nylon  (or 
dacron)  cloth,  which  has  been  pre-impregnated  with  latex  rubber 
(General  Latex  Code  No.  IY-10)  and  cut  to  the  correct  size,  is 
wrapped  around  the  styrofoam  core  and  wooden  root  fitting.  All 
seams  are  sealed  with  neoprene  cement.  Next,  the  drop  chords  are 
put  in  by  sewing  back  and  forth  through  the  top  and  bottom  sur¬ 
faces  with  a  needle  and  thread.  To  insure  that  the  drop  chords 
are  reasonably  perpendicular  to  the  model  surface,  it  ha3  been 
found  necessary  to  pre-drill  the  styrofoam  core  using  an  array 
of  heated  needles,  since  the  sewing  needle  has  a  tendency  to  be¬ 
come  misaligned  when  piercing  the  styrofoam.  For  the  models 
used  in  this  program,  a  drop  chord  spacing  of  one-quarter  inch 
has  been  found  to  be  satisfactory.  Also,  it  is  not  necessary  to 
sew  the  leading-  and  trailing -edge  regions  sinc<®  they  are  rounded 
and  will  automatically  assume  the  desired  cylindrical  shape  upon 
inflation.  After  all  the  drop  chords  are  sewn,  a  coat  of  latex 
is  applied  to  seal  the  holes  made  by  the  drop  chords.  The  final 
atep  is  to  remove  the  styrofoam  core  by  dissolving  it  with  a 
common  solvent  such  as  toluene. 

It  is  evident  that  the  foregoing  is  a  relatively  quick  and 
inexpensive  way  to  produce  inflatable  models.  One  of  its  greatest 
advantages,  which  has  not  yet  been  utilized  to  any  great  extent, 
la  that  the  size  or  shape  of  the  model  depends  on  the  styrofoam 
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core  and  can  therefore  be  very  easily  changed.  In  addition,  the 
model  characteristics  can  be  varied  by  (1)  selecting  a  different 
cloth  material  for  the  model  surface,  (2)  selecting  a  different 
latex,  (3)  varying  the  member  of  layers  of  cloth  imbedded  in  the 
latex,  (4)  varying  the  number  of  coats  of  latex,  and  (5)  varying 
the  orientation  of  the  cloth  threads. 

Twenty-one  delta-wing  models  were  constructed  for  use  in 
the  program.  A  summary  of  the  model  configurations  and  pertlnen* 

model  properties  is  given  in  Table  2.1.  > 

£ 

For  the  first  two  models,  the  nylon  cloth  was  simply  ; 
wrapped  around  the  leading  edge  of  the  styrofoam  mandril;  thus, 
the  warp  orientation  (with  respect  to  the  model  root)  was  zero  f 
degrees  for  the  top  surface  and  50  degrees  for  the  bottom  surface 
In  all  the  other  models,  the  surface  material  was  cut  and  a  seams 
made  along  the  leading  edge,  so  that  the  warp  orientation  was  tht 
same  for  both  top  and  bottom  surfaces.  | 

In  addition  to  these  models,  other  models,  including  | 
square,  rectangular,  and  even  a  model  with  a  large  initial  curva* 
ture,  have  been  constructed  by  students  at  MIT  while  working  on  | 
thesis  or  laboratory  projects.  Photographs  of  some  of  the 
are  jhown  in  Fig.  2.1,  which  also  includes  a  picture  of  a  model 
constructed  and  donated  by  the  Goodyear  Aerospace  Corporation, 
Akron,  Ohio. 

A  system  for  maintaining  any  desired  pressure  differentia] 
p^,  between  model  internal  pressure  and  any  reference  pressure 
(atmospheric  pressure  during  static  and  vibration  testing,  tunne! 
static  pressure  during  wind-tunnel  tests)  Is  connected  to  the 
model  through  the  wooden  root  fitting  (Fig.  2.2). 

The  basic  model  instrumentation  consisted  of  four  strain 
gages  (Baldwin  PA-3  post -yield  type)  for  use  in  static  and 
dynamic  measurements.  The  gages  were  fastened  to  the  outer  mode! 
surface  (Fig.  2.2)  using  neoprene  cement  and  then  covered  by  a 
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protective  coating  of  latex.  Each  of  the  strain  gages  was  used 
as  a  single  active  arm  in  a  four-ara  strain-gage  system,  giving 
four  independent  strain-gage  circuits.*  Eecause  of  the  large  de¬ 
flections  encountered  in  these  models,  care  was  taken  to  allow 
enough  slack  in  the  strain-gage  leads  (such  as  zigzagging  the 
leads  as  shown  in  Fig.  2.2)  to  prevent  lead-wire  breakage  during 
model  tests. 

B.  STATIC  TESTS 

The  response  of  Model  1  to  a  static  uniform  load  was  de¬ 
termined  by  placing  sheets  of  lead  on  the  model  surface  and 
measuring  the  deflection  at  three  points  along  the  leading  and 
trailing  edges.  The  lead  sheet  was  cut  into  small  pieces  (2  x 
2  in.)  which  were  joined  by  tape,  so  that  the  lead's  contribution 
to  the  overall  stiffness  of  the  complete  system  (model  and  weight) 
was  negligible.  The  results  of  the  tests  are  given  in  Fig.  2.3. 

In  Table  2.2,  the  experimental  values  are  compared  with  theo¬ 
retical  predictions  based  on  the  assumption  of  pure  shear  defor¬ 
mations  (see  Appendix  I,  Section  A),  for  the  case  when  the  model 
pressure  is  2  psig. 

From  the  results,  several  important  observations  can  be 

made. 

(1)  The  model  deflection  does  not  vary  inversely  with 
model  pressure,  as  it  would  if  the  deformations 
were  due  entirely  to  shear. 

*  (2)  The  range  of  linearity  of  model  deflection  versus 
load  depends  on  the  model  pressure. 

*The  original  models  were  Instrumented  with  only  the  two  aft 
strain  gages  which  were  wired  to  give  two  Independent  strain- 
gage  circuits.  It  was  later  found  necessary  to  add  two  more 
strain  gages  to  overcome  the  problem  of  frequent  tunnel  shut¬ 
downs  for  replacement  of  gages  ruptured  during  tunnel  starts  and 
model  flutter. 
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(3)  The  theoretical  deflections,  wherein  bending-type  de¬ 
formations  are  not  accounted  for,  are  considerably 
lower  than  the  experimental  deflections. 

From  the  foregoing,  it  is  evident  that  the  deflection  is 
composed  of  bending  and  shear  type  deformations,  and  that  both 
are  equally  important.  This  fact  is  also  verified  by  the  vibra¬ 
tion  results  which  are  presented  in  the  following  section. 

C.  VIBRATION  TESTS 

1.  Test  Procedures 

The  excitation  system  used  for  the  vibration  tests  is  a 
low  frequency  speaker.  The  speaker  cone  is  covered  with  a 
tapered  aluminum  cover,  about  0.020-in.  thick,  which  has  a 
1.25-in.  orifice  to  concentrate  the  air  column  at  the  model  sur¬ 
face.  The  vibration  frequencies  and  damping  coefficients  were- 
obtained  from  vibration  decay  records.  Node  lines*  were  deter¬ 
mined  by  sprinkling  No.  40  grit  carborundum  dust  on  the  model 
surface,  while  the  model  was  vibrating  in  a  natural  mode.  This 
system  was  found  to  be  suitable  for  light  vibrating  models.  The 
only  problems  encountered  were  in  connection  with:  (a)  overexci¬ 
tation,  especially  of  the  higher  modes,  and  (b)  difficulty  in 
properly  locating  the  fourth  mode. 

For  large  excitation  forces,  local  deflections  were  found 
to  be  superimposed  on  the  natural  mode  responses.  This  effect 
was  minimized  by  positioning  the  speaker,  for  each  mode,  such 
that  the  largest  amplitude  was  obtained  for  a  fixed  speaker 
power  input.  Changing  the  speaker  position  had  little  effect  on 
the  node  lines  of  most  modes.  However,  in  a  few  isolated  in¬ 
stances  when  weak  and  not  clearly-defined  modes  were  encountered, 
the  node  line  positions  were  strongly  dependent  on  the  speaker 

*More  specifically,  1-g  lines.  In  the  node  line  presentations 
of  Figs.  2.5  through  2.9,  the  shaded  areas  are  the  regions  where 
the  accelerations  are  less  than  lg. 
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position.  In  these  cases,  the  speaker  was  positioned  so  that  a 
sizable  amplitude  and  clearly-defined  node  lines  were  obtained 
at  the  same  time. 

Another  difficulty,  in  varying  degree,  was  experienced 
when  exciting  the  fourth  mode.  The  speaker  position  was  critical 
for  two  reasons:  (a)  it  influenced  the  model  frequency  for  maxi¬ 
mum  peak,  r’.r.i  (b;  if  somewhat  removed  from  its  optimum,  this  mode 
could  be  easily  overlooked  in  the  frequency  tests.  (The  latter 
actually  occurred  in  some  of  the  earlier  tests.)  Despite  the 
care  exercised  in  testing  of  some  models,  it  was  not  possible  to 
excite  the  fourth  mode  to  a  satisfactory  amplitude  level;  and  the 
vibration  traces  showed  a  beat  frequency  which  was  approximately 
equa*.  to  tk:  of  the  first  mode.  The  damping  coefficient  wa3 
difficult  to  measure  in  any  of  these  cases,  as  the  lower  mode 
decay  overshadowed  that  of  the  higher  mode  following  speaker  shut¬ 
off. 

2-  Test  Results  and  Discussion 

Table  2.3  presents  in  chronological  order,  the  measured 
frequencies  and  damping  coefficients  of  all  21  models.  The  list 
includes  the  following  free-air  vibration  tests. 

(a)  Those  conducted  a  few  days  before  the  wind-tunnel- 
tesu  periods  (at  HIT) . 

(b*  ihos»e  conducted  just  prior  to  the  flutter  tests 
(at  wlr.d-tumel  location). 

(c)  All  repeat  tests  made  between  flutter  runs. 

Node  lines,- which  are  typically  representative  of  most  of 
the  models  tested,  are  presented  in  Figs.  2.5  and  2.6,  for  65-  and 
75-degree-sweep  models,  respectively.  These  figures  apply  for 
the  range  of  Internal  pressures  considered,  since  in  the  majority 
of  cases  the  node  line  locations  were  not  altered  appreciably  with 
changes  in  p^.  Separate  node  line  patterns  are  also  presented  in 
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*'  Figs.  2.7  Co  2.9  for  those  models  which  behaved  otherwise.  These 

cases  are  also  noted  in  Table  2.3.  Whenever  accurate  mode  shapes 
were  needed  for  analysis  of  a  given  model,  the  actual  node-line 
t  data  obtained  during  the  pertinent  vibration  tests  were  used 

(see  Appendix  1,  Section  B). 

As  previously  mentioned,  difficulty  was  experienced  in 
exciting  the  fourth  mode.  The  fourth  mode  was  not  defected  at 
all  in  Models  1  and  2.  This  mode  was  first  observed  in  an  early 
test  made  on  Model  3,  which  has  the  warp  direction  parallel  to 
the  root.  At  that  time,  it  was  felt  to  be  the  third  mode  strongly 
distorted  by  speaker  position.  Because  the  mode  could  not  be 
found  in  Model  4  at  the  low  pressures  (equal  or  less  than  4  psi, 
pressures  used  in  the  initial  tests),  it  was  no  longer  looked 
for.  Later,  when  the  last  of  the  dacron  models  were  introduced 

« 

and  tested,  the  fourth  mode,  which  was  vague  or  nonexistent  in 
the  early  nylon  models,  was  rediscovered  and  found  to  be  ouch 
t  stronger.  It  was  traced  and  recorded  for  all  succeeding  tests. 

In  the  data  presented,  the  mode  numbers  were  assigned  on  the 
basis  of  the  observed  node  lines;  the  fact  that  no  data  are  given 
for  Mode  4  should  be  interpreted,  in  most  instances,  to  mean  that 
the  fourth  mode  was  present  and  simply  overlooked  during  the  vi¬ 
bration  tests.  For  some  special  cases,  as  noted  In  the  table, 
repeated  attempts  to  excite  this  fourth  mode  were  unsuccessful, 
and  thus  it  must  be  assumed  that  no  such  mode  exists. 

The  damping  coefficients  g^,-  .  .  .  gj  are  missing  from 
certain  sections  of  Table  2.3,  because  either  the  procedure  used 
was  later  found  to  be  unsatisfactory  or  the  trace  amplitudes  were 
so  snail  that  the  signal-to-noise  ratio  of  the  recorded  strain- 
gage  output  was  too  low  to  give  an  accurate  reading.  In  the  early 
vibration  tests  on  Models  1  through  10,  the  decay  records  were  ob¬ 
tained  following  cutoff  of  speaker  power  at  the  oscillator.  Thus, 
the  8 low  decay  of  the  oscillator  was  incorporated  into  the  model 
decay,  which  resulted  in  an  apparent  reduction  in  g's.  In  all  the 
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remaining  tests,  the  power  was  cut  off  at  the  speaker.  The 
speaker  cone  had  sufficiently  high  damping  characteristic  and 
did  not  influence  the  model  decay. 

The  node  lines  depicted  in  Figs.  2.5  and  2.6  for  the  first 
four  typical  mode  shapes  of  the  inflatable  model  resemble  those 
of  a  cantilevered  delta  wing  with  approximately  the  same  sweep 
and  made  of  a  uniform  thickness  aluminum  plate  (Ref.  6). 

However,  the  frequency  ratios  between  various  modes  are  quite 
different  for  the  two  types  of  models.  There  is  also  some  simi¬ 
larity  between  the  fifth  mode  of  the  inflatable  model  and  the 
sixth  mode  of  the  aluminum  wing. 

For  ...be  single-layer  nylon  models,  the  effect  of  changing 
the  surface  material  warp  direction  from  zero  and  50  degrees  rela¬ 
tive  to  the  root  chord  on  top  and  bottom  surfaces  (as  in  Models  1 
and  2)  to  zero  degrees  (as  in  all  models  of  this  group)  shows  up 
mainly  as  differences  between  the  fifth  mode  shapes  (see  Figs- 
2.5  and  2.6).  The  first  three  modes  and  frequencies  were  es¬ 
sentially  the  same  for  both  orientations,  with  the  minor  varia¬ 
tions  present  attributed  to  the  slight  differences  in  the  masses 
(due  to  different  thicknesses  of  rubber  coatings).  On  the  other 
hand,  the  warp  orientation  bad  a  marked  effect  on  the  mode  fre¬ 
quencies  of  the  double-layer  dacron  model.  With  the  warp  ori¬ 
entations  at  zero  and  45  degrees  for  the  two  plies  (as  in 
Models  9  and  13),  the  frequencies  were  substantially  higher  than 
those  with  both  warp  directions  at  zero  degrees  (as  in  Models 
13-1  and  19).  Model  20,  which  .was  designed  to  duplicate  Model  13, 
did  not  exhibit  this  type  of  behavior.  In  fact  the  vibration  fre¬ 
quencies  of  Model  20  are  very  close  to  those  of  Models  13-1  and 
19.  This  discrepancy  cannot  be  explained  on  the  basis  of  masa 
differences  alone;  Model  20  is  about  21  per  cent  heavier  than 
Model  13,  but  the  frequency  variations  are  of  the  order  of 
30  per  cent. 
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For  the  double-layer  dacron  models,  the  mode  shapes  are 
similar  to  those  of  the  typical  single-layer  models.  Hcwever, 
one  should  note  the  closeness  cf  the  second  and  third  frequencies, 
especially  at  low  internal  pressures.  Specif ically.  Model  9  at 
Pj  *•  1  psig  has  £2  *  44.. 0  and  *  45.1  cps;  Models  13,  13-1,  19, 
and  even  10  (with  75-aegTee  sweep)  exhibited  similar  behavior, 
but  to  a  lesser  degree.  The  closeness  of  modes  is  not  limited 
to  the  second  and  third-  For  Models  19  and  20,  the  frequency 
ratios  between  modes  four  and  five  are  closer  to  unity  than  that 
between  any  other  two  observed  frequencies.  This  trend  towards 
the  higher  frequencies  is  noted  also  for  the  75-degree  sweep, 
single-layer  models;  and  in  particular  between  the  third  and 
fourth  modes  of  Models  15,  16,  and  17,  and  the  ftxirth  and  fifth 
of  Model  18.  For  the  latter  cases,  this  persists  even  for  the 
higher  pressures.  Such  a  proximity  of  two  modes  presents  po¬ 
tential  flutter  problems,  a  point  that  will  be  further  discussed 
in  connection  with  the  flutter  results. 

Variations  of  frequency,  between  different  tests  of  the 
same  model  inflated  to  equal  internal  pressures,  are  sometimes 
greater  than  what  couldi  be  attributed  to  experimental  error  alone. 
Some  -easons  can  be  given  to  explain  such  variations.  The  most 
■.important  are  the  changes  in  material  properties  with  temperature 
1  humidity,  and  to  a  lesser  degree  with  age.  This  point  is 
-urther  discussed  in  the  section  on  material  properties  (see 
Section  2,  Subsection  F).  In  order  to  obtain  on  estimate  of  the 
temperature  effect  on  the  free-vibration  frequencies  of  the  nylon 
models,  Model  6  was  vibrated  in  a  temperature  environment  which 
could  be  varied  from  room  temperature  (70  °F)  to  about  ISO  °F. 

The  results  of  the  test,  for  the  range  of  temperatures  70  °F  £ 

T  £.  130  °F,  are  presented  in  Table  2.4  for  modes  1,  2,  3,  and  5 
and  a  model  Internal  pressure  of  2  psig.  It  may  be  seen  that 
around  room  temperatures,  the  first  frequency  decreases  by  about 
4  to  5  per  cent  for  every  20-degree  rise  in  temperature.  The  same 
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trend  may  be  noted  for  the  higher  modes,  but  to  a  lesser  degree. 

A  second  reason  for  variations  in  frequencies  between  repeated 
tests  is  the  possibility  of  undetectable  damage  to  the  model.  A 
particular  case  in  point  is  Model  1,  where  the  frequencies  of  the 
last  set  are  lower  chan  in  the  initial  set.  And  yet  another 
reason  may  be  f-'und  in  possible  amplitude  effects.  No  attempt 
was  made  to  vibrate  any  model  with  the  same  amplitude  for  each 
test.  The  material  stiffnesses  are  somewhat  nonlinear  with  re¬ 
spect  to  loading.  Whether  this  alters  noticeably  the  frequencies 
is  not  known  at  the  present. 

In  conventional  flutter  work,  the  configurations  are  rela¬ 
tively  heavy,  and  the  apparent  mass  effects  are  so  small  that  they 
are  usually  ignored.  Because  of  the  lightness  of  the  present 
models,  however,  it  is  anticipated  that  there  will  be  a  noticeable 
apparent  mass  effect  on  the  vibration  frequencies  measured  in  air 
when  compared  with  chose  measured  in  vacuum.*  The  extent  of  this 
effect  on  the  fundamental  mode  of  the  present  models  may  be  esti-  . 
mated  simply  as  follows:  Let  M_  denote  the  mass  of  air  contained 
in  the  seaicone  formed  by  rotating  the  delta  model  about  its 
root  chord;  and  let  it  be  distributed  uniformly  as  added  mass  to 

the  total  mass  m  of  the  model.  If  f  and  f  are  the  fundamental 

v  p 

frequencies  in  vacuum  and  at  any  ambient  pressure,  respectively, 
the  rough  approximation 

4  -  I  «  si  -L 

4  "  ■L~  i  n 

holds  under  the  stated  assumption,  recognizing  that  M aAn  is 
still  a  small  quantity.  Applied  to  the  lightest  model  iNo.  2), 
for  Instance,  this  results  in  s  reduction  of  13  per  cent  in  frs- 
quency,  when  fp  is  compared  to  To  confirm  this  experimentally. 

The  frequencies  and  mode  shapes  in  vacuum  are  the  appropriate 
ones  to  use  In  flutter  work. 
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the  model  was  vibrated  in  a  vacuum  chamber  at  several  valuet;  of 
reduced  ambient  pressures  and  at  two  model  pressures  (2  and 
4  psig).  The  experimental  values  are  given  in  Fig.  2.10,  where 
a  straight  line  is  passed  th:\'Ugh  the  points,  since  f  /f  is  ex- 
pected  to  have  linear  variation  with  M_  (or  equivalently  with 
ambient  pressure)  a  cording  to  the  above  equation.  For  Model  2, 
the  predicted  reduction  in  the  fundamental  frequency  due  to  ap¬ 
parent  mass  effects  is  somewhat  higher  than  the  corresponding 
value  of  7.1  per  cent  from  the  observations.  For  the  higher  fre¬ 
quencies,  one  may  anticipate  smaller  influences;  this  point  is 
confirmed  by  the  experimental  results.  The  fact  that  experiment 
shows  about  half  the  predicted  effect  may  be  easily  explained. 

The  assumption  in  the  analysis  would  be  appropriate  if  all  the 
points  of  the  wing  were  vibrating  at  the  same  amplitude.  For  the 
first  cantilever  frequency,  with  parabolic  spanvise  variation  in 
mode  shape,  a  more  realistic  distribution  of  apparent  air  mass 
would  have  included  some  weighting  factors  based  on  the  mode 
shape. 

Some  of  the  experimental  frequencies  are  compared  in 
Table  2.5  with  those  from  the  simplified  calculations  based  on 
the  assumption  that  the  deformations  are  only  of  the  pure  shear 
type  (see  Appendix  I,  Section  A).  The  comparison  indicates  that 
this  sort  of  simplification  is  not  valid.  For  Models  1  and  8t 
the  experimental  frequencies  are  substantially  lower  than .those 
j  edicted,  even  with  due  allowance  for  apparent  mass  effects,  in¬ 
dicating  appreciable  bending  deformations.  The  comparison  be¬ 
comes  worse  when  one  considers  that  the  additional  shear  stiffnesses 
from  the  leading  and  trailing  edges  were  unaccounted  for  in  the 
analysis.  For  Model  9,  the  results  are  reversed  in  that  the  ex¬ 
perimental  values  sre  somewhat  higher  than  the  predicted  ones. 

This  is  mainly  due  to  two  causes.  First,  bending  effects  should 
be  such  lower  in  Model  9  with  its  double  layer  of  dacron  material. 
Second,  Model  9  actually  has  four  layers  of  sloth  material  along 
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the  leading  and  trailing  edges,  due  to  overlap  at  the  sea ms,  and 
the  transverse  shear  rigidity  of  these  edges  adds  a  large  amount 
of  stiffness  to  the  model. 

The  effect  of  model  pressure  differential  on  the  vibration 
frequencies  can  be  determined  from  Table  2.3.  The  frequencies  do 
not  vary  as  the  square  root  of  the  model  pressure  differential, 
as  would  be  expected  if  the  model  stiffness  was  totally  derived 
from  the  internal  pressure. 

D.  VIBRATION  MODE  SHAPES  USING  THE  MIRROR  TECHNIQUE 

1.  Test  Procedures 

An  accurate  knowledge  of  the  vibration  mode  shapes  of  in¬ 
flatable  models  would  be  extremely  useful  for  checking  the  re¬ 
sults  of  theoretical  vibration  calculations  and  also  for  calcu¬ 
lating  the  generalized  mass  and  aerodynamic  forces  needed  in  the 
flutter  analyses.  In  this  manner,  the  mode-shape  data  can  be 
used  in  the  final  steps  of  one  problem  area  (vibration  analysis) 
and  the  initial  steps  of  another  problem  area  (flutter  analysis), 
making  it  possible  to  proceed  on  both  problems  independently  of 
one  another. 

As  an  initial  effort  toward  providing  this  mode-shape 
data,  an  attempt  was  made  to  measure  vibration  mode  shapes  using 
an  optical  technique  which  has  been  successfully  applied  to  mode- 
shape  measurements  on  rectangular-  and  delta-wing  models  (Refs. 

7,  8,  9).  For  the  tests,  optically  flat  mirrors  (3/16. x  3/16  x 
0.02C  inches)  were  glued  to  one  side  of  the  wing  at  each  of  37 
stations,  using  bee's  wax  as  an  adhesive  (Fig.  2.11).  As  shown 
in  Fig.  2.12,  the  wing  was  mounted  vertically,  with  the  root 
parallel  to  the  floor,  and  with  the  mirrors  facing  a  4  x  8-foot 
board  located  about  10  feet  away  and  covered  with  grid  lines 
spaced  at  0.40-inch  intervals.  Light,  from  a  point  Source  placed 
under  the  4  x  8-foot  board,  was  passed  through  37  lenses  (1/4-in. 
x  1-in.  diameter)  mounted  in  a  rigid  board  placed  just  in  front 
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*  of  the  wing.  The  light  rays  were  made  parallel  by  the  lenses 

and  then  reflected  from  the  mirrors  back  through  the  lenses, 
forming  an  image  of  the  point  source  (called  a  dot  from  here  or.) 

‘  on  the  4  x  8-foot  board,  for  each  mirror  on  the  wing. 

When  the  wing  was  vibrated  in  one  of  its  natural  modes, 
the  dots  traced  out  lines  on  the  4  x  8-foot  board  which  gave  a 
direct  measure  of  the  relative  magnitude  and  direction  of  the 
slopes  of  the  wing,  during  each  cycle  of  vibration,  at  each 
mirror  station.  The  slope  in  the  horizontal  (x  or  streamwise) 
and  vertical  (y  or  spanwise)  directions  was  then  determined  by 
measuring  the  appropriate  components  of  the  lines  traced  out  by 
the  moving  dots.  Integration  of  these  slopes  along  any  given 
direction  yielded  the  modal  deflection  amplitude  at  the  various 
stations. 

< 

To  facilitate  data  reduction  and  insure  that  all  the  data 
were  recorded  at  the  same  vibration  amplitude,  a  photograph  was 
»  taken  of  all  the  lines  on  the  board  while  the  model  was  vibrating 

in  a  natural  mode.  No  particular  effort  was  made  to  initially 
place  the  dots  at  a  particular  spot  on  the  board.  The  dots  were 
simply  identified  prior  to  taking  the  photograph  and  an  identi¬ 
fication  number  was  written  beside  each  dot  on  the  board.  This 
procedure  required  a  double  exposure;  the  first  exposure  was 
made  with  room  lights  on  to  record  the  grid  lines  and  identifying 
numbers,  while  the  second  exposure  was  made  in  complete  darkness 
to  record  the  dot  traces.  Figure  2.13  is  typical  of  the  photo¬ 
graphs  obtained;  for  data-reduction  purposes,  the  photographs 
were  enlarged  to  13  x  24  inches'. 

c  2.  Results  and  Discussion 

Mode-shape  data  was  taken  for  the  first  six  modes  of 
Models  7  and  8  at  model  pressure  differentials  of  1,  2,  3,  and 
4  psig.  Only  the  data  for  the  first  three  modes  of  Model  8  at  a 
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model  pressure  of  2  psig  have  been  reduced  to  provide  data  for 
the  two-  and  three-degrees-of-freedoo  flutter  calculations  dis¬ 
cussed  in  Appendix  I.  The  mode  shapes  are  drawn  in  Figs.  2.14 
through  2.16. 

As  previously  mentioned,  the  data  was  reduced  ly  measuring 
the  streamwise  and  spanwise  components  of  the  trace  for  each 
mirror.  The  data  was  plotted  in  both  directions  and  along  di¬ 
agonals  as  well,  and  the  curves  were  used  to  find  values  of  the 
slope  at  the  root  and  at  intermediate  stations  where  there  were 
no  mirrors.  It  should  be  noted  that  the  algebraic  signs  of  the 
slope  cannot  be  obtained  from  the  photographs,  but  oust  be  de¬ 
termined  by  recourse  to  node  lines  obtained  from  previous  vibra¬ 
tion  tests  using  the  techniques  discussed  in  Section  2C;  this 
presents  a  difficult  problem  for  the  higher  modes.  The  vibration 
amplitude  was  then  fctmd  by  numerical  integration  of  the  curves 
obtained  in  the  previous  step.  To  find  the  vibration  amplitude 
at  any  point,  there  are  two  main  paths  along  which  one  may  pro¬ 
ceed.  Taking  point  22  In  Fig.  2.11,  for  example.  It* is  possible 
to  find  its  acqplitude  by  (!)  integrating  along  the  spanwise  di¬ 
rection  starting  from  point  20R,  or  (2)  integrating  along  another 
spanwise  station,  say  the  rearmost  station,  from  point  1R  to 
point  3  and  from  there  integrating  in  a  streamwise  direction  to 
point  22.  In  practice  it  has  been  found  necessary  to  carry  out 
such  a  dual  process  for  each  point,  compare  the  results,  and 
then  reiterate  on  the  original  data  to  make  the  overall  agreement 
better.  From  the  type  of  agreement  obtained  for  the  Inflatable 
models  by  this  method.  It  Is  recognized  that  the  results  shown  in 
Figs.  2.14  through  2.16  are  not  nearly  as  accurate  as  previous 
results  obtained  by  applying  essentially  the  same  techniques  to 
solid  wing  models  (Refs.  7,  8).  The  main  source  of  error  lies  In 
determining  the  spanwise  slope  at  the  model  root.  This  slope  la 
zero  for  solid  cantilever  models,  but  such  is  not  the  case  for 
inflatable  cantilever  models,  as  indicated  by  Figs.  2.14  through 
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•  2.16.  In  fact,  the  results  showed  that  there  are  very  marked 
fluctuations  in  the  spanviss  slope  near  the  root,  and  that  the 
mirrors  should  have  been  placed  much  closer  to  the  root  than  the 

*  one-inch  spanwise  spacing  shown  in  Fig.  7.11. 

Another  problem  encountered  during  the  tests  was  that  the 
mirrors  were  found  to  bend  as  the  model  pressure  was  increased. 
This  was  caused  by  local  stretching  of  the  model  surface  between 
drop  chords  and  could  be  alleviated  by  using  slightly  stronger 
(and  heavier)  mirrors.  The  bending  of  a  mirror  manifested  itself 
by  turning  the  dot  on  the  board  into  a  line.  In  an  extreme  case, 
the  dot,  which  is  normally  about  3/32  inches  in  diameter,  was 
expanded  into  a  3/16  x  3/4-inch  line.  When  this  occurred,  it 
was  necessary  to  remount  the  mirrors.  The  effect  was  found  to 
,  be  minimized  by  installing  the  mirrors  at  the  higher  model  pres¬ 

sures,  the  bending  of  the  mirrors  being  less  noticeable  with  a 
decrease  in  model  pressure. 

i  With  regard  to  the  effect  of  the  mirrors  on  the  vibration 

results,  the  experimental  frequencies  and  node  lines  were  de¬ 
termined  at  a  model  pressure  of  2  psig,  both  with  and  without 
mirrors  mounted  on  the  model,  using  the  vibration  techniques 
described  previously.  There  was  no  measurable  difference  between 
the  results  for  both  cases. 

As  a  check  on  the  orthogorality  of  the  measured  modes,  the 
first  three  mode  shapes  were  useo  to  compute  the  generalized 
masses  (relative  values): 

-  1197  Hu  -  34.4 

H 22  -  582.4  -  39.0 

1*33  -  141.3  -  18.3 

Since  the  cross  terms  are  not  zero,  the  modes  are  not  normal  to 
each  other  and  ire  therefore  not  the  true  vibration  modes.  An 
estimate  of  the  error  can  be  determined  by  comparing  Mu  with  the 
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square  root  of  x  .  In  the  worst  case  (M^  compared  with 

v/’  HI1M33)  '  amounts  to  a  little  less  chan  ten  per  cent.  It 

should  be  noted  that  the  generalized  masses  presented  above  have 
been  corrected  for  the  effect  of  a  nonuniform  mass  distribution 
caused  by  material  lap-over  at  the  leading-  and  trailing-edge 
seams.  In  some  cases  this  correction,  which  was  only  approximate 
in  the  sense  that  the  exact  amount  of  lap-over  was  not  known,  was 
of  the  same  order  of  magnitude  as  the  cross  terms  and  could  ac¬ 
count  for  an  appreciable  part  of  the  error. 

E.  INFLUENCE  COEFFICIENTS 

1.  General  , 

Although  the  techniques  for  measuring  influence  coefficients 
and  calculating  the  natural  vibration  modes  and  frequencies  are 
well  known  (Ref.  10),  several  interesting  observations  pertinent 
to  inflatable- type  models  are  noteworthy. 

In  the  measurement  of  influence  coefficients  on  solid 
models,  it  is  customary  to  use  a  "point  force"  in  order  to  ac¬ 
curately  locate  the  point  of  application.  This  point  load  es-’ 
sentially  acta  as  a  distributed  force,  with  the  centroid  of  force 
acting  through  the  point  of  application.  In  measurements  on 
actual  aircraft  it  is  sometimes  necessary  to  use  re-enforcing 
pads  at  the  loading  points  to  prevent  local  buckling.  For  in¬ 
flatable  models,  locally  applied  load  results  in  a  large  local 
deformation  which  is  not  at  all  representative  of  the  type  of 
loads  and  deformations  that  will  normally  be  experienced  by  the 
model.  It  was  therefore  decided  to  apply  the  loads  to  the  model 
with  lead  pads  placed  on  the  top  surface  of  the  model.  The  pads 
had  small  holes  in  the  centers  to  enable  relatively  accurate 
placement  over  the  desired  loading  point,  and  additional  weights 
could  be  added  on  top  of  the  pads  as  required.  The  effects  of 
the  loading-pad  diameter  were  investigated  using  three  lead  pads 
with  diameters  of  1,  2,  and  4  inches.  Three  representative  in- 
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fluence  coefficients  were  measured  using  the  different  pads  and 
the  results  are  given  in  Table  2.6.  It  is  evident  that  the  local 
deflection  (C55)  increases  as  the  load  approaches  a  point  appli¬ 
cation,  while  pad  diameters  between  2  inches  and  4  inches  give 
essentially  the  same  results.  Therefore,  all  subsequent  tests 
were  carried  out  using  the  2 -inch  pads. 

In  influence-coefficient  measurements  it  is  important  that 
the  transducers  used  in  measuring  deflections  do  not  introduce 
any  extraneous  forces  and  restraints  into  the  system.  Because  i 
of  the  large  deformations  experienced  with  inflatable  models, 
under  even  low  loadings,  it  was  possible  to  use  glass  vernier  ^ 
scales  hung  from  the  bottom  surface  of  the  model.  The  movement  • 
of  the  scales  was  observed  through  a  surveyor's  transit.  With 

if 

this  system,  which  satisfies  the  criterion  of  net  adding  ex-  f 
traneous  forces  and  restraints  to  the  system,  it  was  possible  to  t 
estimate  deflections  to  within  +  0.001  inch.  |' 

In  order  to  keep  the  deflection  resolution  error  to  within 
2  per  cent,  the  applied  loads  were  chosen  to  give  a  deflection  of 
over  0.100  inch  at  the  measuring  station.*  Where  this  was  not  tp 
possible,  a  large  number  of  repeated  measurements  were  made  and  Ig 
an  average  deflection  was  determined.  Gl 

Another  problem  encountered  was  that  of  hysteresis,  which  B 
was  minimized  by  repeating  the  load-no-load  cycle  until  a  numbers 
of  repeatable  deflection  measurements  were  obtained.  In  general JS 
this  required  5  to  6  cycles.  V 

In  the  iteration  procedures  used  to  determine  normal  modes! 
and  frequencies  analytically,  it  is  customary  to  impose  a  converJI 
gence  criterion  based  on  successive  eigenvalues.  In  the  present  11 
example,  a  more  restrictive  criterion  had  to  be  adopted  in  order  11 

*Each  measurement  Involves  two  readings,  load  on  and  load  off,  H 
each  with  a  possible  errur  of  0.001  inch.  Sg 
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to  obtain  a  satisfactory  third  node  shape.  The  convergence 
criterion  used  fcr  the  first  and  second  node  shapes  required 
that  no  componenc  of  the  eigenvector  eotrputed  during  successive 
iterations  differed  by  nore  than  1  x  10  ^ \  this  was  with  refer¬ 
ence  to  a  normalized  mcdc  shape  with  maximum  cocponent  equal  to 
one.  For  the  third  mode,  the  convergence  criterion  required  that 
the  difference  in  any  eigenvector  component  be  less  than  1  x  10  ; 

again  with  maxima  component  equal  to  one.  The  use  of  the  con¬ 
vergence  criterion  on  the  eigenvectors,  rather  than  on  the  eigen¬ 
values,  resulted  in  marked  changes  in  some  of  the  components  of 
the  third  mode  shape,  with  differences  as  high  as  30  per  cent  in 
some  cases;  the  final  eigenvalues  were  the  same,  irrespective  of 
the  convergence  criterion  used. 

2.  Results  and  Discussion 

The  first  three  mode  shapes  of  Model  5  were  calculated 
for  a  model  pressure  differential  of  2  psig,  using  measured 
flexibility  influence  coefficients;  influence  coefficient  stations 
and  mass  data. are  given  in  Fig.  2.17.  The  influence  coefficient 
matrix,  which  has  been  made  symmetrical  by  averaging  reciprocal 
elements  and  Cj^,  to  correspond  with  Maxwell's  reciprocity 
relation,  is  given  in  Table  2.7.  The  calculated  mode  shapes  for 
the  first  three  modes  are  given  in  Table  2.8.  The  calculated 
and  experimental  frequencies  and  node  lines  are  compared  in 
Table  2.9  and  Fig.  2.18,  respectively.  In  the  figure,  the  ex¬ 
perimental  nod  a  lines  at  the  root  are  indicated  as  extending 
over  an  area,  since  it  is  not  possible  to  locate  them  exactly 
within  this  region  by  the  1-g  acceleration  technique. 

The  comparison  between  theory  and  experiment  is  reasonable, 
and  It  appears  that  carefully  measured  influence  coefficients  can 
be  used  to  determine  the  natural  modes  and  frequencies  of  in¬ 
flatable  wings. 
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*  F.  MATERIAL  PROPERTIES 

1-  Ceneral 

t  Comparison  between  experimental  model  vibration  results 

and  theoretical  calculations,  based  on  an  assumption  of  pure 
shear-type  deformations  only,  suggests  that  the  bending  deforma¬ 
tions  may  be  of  the  same  order  of  magnitude  as  the  shear  de¬ 
formations.  The  bending  characteristics  are  dependent  upon  the 
properties  of  the  model  surface  material.  Accordingly,  a 
limited  amount  of  testing  was  performed  to  determine  the  struc¬ 
tural  properties  of  the  nylon  cloth  used  in  the  construction  of 
the  models. 

Three  test  procedures  were  used  to  measure  the  structural 
properties:  (1)  uniaxial,  (2)  biaxial  on  cruciform-shaped  samples, 
and  (3)  biaxial  on  cylinders.  The  uniaxial  and  biaxial  cruciform 
tests  suffered  from  size  limitations  and  severe  creep  effects. 
Their  results  are  inconclusive  and  thus  are  not  reported  here. 

The  biaxial  tests  on  the  cylinder  are  described  below  in  detail. 

2.  Description  of  Models 

Two  sets  of  rubberized  nylon*  cylinders  were  constructed 
and  tested.  The  initial  set,  used  early  in  the  investigation, 
consisted  of  two  4-inch  diameter,  17-inch  length  cylinders.  The 
cloth  warp  was  oriented  in  the  axial  direction  for  one  of  the 
cylinders,  and  in  the  circumferential  direction  for  the  other. 

The  other  set  consisted  of  three  4-inch  diameter,  20-inch  length 
cylinders.  These  cylinders  differed  from  one  another,  in  that 
the  warp  directions  were  axial,  circumferential  and  at  45  degrees 
to  the  cylinder  axis,  respectively.  Each  cylinder  was  sealed  at 
both  ends  with  4-inch  diameter  wooden  disks,  and  was  inflated 
to  a  pressure,  pQ,  which  made  the  hoop  (circumferential)  stress 
* 

This  was  the  same  material  as  was  used  for  Models  1  through  8 
inclusive,  having  98  x  98  thread  counts  per  inch.  All  cylinders 
were  coated  with  two  layers  of  latex. 
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equal  Co  Che  stress  in  the  cover  of  the  wing  model  under  a  pres¬ 
sure  differential,  pd-  Testing  at  the  same  stress  levels  in  the 
cylinders  and  in  the  models  is  retired  since  material  properties 
are  known  to  be  dependent  on  the  stress  conditions.  The  cylinders 
were  inflated  to  various  pressures,  with  the  objective  of  cover¬ 
ing  the  contemplated  p^-range  of  the  models. 

To  find  the  required  pressure  in  a  cylinder  as  a  function 
of  the  model  pressure,  consider  the  typical  inflated  delta  mode. 
(Fig.  2.1).  Except  for  parts  near  the  edges,  the  stress  level, 
crs  ,  in  the  covers  is  essentially  equal  in  all  directions,  and 
is  given  approximately  by 


where  h,  t  are  the  model  depth  and  the  cover  thickness,  respec¬ 
tively.  On  the  other  hand,  the  hoop  stress  in  the  cylinder,*^  * 
is  given  by 


where  D  is  the  diameter  of  the  cylinder.  The  thickness,  t, 
is  the  same  for  cylinder  and  model,  as  the  same  coated  material 
was  used  for  both.  The  requirement  that  the  model  and  the 
cylinder  have  the  same  stress  conditions  yields 

b  s  -b 
%  ])  U 

Since  the  axial  stress  in  the  cylinder  is  equal  to  one  half  of  the 
hoop  stress,  an  additional  axial  load  must  be  added  to  increase 
the  axial  stress  level  to  that  of  the  hoop  stress.  This  additional 
load  turns  out  to  be 


(Load) 
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The  tes  :s  were  run  with  the  icr-gitudiiv'l  axis  of  the  cylinder 
in  a  vertical  direction,  so  that  a  weight,  W  ,  could  be  applied 
by  simply  hanging  it  from  the  bottom  face  of  the  cylinder.  In 
this  position,  the  wcignt  of  the  cylinder  end  plates  contributes 
to  the  load.  Thus,  the  added  weight,  W  ,  was  taken  to  be  the 
required  load  minus  one  half  the  total  weight  of  the  cylinder 
unit.* 

The  tests  on  the  cylinders  may  be  split  into  two  phases: 
the  static  and  the  dynamic  series.  The  dynamic  tests  were  ne^ed 
in  view  of  the  sizable  creer  observed  during  the  static  tests. 

In  each  case,  the  material  stiffnesses,  Et  and  Gt  ,  were 
measured  for  each  cylinder  and  for  various  levels  of  internal 
pressure.  Static  tests  were  also  performed  to  measure  the 
Poisson's  ratios,  fX.  .  No  attempt  was  made  to  measure  p,  dy- 
nan.ically,  for  this  presents  a  mere  difficult  and,  as  yet,  un¬ 
resolved  problem. 

3.  Static  Tests 

In  the  static  tests,  each  cylinder  wa3  deformed  in  a  par¬ 
ticular  way  with  application  of  a  known  stress,  and  the  deforma¬ 
tion  was  measured  to  find  the  corresponding  strain.  For  the  de¬ 
termination  of  the  effective  moduli  E„5  ,  extensional 

loads  were  applied.  Here  the  subscripts  (■  ,  yr  ,  and  45  denote 
the  fill,  the  warp  and  the  45  degrees  to  the  warp  directions,  re¬ 
spectively.  The  last  of  these  is  referred  to  as  the  45-degree 
direction  hereafter.  For  the  shear  moduli  ,  and  G  , 

torques  were  imposed  at  the  base  of  the  cylinders.  Q.  and  G 
denote  the  values  obtained  from  the  cylinders  with  axes  parallel 
to  the  fill  and  the  warp,  respectively.  Since  the  warp  and  fill 
are  at  90  degrees,  the  values  of  and  should  he  the  same, 

theoretically.  G*s  is  that  value  of  Q  associated  with  shsar 
imposed  on  a  cylinder  with  its  axis  at  45  degrees  to  the  warp. 

The  major  contribution  to  the  cylinder  mass  comes  from  the 
cylinder  end  plates. 


fiyfW qmwwww m wmmmm y « j  ■■»■■  j  ■  ■  "  !  U  I- 

«+*?#%&  w »W»*  •  ■■•«  «-  ■«'•'•■*»■'• 

\-.-V~,  •„  -;v  vV-w,.!*  ••-  -  .^.-;».^M';-'vy,.  1  -  . ...B^R, ?£■'**'••<* >U;*ij, j^***-1 


MWI^IWWjyW 


<**&« 


$ 


FDI.  TDR  64-147 


The  -test  for  finding  Et  consisted  of  setting  up  the 
equilibrium  conditions  in  the  cylinder  corresponding  to  some 
value  ot  model  pressure  and  of  recording  the  distance,  t  ,  be¬ 
tween  two  reference  points  on  the  cylinder,  one  point  near  each 
end.  The  applied  load  was  then  varied  through  a  range  which 
bracketed  the  original  load  by  removing  or  adding  weights.  The 
elongated  distance  between  the  reference  points,  /  +  aI  ,  was 
recorded  after  each  loading.  The  cylinder  was  allowed  ample 
time  between  loadings  to  return  to  its  preloaded  position.  All 
readings  were  taken  as  quickly  as  possible  to  minimize  creep  ef¬ 
fects.  In  terms  of  the  strain,  Al/ 1  ,  the  exteneicnal  stiff¬ 
ness  is  given  by 


Li 


aP 

ncD(Ai/£) 


where  AP  ic  the  added  or  subtracted  weight. 

To  measure  Qt  ,  equilibrium  conditions  (median  state) 
were  again  established  as  before.  The  cylinder  was  subjected  to 
a  torque  at  the  bottom  disk,  and  the  resultant  rotation  was 
measured  with  a  double-ended  pointer  attached  to  the  bottom  disk. 
Angular  deflections  were  read  from  a  protractor  mounted  on  a 
shelf  below  the  rotating  cylinder.  As  was  done  with  the  ex¬ 
tensions!  tests,  readings  were  taken  as  quickly  as  possible  after 
load  applications,  and  after  each  reading  the  torque  was  removed 
to  allow  the  cylinder  to  return  to  its  zero  position. 

The  contributions  to  the  effective  (or  total)  shear  stiff¬ 
ness,  Qt  ,•  are  primarily  from  two  sources:  (1)  from  the 
coated  nylon,  and  (2)  from  the  internal  pressure.  It  can  be 
shown  that  the  relation  between  the  effective  shear  stiffness, 
and  the  measured  deflection  angle,  Q  ,  is  approxisutely 

‘The  parameter  (  Gt  ),  corresponding  to  A«3  of  Ref.  11,  re¬ 
quire!  for  vibration  calculations  which  include  bending  deforma¬ 
tions  . 
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where  T  is  the  applied  torque  and  i,  is  the  length  of  the 
cylinder  (see  sketch  below).  Under  ideal  conditions  of  no  fric¬ 
tion  between  junctions  of  fill  and  warp  filaments,  etc . 

the  effect  of  the  internal  pressure  can  be  estimated.  Consider 
the  case  of  the  cylinder  with  the  warp  parallel  to  the  axis. 

When  a  torque  is  applied,  the  warp  filaments,  which  are  assumed 
to  act  as  the  vertical  strands  of  a  bifilar  pendulum,  impose  a 
restoring  moment  of  magnitude 


where  is  the  total  suspended  weight,  and  the  (factor)  cor¬ 

responds  to  that  fraction  of  the  total  load  carried  by  the  ver¬ 
tical  strands.  Since  the  rubber  matrix  supports  only  a  small 
percentage  of  the  vertical  load,  the  (factor)  was  taken  to  be 
unity.  This  restoring  torque  resists  a  part  of  the  applied 
torque  V  .  Thus  the  contribution  to  Gt  from  the  material 
alone,  in  the  absence  of  initial  stresses,  would  be  approximate! 
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since  \JT  has  to  equal  for  equal  hoop  and  axial  stresses. 

If  the  effect  of  internal  pressure  is  similar  to  that  of  a  bi- 
filar  pendulum  a9  described  above,  then  (Gt  should  be  in¬ 

dependent  of  pfi.  As  will  be  seen  later  from  the  experimental 
results,  the  dependence  of  (  Gt"  on  is  considerably 

less  than  that  for  Gt-  Had  (  Gt  )m^  exhibited  independence 
on  £  ,  it  would  have  been  possible  to  test  these  types  of  ma¬ 
terials  at  one  pressure,  and  calculate  the  effective  Gt  at  any 
other  pressure  according  to 

Gi  ?._t  1  k 

The  magnitude  of  the  term  Q  fe  ,  obtained  by  assuming  the  bi- 
filar  pendulum  effect  is  comparable  to  (  ,  when  the 

warp  or  fill  a^e  along  the  axis  of  the  cylinder.  For  the  tests, 
with  the  warp  filaments  at  45  degrees  to  the  cylinder  axis,  this 
sort  of  estimation  of  the  effects  of  internal  pressure  is  not 
necessarily  valid.  However,  as  will  be  seen  later,  Gt¥S "t  • 
is  considerably  larger  than  Gf-t  or  ,  and  exhibits  smaller 

percentage  dependence  on  the  internal  pressure  than  either  QJr 
or  G*t  • 

Several  static  methods  were  attempted  to  measure  Poisson's 
ratio  jjl  .  These  were: 

(1)  An  optical  method,  using  a  transit  and  a  scale 

(2)  The  usual  technique  of  strain  gages  mounted 
lengthwise  and  circumferentially 

(3)  A  "null-gage"  method,  which  is  described  below. 
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Method  (1)  was  abandoned  when  the  extension  in  the  circum¬ 
ferential  direction  was  found  to  be  too  snail  to  be  measured  with 
any  degree  of  accuracy.  If  touch  larger  diameter  cylinders 
(  D  ■  2  ft  or  more)  had  been  used,  useful  data  might  have  been 
obtained,  provided  the  cylinders  were  made  correspondingly  longer 
to  minimize  end-plate  effects.  Mathod  (2)  was  also  abandoned 
since  no  reliable  method  could  be  found  for  calibrating  the  gages, 
particularly  in  the  circumferential  direction  where  the  extensions 
were  very  small.  Method  (3)  afforded  seme  promise  for  determining 
yt4  .  Let  (Tand  £  ,  denote  the  stress  and  the  strain,  respectively 


and  further,  let  subscripts  H  and  A  designate  the  hoop  and 
axial  directions,  respectively.  Assuming  linearity,  one  of  the 
stress-strain  equations  is 


Ea 


SSL 


where  is  the  Poisson's  ratio  relating  strain  in  the  hoop 

direction  due  to  stress  in  the  axi:/.’  direction.  If  one  sets 
€h  -  0,  then 


Pah  ~ 


‘‘h  I 

rA  e* 


If  <£*  is  changed  after  balancing  the  strain  gage  at  the  refer¬ 
ence  pressure  and  loading,  oust  also  be  changed  to  keep 

-  0.  This  was  done  as  follows:  the  strain  gage  was  first 
nulled  at  the  reference  pressure  and  loading  condition  (median 
position).  The  axial  load  was  altered  by  adding  or  subtracting 
known  weights.  This  unbalanced  the  strain  gage  in  the  hoop  di¬ 
rection.  The  strain  gage  was  then  nulled  again  by  increasing 
ot  decreasing  the  internal  pressure  as  needed.  Assuming  linear¬ 
ity  about  the  median  state,  it  is  easily  deduced  the  above 

*At  the  median  position.  O'  ■  ^  *0,  i.e.,  0“ and  €  represent 
the  incremental  stress  and  rtrain. 
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stress-strain  relation  that 


ft 


AH 


'  H 


(* 


/  ) 

7TJ)Z  A&7 


where  Afe  is  the  change  in  internal  pressure  required  for  null¬ 


ing,  and  AW  is  the  added  weight  to  alter  the  axial  stress. 
The  ratio  ^E^/f^'has  to  be  taken  from  separate  tests. 


4.  Dynamic  Tests 

The  dynamic  method  for  measuring  Et  of  the  inflated  cyl¬ 
inders  consisted  of  exciting  the  fundamental  axial  vibration  mode 
of  the  cylinder  while  it  was  in  its  desired  biaxial  stress  con¬ 
dition,  and  recording  the  vibration  trace  using  the  output  of  an 
axially-aligned  strain  gage.  The  axial  mode  was  excited  by  sup¬ 
porting  a  small  weight  (about  1  pound)  with  an  electromagnet 
attached  to  the  original  equilibrium  vertical  load.  When  the 
magnet  was  switched  off,  the  weight  fell,  exciting  the  longi¬ 
tudinal  mode.  Knowing  the  frequency  from  analysis  of  the  vibra¬ 
tion  records,  it  is  a  simple  matter  to  obtain  the  stiffness,  Et  , 
from  the  relation 

Ft  a, 


where  m  is  the  total  suspended  mass  (including  one  third  the 
material  mass,  which  corresponds  to  the  well-known  approximation 
for  the  effective  mass  of  a  spring,  and  the  mass  of  one  end  plate), 
and  y  is  the  frequency  in  cps. 

Gt  was  measured  dynamically  by  using  each  cylinder  as 
the  spring  of  a  torsional  pendulum.  The  cylinder  had  one  of  its 
ends  clamped  to  a  horizontal  strut  from  which  the  cylinder  hung 
vertically.  A  three- foot  bar  was  fixed  horizontally  to  the  cyl¬ 
inder's  lower  free  end.  Equal  known  masses  were  placed  near  the 
ends  of  the  bar,  equidistant  from  the  center  of  the  cylinder.  A 
sketch  of  the  setup  is  given  below.  The  added  mass,  required  to 
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equalize  the  hoop  and  axial  stresses,  was  attached  to  the  bottom 
of  the  cylinder  from  the  end  of  a  long  wire.  This  allowed  the 


cylinder  to  twist  independently  cf  the  added  mass.  After  the 
required  weights  were  added  to  the  bottom  of  the  cylinder,  the 
cylinder  was  inflated  to  the  required  pressure,  i  The  adjustable 
(sliding)  masses  (each  of  ffl  slugs)  were  positioned  at  a  known 
distance,  dt  ,  from  the  center  of  the  rod.  The  natural  torsion 
frequency,  J-  (cps),  of  the  cylinder-weight  combination  was  re¬ 
corded  using  the  output  of  a  strain  gage  mounted  on  the  cylinder 
at  an  angle  of  45  degrees  to  the  axial  direction.  The  masses 
were  moved  to  new  positions  on  the  rod,  at  a  distance  from 
the  center,  thus  changing  the  mass  moment  of  Inertia  by  an 
amount 

4J  2 m.  (d 

The  natural  frequency,  f  ,  corresponding  to  the  new 
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mass  positions  was  recorded  as  before,  and  the  effective  torsional 
"stiffness”  was  computed  according  to 


Gt 


3zrLtn  Zf* 


5.  Results  and  Discussion 

The  experimental  results  from  the  cylinder  tests  are  sum¬ 
marized  in  Tables  2.10  through  2.14.  They  include  three  series 
of  tests  which  were  conducted  at  different  times.  In  the  first 

4 

series  the  initial  17-inch  long  cylinder  models  were  used;  In 
the  last  two  series  these  models  were  replaced  by  the  20-inch 
cylinder  models.  Within  each  series,  repeated  tests  were  tun  in 
the  majority  of  cases,  and  the  indicated  values  are  the  computed 
averages . 

In  Table  2.10,  the  extensional  stiffnesses,  Et  ,  along 
the  directions  of  warp,  fill,  and  45  degrees,  are  presented  for 
cylinder  Internal  pressures  corresponding  to  model  pressure  of 
•  1.0  -  6.0  lbs/in.*  in  1-lb/in.2  intervals. 

The  dynamic  and  static  test  results,  generally,  agree 
within  the  experimental  accuracy. 

There  exist  substantial  differences  between  dynamic  tests 
2  and  3;  in  each  series,  the  repeatability  was  within  a  such 
narrower  band  than  this  comparison  indicates.  Since  series  2* 
was  conducted  on  hot  days  in  the  summer  and  series  3  on  cool  days 
in  late  fall  (four  months  apart),  it  was  originally  thought  that 
aging  of  the  rubber  was  responsible  for  this  change,  with  possibly 
some  temperature  and  humidity  effects.  Additional  tests  revealed,, 
however,  that  both  temperature  and  humidity  are  strong  influencing 
factors.  Table  2.11  suoaurlzes  some  repeat  dynamic  tests  for 
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The  presented  values  are  for  internal  pressure  ji  «  0.688  psig 
(corresponding  to  ■  2  psig).*  To  assess  the  influence  of 
humidity,  the  surface  was  wetted  with  water  at  room  temperature 
using  a  sponge.  The  effect  of  temperature  change  was  assessed 
by  heating  the  model  with  a  100-watt  light  bulb  and  reflector 
placed  about  3  inches  from  the  model.  The  results  indicate  sub¬ 
stantial  reductions  in  with  higher  temperature  and  humidity. 

The  significant  variations  of  model  frequencies  with  temperature 
as  noted  in  Table  2.4  further  substantiate  this  point. 

The  values  of  the  extcnsional  stiffnesses  along  the  warp, 
t  are  higher  than  the  corresponding  quantities  along  the 
fill,  EjJ:  •  As  the  mean  stress  levels  are  raised  (i.e.,  with 
larger  internal  pressures),  both  and  E^t  increase,  but  the 
ratio  E„t / E*t  decreases.  This  can  be  explained  by  postulating 
that,  at  low  stress  levels,  the  fabric  stretches  by  thread 
straightening  (also  called  "crimp  interchange"),  whereas  at 
higher  pressures,  thread  stretching  i3  the  dominant  mechanism 
with  a  resultant  increase  in  stiffness. 

E^t  is  considerably  less  than  either  E^E  or  E^ t  ; 
this  is  as  one  might  expect,  for  the  contribution  of  the  nylon 
cloth  in  this  direction  is  considerably  smaller  than  in  the  fila¬ 
ment  directions. 

Parallel  test  results  for  the  effective  shear  stiffnesses 
G^t  ,  Gf-t  »  and  are  given  in  Table  2.12,  for  essentially 

the  same  cylinder  pressures  as  before.  For  the  test  series  2, 
the  following  procedure  was  used.  One  torsional  frequency  was 
observed  for  each  of  the  two  positions  of  mass  m. ,  <£ ,  ,  and  . 
Let  the  value  of  Gt  obtained  from  these  two  frequencies  (  ) 

be  denoted  by  (  G  t  )u  .  The  test  was  repeated  for  a  new  mass 
position  c/j  ,  and  a  third  frequency,  ,  was  recorded.  In  con- 

Although  not  reported  here,  a  few  additional  test  points  showed 
similar  behavior  at  higher  internal  pressures. 
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junction  with  the  data  for  oosition  cLz  ,  another  value  (Gt  )t3 
was  calculated.  Ideally,  identical  results  should  be  obtained 
for  (6^  ),  and  (£  £  ),3  •  This  was  not  the  situation  however, 

especially  in  cases  involving  small  changes  in  frequencies  be¬ 
tween  £  and  £  .  The  calculated  result  for  &t  ,  being  inversely 
proportional  to  the  difference  of  two  observed  periods  of  oscil¬ 
lation,  requires  good  accuracy  in  the  frequency  measurements.  In 
cases  where  the  vibration  traces  were  clear  and  of  sufficient 
length,  reasonable  agreement  was  obtained  between  (Cft'  )/JL  and 
(  (=j-b  )fJ  .  The  entries  given  in  the  table  for  these  series  of 
tests  are  the  averages  of  ( Gt  )  and  (<S£  )zJ  .  For  the  test 
series  3  results,  the  presented  values  are  for  »  with 

positions  cL,  and  as  the  extreme  locations  of  mass  m  .  For 
test  series  1,  2,  and  for  Gjt  and  ,  two  values  are  entered 
in  each  box,  with  the  top  value  for  the  effective  stiffness,  and 
the  lower  for  .  The  latter,  which  is  due  solely  to  the 

material,  is  computed  by  subtracting  from  the  effective  value  an 
amount  due  to  the  bifilar  effect,  as  discussed  previously. 


The  dynamic  and  static  test  results  agree  within  experi¬ 
mental  accuracy. 

» 

Like  the  extensional  stiffnesses,  the  shear  stiffnesses 
are  temperature  dependent  as  evidenced  by  differences  in  tha  dy¬ 
namic  tests  2  and  3,  but  this  dependency  appears  to  a  lesser 
degree. 


The  experiments  confirm  the  expected  results  that 

Mr 

and  that  both  are  raised  with  increasing  pressure. 


The  bifilar  effect  does  not  account  for  the  total  increase 
in  Gt  with  increasing  pressure.  The  additional  and  unaccounted 
Increase  may  be  possibly  explained  by  postulating  that  the  fric¬ 
tion  at  each  crossing  point  of  the  threads  is  larger  at  the  higher 
pressures  and  the  threads  cannot  move  as  freely  relative  to  oae 
another. 
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(G¥S t  shows  a  smaller  percentage  dependence  on  internal 
pressure  and  is  an  order  of  magnitude  larger  than  either  Gut  or 

Gift  . 

The  static  test  results  for  the  Poisson's  ratios  , 

and  JU-¥S  ,  are  presented  in  Table  2.13.  Large  differ¬ 
ences  exist  between  series  2  and  3  values,  with  the  most  marked  , 

ones  at  the  lower  pressures.  Since  the  latex  is  a  major  contribu-  i 
ting  factor  for  and  ,  it  is  not  surprising  to  see  a 

large  influence  of  temperature  and  humidity. 


and  decrease  rapidly  with  increasing  pressure. 

At  low  pressures,  Poisson's  ratios  are  large  since  crimp  inter¬ 
change  dominates.  At  higher  pressures,  the  effect  of  thread 
stretching  dominates,  and  smaller  jm' s  may  be  expected. 

The  limiting  value  of  unity  is  obtained  for  the  re¬ 

sults .  This  suggests  that  the  material  is  deforming  as  a  pin- 
jointed  trellis:  At  higher  pressures,  thread  stretching  and 
friction  between  warp  and  fill  threads  cause  to  decrease 
somewhat. 

A  comparison  of  the  series  2  results  with  those  of  series 
in  which  slightly  thicker  latex  coating  was  used,  suggests  little 
change  in  warp  and  fill  extensional  stiffnesses.  Although  the 
cruciform  tests  gave  no  reliable  quantitative  results,  they  indi¬ 
cated  a  strong  influence  of  latex  thickness  on  A  similar 

influence  on  and  G^t  may  be  anticipated. 

The  effect  of  the  number  of  coatings  on  the  vibration  fre¬ 
quencies  of  rectangular  cantilever  models  was  investigated  in 
Ref.  12.  Adding  extra  coats  of  latex  decreases  the  model  fre¬ 
quencies,  but  also  adds  some  additional  stiffness,  as  evidenced 
by  the  fact  that  the  changes  in  frequencies  were  not  inversely 
proportional  to  the  square  root  of  the  total  masses. 

In  the  absence  of  internal  losses,  one  would  expect  the 
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relation 
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to  be  satisfied,  if  the  material  is  assumed  to  behave  linearly. 
Table  2.14  presents  the  ratio  ///^  The  closer  this 

ratio  is  to  one,  the  better  the  experimental  results  satisfy  the 
above  relation.  Except  for  the  isolated  starred  point  in  the 
table  which  is  questionable,  this  ratio  stays  within  12  per  cent 
of  1.  This  is  well  within  the  expected  experimental  accuracy 
considering  the  fact  that  the  static  ^-measurements  are  not  be¬ 
lieved  to  be  very  accurate. 
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Figure  2.1  Sea-pits  of  Inflatable  Models. 
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Figure  2.2  Schematic  of  HodelPressure  Regulating 
System  end  Strain  Gages. 
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Figure  2.*.  Excitation  System  Used  for  Vibration  Testa* 
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Mode  2 


Mode  3 
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Figure  2.8  Typical  Free  Vibration  Node  Llnea  for  Model  10. 


Figure  2.9  Free  Vibration  Node  Lines  for  Model  18  • 
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Figure  2.12  Schematic  of  Kirror  Test  Setup. 
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Figure  2.13  Typical  Photograrah  for  Mirror  Mode  Shape  Data 
For  Data  Reduction,  a  15  x  24  in.  Print  is 
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Figure  2.14  Experimental  Mirror  Mode  Shape  for  First  Vibration 
Mode  of  Model  8  at  pj  »  2  pslg. 
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Lgure  2.13  Comparison  of  Experimental  Nod®  L*"** ^d(w! 

lines  Calculated  from  Measured  Influence  Ccei 
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t  WsrP  orientation  waa  0°  for  top  aurface  and  50°  for  bottoa  surface 


TDR  64-147 


TABLE  2.2  COMPARISON  OF  THEORETICAL  AND  EXPERIMENTAL 
DEFLECTIONS  OF  MODEL  1  UNDER  A  UNIFORM, 
STATIC  LOO.  MODEL  PRESSURE  IS  2  PSIG 

!  Station* 


Total  Load'  _  _ 'i  1 _ 2 _ 3 


2.33  lbs 

Theory  i 

Experiment 1 

0.164  in. 
0.22 

0.290  in. 
0.60 

WuEMm 

4.66 

Theory  ! 

Experiment  i 

0.328 

0.56 

0.580 

‘1.32 

0.484 

0.78 

6.99 

Theory  j 

Experiment  | 

0.492 

0.94 

0.870 

2.14 

0.726 

1.34 

#See  Fig.  2.3  for  location  of  stations 
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TABLE  2.4  EFFECT  OF  TEMPERATURE  ON  THE  FREE  VIBRATION 
FREQUENCIES  OF  MODEL  6  AT  fa  •  2  PSIG 


Mode  1  ] 

Mode 

2 

Mode  3 

■ 

Mode  5 

Temperature 

tx 

mm 

mm 

T(in.°F) 

(in.  cps) 

mm 

f2 

■9 

f3 

mm 

71 

70.5 

32.3 

70.5 

56.2 

70.5 

MMJ 

89 

87 

31.7 

86 

55.2 

85 

115 

mSBm 

109 

30.7 

112 

52.6 

111 

134 

15.4 

129 

30.8 

134 

51.8 

130 

74.0* 

^Questionable  value 

TABLE  2.5  COMPARISON  OF  EXPERIMENTAL  FREQUENCIES  AND 
ANALYTICAL  FREQUENCIES  BASED  ON  PURE-SHEAR- 
TYPE  DEFORMATIONS,  FOR  THREE  MODELS  AT  ft  ■ 
2  PSIG 


Experiment  | 

Simplified  Analysis  | 

Model  1 

8 

■BE 

Model  1 

8 

14.3 

16.7 

26.4 

18.2 

20.8 

19.3 

EH 

30.7 

33.7 

52.5 

36.6 

41.7 

38.8 

47.9 

52.6 

59.2 

50.5 

57.5 

53.4 

TABLE  2.6  EFFECT  OF  PAD  SIZE  ON  INFLUENCE 
COEFFICIENT  MEASUREMENTS 


c, . 

Lead  Pad  Diameter 

ij 

1  inch 

2  inch 

4  inch 

C54 

0.059 

w 

0.055 

C55 

0.228 

B 

0.212 

C56 

0,173 

0.174 

0.170 
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Station 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 


Mode  Amplitude 
(Normalized  with  Respect  to 


Station  10) 


Mode  1 

Mode  2 

0.00777 

-0.71470 

0.0999 

0.06226 

-1.20041 

0.0248 

0.14429 

-2.94835 

0.1623 

0.11510 

-0.86345 

-0.3390 

0.28120 

-2.64691 

-0.3287 

0.47852 

-2.73049 

0.1257 

0.15755 

0.11175 

-1.2230 

0.48997 

0.59971 

-1.4191 

0.74594 

1.61580 

-0.1107 

1.00000 

1.00000 

1.0000 

See  Fig.  2.17  for  station  locations 


TABLE  2.9 


COMPARISON  OF  EXPERIMENTAL  FREQUENCIES  AND 
FREQUENCIES  CALCULATED  FROM  MEASURED  IN¬ 
FLUENCE  COEFFICIENTS,  FOR  MODEL  5  AT  A 
PRESSURE  OF  2  PSIC. 


Frequency  (cps) 

Model  1 

Mode  2 

l  ode  3 

Experiment 

16.7 

32.3 

54.0 

Calculated 

15.9 

32.5 

46.4 

Experiment 

Calculated 
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TABLE  2.10  EXTENSIONS  STIFFNESSES  f  S.t  ,  AND 

ALONG  DIRECTIONS  OF  WAR?,  FILL  AND  45  DEGREES 
TO  WARP,  RESPECTIVELY 


Test 

Series 


Questionable  point 


TABLE  2.11  INFLUENCE  OF  TEMPERATCRE  AND  HUMIDITY  ON  THE 
EXTENSIONS  E"  IFFNESS  ,  AT  CYLINDER 

PRESSURE  OF  0.688  PSIC^ 


Humidity 


381 

V 


\wm 

Remarks 

133 

124 

Surface  watted 

118 

Vetted  and  heated  for  a  short 
period 

109 

Vetted  and  heated  for  a  longer 
period 

75 
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TABLE  2.13  POISSON'S  RATIOS  /4,  .  fj.  #  U  m 

>. f*'  i  •  r *s 


/W 

Test 

Series 

Method  of 
Test 

4  "2 

(psig) 

4 

6 

2 

Static 

0.343 

0.289 

0.279 

3 

Static 

0.545 

0.488 

0.439 

2 

4 

6 

2 

Static 

0.578 

0.412 

18SI 

3 

Static 

0.802 

0.641 

/ 

2 

4 

6 

2 

Static 

1.005 

0.947 

0.681 

3 

Static 

0.942 

0.815 

0.740 

TABLE  2.14  THE  RATIO  /V  ^ FROM  THE 
EXPERIMENTAL  MEASUREMENTS 


Test 

Series 

••  Method  of 

Test 

pd  -  2  psi 

4 

6 

2 

Static 

1.06 

1.12 

EH 

2 

Static  for  fU 
Dynamic  for  £ 

1.33* 

1.07 

■ 

3 

Static  for  ft* 
Dynamic  for  £ 

1.12 

1.07 

D 

* Questionable  point 
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SECTION  3 
WIND-TUNNEL  TESTS 


A.  WIND  TUNNELS 

The  supersonic  wind-tunnel  tests  were  conducted  in  Tunnel  A, 
Gas  Dynamics  Facility,  Arnold  Engineering  Development  Center 
(AEDC),  Tullahoma,  Tennessee.  This  is  a  40-  by  40-inch  continuous, 
closed-circuit,  variable-density  wind  tunnel  with  a  Mach  number 
range  of  1.5  to  6.0  at  maximum  stagnation  pressures  from  29  to 
200  psia,  respectively,  and  stagnation  temperatures  up  to  300  °F 
(M  -  6). 

The  transonic  wind-tunnel  tests  were  run  in  the  8-foot 
transonic  wind  tunnel  at  Cornell  Aeronautical  Laboratory,  Inc. 
(CAL),  Buffalo,  New  York.  This  is  a  pressurized  wind  tunnel 
capable  of  operating  from  1/6  to  2-1/2  atmospheres  and  over  a 
Mach  number  range  from  zero  to  approximately  1.4. 

A  single  incompressible  flutter  point  was  obtained  in  the 
7-1/2-ft  x  10-ft  Wright  Brothers  Wind  Tunnel  at  MIT.  This  is  a  • 
low-speed,  variable-density  wind  tunnel  with  a  maxlana  speed  of 
about  140  mph. 

B.  TESTING  PROCEDURES 

1.  Supersonic  Tests 

For  the  supersonic  tests,  the  models  were  mounted  on  a 
30-inch  circular  splitter  plate  which  was  designed  and  built  as 
part  of  this  program  (Fig.  3.1). 

It  was  originally  anticipated  that  some  difficulty  would 
be  experienced  from  the  dynamic  loads  imposed  on  the  models  by 
the  starting  and  stopping  shocks;  however,  the  models  were  flex¬ 
ible  enough  so  that  they  were  not  destroyed  by  these  starting 
loads  and  special  equipment  to  protect  the  models  was  not  re- 
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quired.  Some  difficulty  was  experienced  with  the  breaking  of 
strain  gages  during  wind-tunnel  starts,  but  thi3  problem  was 
alleviated  by  making  the  wind-tunnel  starts  at  the  minimum  stag¬ 
nation  pressures  and  by  increasing  the  number  of  strain  gages 
mounted  on  the  model  from  2  to  4,  so  as  to  increase  the  proba¬ 
bility  of  having  at  least  one  usable  gage  circuit  after  the  tun¬ 
nel  start. 

.  The  normal  test  procedure  was  to  establish  the  desired 

st'sady-state-test  conditions  in  the  tunnel  with  the  model  in¬ 
flated  to  a  known  pressure  differential'  above  tunnel  static  pres¬ 
sure.  The  model  was  then  excited  by  the  plucking  mechanism 
(Fig.  3.1);  if  flutter  occurred,  the  stagnation  pressure  was 
lowered  until  flutter  stopped;  if  the  model  did  not  flutter,  the 
stagnation  pressure  was  increased  to  obtain  a  new  steady-state- 
flow  condition  and  the  model  was  again  excited  with  the  plucking 
mechanism.  This  procedure  was  repeated  until  cither  flutter  oc- 
.  curred,  or  the  maximum  tunnel  stagnation  pressure  was  achieved. 

In  some  tests,  flutter  was  initiated,  or  checked,  by 
varying  the  model  pressure  while  maintaining  steady-state-tunsel 
conditions. 

A  record  of  the  outputs  of  the  strain  gages  and  the  pres¬ 
sure  differential  between  the  model  pressure  and  tunnel  static 
pressure  was  obtained  on  a  direct  writing  oscillograph.  During 
the  first  teats,  some  flutter  frequencies  were  not  obtained  be¬ 
cause  the  oscillograph  was  not  running  at  a  high  enough  speed  at 
the  time  flutter  occurred.  Therefore,  for  the  later  tests,  this 
data  was  also  continuously  recorded  on  magnetic  tape  to  provide 
a  aora  reliable  frequency  reading  capability. 

The  model  pressure,  the  tunnel  static  and  stagnation  pres¬ 
sures,  the  pressure  differential  between  the  model  pressure  and 
tunnel  static  pressure,  the  stagnation  temperature,  and  the  Mach 
number  were  all  recorded  and  used  in  conjunction  with  the  VKF 
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ERA  1102  computer  to  tabulate  tunnel  parameters  and  the  outputs 
of  the  pressure  transducers.  The  continuous  oscillograph  record 
was  tied  in  with  the  tabulated  data  (which  was  taken,  tabulated, 
and  printed  out  every  three  seconds)  through  the  use  of  a  manual 
"flutter-on-off"  switch  which  (1)  controlled  a  two-position  trace 
on  the  oscillograph  record,  and  (2)  caused  a  "0"  (flutter-off), 
or  a  "1"  (flutter-on)  to  be  printed  in  one  column  of  the  tabulated 
data.  Another  two-position  trace  on  the  oscillograph  recorded 
whenever  data  was  being  taken.  In  this  manner  it  was  possible  to 
determine  the  onset  or  termination  of  flutter  from  the  output  of 
the  stra in-gage* traces  on  the  oscillograph  records,  and  then  de¬ 
termine  the  corresponding  tunnel  and  model  conditions  from  the 
tabulated  data. 

Visual  records  of  the  flutter  motion,  which  provided  in¬ 
formation  on  flutter-mode  shapes,  were  obtained  using  high-speed 
motion  picture  cameras. 


2.  Transonic  Tests 

For  the  transonic  tests,  the  models  were  attached  to  an 
aluminum  base  plate  which  was  installed  on  the  reflection  plane 
cart  of  the  transonic  wind  tunnel  (Fig.  3.1).  A  rotating  mecha¬ 
nism  was  used  to  excite  the  model  in  the  transonic  tests. 

The  testing  procedures  were  essentially  the  same  as  those 
used  for  the  supersonic  tests. 


3.  Low-Speed  Tests 

For  the  low-speed  test  at  HIT  the  model  was  mounted  on  a 
plywood  board  which  was  placed  vertically  in  the  center  of  the 
tunnel  and  spanned  the  entire  height  of  the  tunnel.  For  these 
tests,  the  tunnel  density  was  atmospheric,  and  the  flutter  point 
was  determined  by  varying  tunnel  speed  while  maintaining  constant 
model  pressure.  Essentially  the  same  data  was  recorded  as  In  the 
supersonic  and  transonic  tests. 
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ERA  1102  computer  to  tabulate  tunnel  parameters  and  the  outputs 
of  the  pressure  transducers.  The  continuous  oscillograph  record 
was  tied  in  with  the  tabulated  data  (which  was  taken,  tabulated, 
and  printed  out  every  three  seconds)  through  the  use  of  a  manual 
"f lutter-on-off"  switch  which  (1)  controlled  a  two-position  trace 
on  the  oscillograph  record,  and  (2)  caused  a  "0"  (flutter-off), 
or  a  "1"  (flutter-on)  to  be  printed  in  one  column  of  the  tabulated 
data.  Another  two-position  trace  on  the  oscillograph  recorded 
whenever  data  was  being  taken.  In  this  manner  it  was  possible  to 
determine  the  onset  or  termination  of  flutter  from  the  output  of 
the  strain-gage ‘traces  on  the  oscillograph  records,  and  then  de¬ 
termine  the  corresponding  tunnel  and  model  conditions  from  the 
tabulated  data. 

Visual  records  of  the  flutter  motion,  which  provided  in¬ 
formation  on  flutter-mode  shapes,  were  obtained  using  high-speed 
motion  picture  cameras. 

2 .  Transonic  Tests 

For  the  transonic  tests,  the  models  were  attached  to  an 
aluminum  base  plate  which  was  installed  on  the  reflection  plane 
cart  of  the  transonic  wind  tunnel  (Fig.  3.1).  A  rotating  mecha¬ 
nism  was  used  to  excite  the  model  in  the  transonic  tests. 

The  testing  procedures  were  essentially  the  same  as  those 
used  for  the  supersonic  tests. 

3.  Low-Speed  Tests 

For  the  low-speed  test  at  MIT  the  model  was  mounted  on  a 
plywood  board  which  was  placed  vertically  in  the  center  of  the 
tunnel  and  spanned  the  entire  height  of  the  tunnel.  For  these 
tests,  the  tunnel  density  was  atmospheric,  and  the  flutter  point 
was  determined  by  varying  tunnel  speed  while  maintaining  constant 
model  pressure.  Essentially  the  same  data  was  recorded  as  in  the 
supersonic  and  transonic  tests. 

81 

CONFIDENTIAL 

l«l  ‘TiiWUwXl  f  Xkbp) 


FDL  TDR  64-147 


CONFIDENTIAL 


C.  FLUTTER  RESULTS  AND  DISCUSSION 

1.  General 

The  plucking  mechanisms  used  Co  excite  the  flutter  models 
in  the  transonic  and  supersonic  tests  (Fig.  3.1)  were  a  very  im- 
portant  part  of  the  test  procedure.  If  the  model  was  not  ex¬ 
cited,  it  was  possible  to  penetrate  deeply  into  the  unstable 
region  without  encountering  flutter. 

For  the  supersonic  tests,  once  flutter  was  initiated  by 
exciting  the  model  with  the  plucking  mechanism,  the  flutter 
motion  would  continue  for  many  cycles  (at  limited  amplitude)as 
the  stagnation  pressure  was  decreased,  even  in  regions  where  the 
model  had  previously  been  stable.  Figure  3.2  presents  photo¬ 
graphs  of  typical  strain-gage  responses  showing  (1)  response  to 
plucking  excitation  in  a  stable  region,  (2)  response  to  plucking 
excitation  in  an  unstable  region,  and  (3)  limited  amplitude 
flutter  continuing  while  the  stagnation  pressure  decreases,  until 
finally  the  flutter  motion  suddenly  subsides.  It  is  this  latter 
point,  where  the  motion  suddenly  ceases,  which  was  taken  as  the 
flutter  boundary.  For  the  supersonic  tests  the  model  was  ex¬ 
cited  by  an  up  and  down  movement  of  the  plucking  mechanism.  Be¬ 
cause  of  the  placement  of  the  model  with  respect  to  the  plucking 
mechanism,  the  downward,  or  second  excitation  gave  a  larger 
forcing  amplitude,  as  can  be  seen  from  the  strain-gage  responses. 
In  general,  flutter  was  Initiated  by  this  second  excitation,  with 
a  rapid  buildup  to  a  steady-state  amplitude  after  a  few  initial 
cycles  of  erratic  motion,  it  should  be  noted  that  the  shape  of 
the  decaying  portions  of  the  strain-gage  responses  are  not  ex¬ 
ponential  as  is  usually  the  case  for  linear  systems. 

The  inflatable  models  used  in  this  program  were  capable 
of  withstanding  many  cycles  of  limited  amplitude  oscillation 
without  damage.  For  the  supersonic  tests,  it  was  almost  always 
possible  to  enter  into  the  unstable  region,  initiate  flutter  by 
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exciting  the  model,  and  return  to  the  stable  region  without  de¬ 
stroying  the  model,  so  that  more  than  one  flutter  point  could  be 
obtained  from  one  model.  In  fact.  Model  7  alone  accounted  for 
25  flutter  points  before  it  was  finally  destroyed.  The  model 
destruction  usually  occurred  after  some  local  rupture  in  the 
model  surface  resulted  in  a  complete  loss  of  internal  air  pres¬ 
sure,  and  therefore  model  stiffness.  The  model  would  then  buf¬ 
fet  against  the  splitter  plate  and  be  destroyed.  Figure  3.3, 
taken  from  the  high-speed  motion  picture  film,  shows  a  model 
during  flutter  and  after  losing  pressure. 

Unlike  the  supersonic  tests  at  AEDC,  flutter  in  the  CAL 
transonic  tunnel  was  most  often  catastrophic.  The  models  usually 
failed  after  a  few  cycles  and  it  was  only  possible  to  get  eight 
flutter  points  from  six  models.  This  may  be  attributed  to  the 
fact  that  the  transonic  tests  were  generally  made  ar  higher 
values  of  model  pressure  differential  and  dynamic  pressures. 

Figure  3.4  shows  some  typical  models  after  failure  oc¬ 
curred  and  the  tunnel  had  been  shut  down.  Incidentally,  it  may 
be  noted  that  all  model  failures  were  not  necessarily  cata¬ 
strophic,  as  evidenced  by  Model  1  which  developed  a  large  bubble 
at  the  wing  tip  caused  by  drop-chord  failure  during  flutter. 

2.  Results  and  Discussion 

A  complete  summary  of  all  the  experimental  flutter  data 
is  given  in  Table  3.1,  which  also  includes  calculated  flutter 
parameters  based  on  model  properties  and  tunnel  conditions  at 
flutter. 

As  previously  mentioned,  the  vibration  characteristics 
of  the  inflatable  models  were  somewhat  variable;  vibration  tests 
made  on  the  same  model  at  different  times  gave  frequencies  which 
differed  by  wore  than  the  experimental  error.  The  flutter  pa¬ 
rameters  given  in  Table  3.1  arc  based  on  the  last  frequencies 
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measured  before  Che  actual  flutter  test.  In  some  Instances  the 
models  fluttered  at  values  of  model  pressure  differential  for 
which  no  vibration  data  exists.  For  these  cases,  the  vibration 
frequencies  were  estimated  by  interpolation  and  extrapolation  of 
the  available  vibration  data.  The  vibration  results  used  to 
reduce  the  flutter  data  for  particular  flutter  points  is  indi¬ 
cated  in  Table  2.3. 

Mo  corrections  have  been  made  in  the  tabulated  data  to 
account  for  temperature  or  apparent  mass  effects  through  their 
influence  on  the  measured  vibration  frequencies  (see  Section  2, 
Subsection  C) .  The  apparent  mass  effects  can  be  accounted  for 
by  appropriate  modification  of  the  still-air  vibration  tests, 
based  on  the  results  of  the  vibration  tests  carried  out  on 
Model  2  under  reduced  ambient  pressures.  This  effect  is  small 
and  is  considered  to  be  negligible  for  the  present  program.  A 
theoretical  method  of  correcting  for  temperature  effects,*  ac¬ 
counting  for  the  temperature  variations  occurring  over  the  model 
surfaces  during  the  wind-tunnel  tests,  is  beyond  the  state  of 
the  art.  This  approach  is  further  complicated  by  the  fact  that 
it  would  require  accurate  knowledge  of  material  structural  and 
thermal  properties  as  a  function  of  temperature .  A  more  fea¬ 
sible  approach  would  be  to  determine  the  actual  temperature  dis¬ 
tributions  by  appropriate  wind-tunnel  tests,  and  then  measure 
the  model  frequencies  through  vibration  tests  performed  under 
tl*  real  conditions  sinilating  those  measured  in  the  wind  tunnel. 
The  main  problem  encountered  here  would  be  that  of  simulating 
the  required  temperature  gradients  with  the  relatively  small- 
models  iiged  in  the  present  program.  While  the  uniform  tempera- 

Por  models  where  only  shear-type  deformations  are  of  importance, 
the  temperature  effects  will  be  nonexistent,  since  all  the  stiff¬ 
ness  is  derived  from  the  model  pressure  differential.  The  fact 
that  the  models  exhibit  s  marked  dependence  on  temperature  is 
further  indication  that  bending -type  deformations  must  be  con¬ 
sidered  in  any  analysis  carried  cut  on  the  present  type  models. 
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Cure  tests  carried  out  under  this  program  are  not  useful  for  any 
quantitative  evaluation  of  the  effects  of  temperature  on  the 
flutter  results,  they  do  serve  to  indicate  that  these  effects 
are  important.  The  effect  of  increased  temperature  Is  to  lower 
the  natural  vibration  frequencies.  Thus,  the  experimental 
flutter  parameters  Vf  /_J>  and  ,  which  are  based  on 

's  measured  at  room  temperature,  would  be  lower  than  the 
actual  values  for  stagnation  temperatures  considerably  above 
room  temperature.  This  particular  point  will  be  brought  up 
again  during  subsequent  discussions  of  the  flutter  data. 

Eight  of  the  flutter  points  given  in  Table  3.1  correspond 
to  repeated  tests  for  identical  model  configurations.  Except 
for  flutter  points  17,  22,  and  24,  each  pair  of  duplicate  tests 
was  made  during  the  same  tunnel  run.  The  results  are  repeatable 
to  within  approximately  1  to  1-1/2  per  cent  in  the  flutter  index, 
/b'Uii'TjZ  *  *nd  2  to  2-1/2  per  cent  in  the  flutter  frequency 
parameter,  o^/^,  which  is  within  the  accuracy  of  the  measure¬ 
ments.  For  flutter  points  17,  22,  and  24  which  all  correspond 
to  the  same  test  configuration,  the  maximum  variations  between 
the  flutter  parameters  are  somewhat  larger  (3.7  per  cent  for 

and  6.2  per  cent  for  Flutter  points  17  and 

22  were  obtained  during  the  same  tunnel  run  (at  the  beginning 
and  end  of  the  series  of  tests  made  to  determine  the  effects  of 
model  pressure  differential).  Flutter  point  24  was  obtained  at 
a  ;~h  later  time.  Model  7  having  been  removed  from  the  tunnel 
while  other  models  were  tested.  The  variation  in  the  flutter 
parameters  can  be  attributed  mainly  to  the  fact  that  the  stag¬ 
nation  temperature  for  flutter  point  24  is  40  degrees  higher 
than  that  for  points  17  and  22. 

In  order  to  determine  the  variation  of  the  flutter  index, 

JjZ *  and  the  flutter  frequency  parameter,  ,  with 

Mach  number,  a  large  number  of  tests  were  made  on  a  basic  model 
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configuration,  consisting  of  a  65-degree-sweep  delta,  canti¬ 
levered  model  with  a  model  pressure  differential  of  2  psig  and 
an  angle  of  attack  of  zero  degree.  The  value  of  2  psig  for 
model  pressure  differential  was  selected  cm  the  basis  of  pre¬ 
liminary  calculations  which  indicated  that  this  value  would  pro¬ 
vide  the  largest  Mach  number  range  over  which  model  instabili¬ 
ties  would  be  encountered,  within  the  operating  range  of  the 
wind  tunnel  used  in  the  supersonic  testa.  The  results  for  this 
basic  configuration,  in  terms  of  the  parameters  and 

vs  M,  are-plotted  in  Figs.  3.5  and  3.6.*  With  a  few  ex¬ 
ceptions  which  will  be  covered  later,  the  flutter  mode  for  the 
basic  configuration  was  a  combination  of  the  first  and  second 
vibration  modes,  and  only  this  type  of  flutter  mode  will  be 
considered  here.  Henceforth  this  will  be  referred  to  as  low-  • 
frequency  flutter. 

An  alternate  choice  of  dimensionless  parameters,  which 
arise  in  theoretical  calculations  based  on  piston  theory  or 
s lender-body  theory,  are  the  flutter  parameter  /boat  ,  the 
flutter  frequency  parameter,  u*  /u3l  ,  and  the  mass -ratio -Mach- 
number  product,  JL M.  '  Figures  3.7  and  3.8  present  the  results 
of  the  experiments  on  the  basic  model  configuration  in  terms  of 
these  parameters.  Again,  only  the  low  frequency  flutter  points  , 
ere  considered. 

Theoretical  curves,  based  on  the  first  two  vibration 
nodes  of  a  typical  model,  are  also  Included  in  Figs.  3.5  through 
3.8.  The  calculations  are  for  sero  structural  damping  in  all 
nodes  and  for  frequency  ratios  0.470  (for  ends  shapes 

A  and  B)  and  (  cJ,  /u^)  m  0.478  (for  mode  shapes  C-l  and  C-2). 
They- are  taken  fron  the  lower  boundaries  in  Figs.  1.7  through 
1.12  and  the  boundary  in  Fig.  1.16  of  Appendix  I. 

*  Refer  to  Table  2.1  for  identification  of  synbols  assigned  to 
nodels  for  use  in  plotting  flutter  data. 
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Theoretically,  the  variation  of  the  flutter  index, 

,  with  Mach  number  is  not  independent  of  jx.  ,  so 
that  one  cannot  strictly  draw  a  single  theoretical  curve  of 
vs  ***  as  *s  ^one  *n  Figs .  3.5  and  3.6.  Owing  to 
the  fact  that  the  anb'ent  density  becoaes  progressively  lower  as 
M  is  raised,  the  mass  ratio,  y.  ,  at  flutter  increased  by  a  factor 
exceeding  8  at  the  supersonic  and  transonic  ends  of  the  speed 
range.  Since  this  behavior  seems  unavoidable  in  such  experi¬ 
ments,  it  is  fortunate  that  a  large  amount  of  theoretical  and 
experimental  evidence  exists  that  shows  that  the  dependency  of 
Vp/btJi'fF’  uPon  /-*•  very  weak,  particularly  at  supersonic 
speeds  and  for  the  range  of  y- ' s  considered  here.  This  point 
is  further  illustrated  by  the  closeness  of  the  calculations 
shown  in  Figs.  3.5  and  3.6  for  JjLm  71.7  and  y-  •  317. 

The  analysis  using  mode  shape  sets  A  and  B  were  conducted 
early  in  the  program  and  preceded  the  initial  flutter  tests  on 
Models  3,  4,  7,  and  8.  As  pointed  out  in  Appendix  I,  the  first 
mode  is  the  same  for  mode  shape  sets  A  and  B.  The  second  modes 
are  different,  but  they  both  satisfy  the  node  line  locations 
obtained  from  Model  1  vibration  data,  and  variations  between 
these  two  second  modes  are  seemingly  Mall  (see  Figs.  1.4  and 
1.5).  However,  the  theoretical  flutter  boundaries  are  quite 
separated,  that  of  set  A  differing  from  that  of  set  B  by  as 
uich  as.  15  per.  cent. 

The  experimental  flutter  parameters  /b, ,/T  and 
for  Models  3,  4,  and  8  are  bracketed  by  the  two  theo¬ 
retical  boundaries,  obtained  using  node  shape  sets  A  and  B,  each 
boundary  giving  better  correlation  at  a  different  range  on  the  M 
or  yM  scale;  set  A  is  better  at  the  higher  end  (H  ^  5,  1200) 

and  set  B  is  better  for  the  middle  range  (M  *  3,  y M  >  800) .  The 
experimental  frequency  ratios,  C±/toL  •  compare  more  favorably 
with  calculations  based  on  mode  shape  set  B. 


CONFIDENTIAL 


FDL  TDR  64-147 


CONFIDENTIAL 


As  Che  experimental  vibration  code  shape  obtained  froa 
the  mirror  tests  on  Model  8  became  available,  two  additional 
sets  of  calculations  were  performed  (sets  C-l  and  C-2).  In  these 
calculations  the  nonuniformity  of  the  mass  distribution  was  ac¬ 
counted  for  approximately  (see  Appendix  X,  Section  B).  Addi¬ 
tionally,  in  the  analysis  with  set  C-2,.  an  attempt  was  made  to 
estimate  the  influence  of  the  blunt  leading  and  trailing  edges 
on  flutter,  by  simply  omitting  the  aerodynamic  contributions 
from  these  regions.  Although  such  an  empirical  correction  is  a 
rough  approximation,  particularly  when  the  leading  and  trailing 
edges  are  not  circular  in  planes  parallel  to  the  free  stream, 
one  would  have  expected  that  the  experimental  results  for  Model  8 
would  lie  somewhere  between  the  boundaries  of  C-l  and  C-2.  How¬ 
ever,  this  is  not  the  case. 

The  results  for  Model  7  are  not  consistent  with  those  for 
the  other  models.  This  lack  of  consistency,  especially  between 
Models  7  and  8,  is  somewhat  surprising.  Models  7  and  8  were 
built  to  be  identical;  they  differed  slightly  in  mass  and  in 
frequencies,  but  displayed  identical  node-line  patterns.  The 
fact  that  the  Model  7  tests  were  run  at  slightly  lower  stagna¬ 
tion  temperatures  cannot  account  for  the  difference  in  the  level 
of  the  boundaries  for  Models  7  and  8.  Because  of  the  strong  de¬ 
pendence  of  the  flutter  results  on  mode  shapes  as  indicated  by 
the  theoretical  calculations,  the  only  plausible  explanation 
that  can  be  offered  here  is  that  Models  7  and  8  did  have  slightly 
different  mode  shapes  and  that  this  critically  affected  the  flutter 
results. 

A  comparison  can  also  be  made  as  to  hw  veil  theory  predicts 
the  trends  of  the  various  flutter  parameters  with  respect  to  M  and 
/-CM.  Recalling  that  (1)  the  experimental  flutter  parameters 

•  •ad  *£/*>*.  »  Which  are  based  on 
measured  at  room  temperature,  would  be  lower  than  the  true  values 
for  stagnation  temperatures  considerably  above  room  temperature 
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and  would  be  higher  Chan  Che  true  values  for  stagnation  tesapera- 
Cures  ac  or  below  room  temperature,  (2)  the  high  Mach  number 
tests  were  conducted  at  stagnation  temperatures  well  above  rooa 
temperature,  and  (3)  the  lower  Mach  number  tests  were  conducted 
at  stagnation  temperatures  at  or  below  rooa  temperatures,  the 
experimental  results,  corrected  for  this  temperature  effect, 
would  be  slightly  lower  at  the  lower  Mach  numbers  and  slightly 
higher  at  the  higher  Mach  numbers.  This  would  tend  to  Increase 
the  slope  of  the  experimental  data  presented  in  Figs.  3.5  through 
3.8.  Therefore,  the  ability  of  the  theory  to  correctly  indicate 
the  trends  is  ^  Vf/b,uL  w^t*1  M  is  ex¬ 

pected  to  be  somewhat  better  than  that  indicated  by  Figs.  3.5  and 
3.7. 

As  shown  in  Figs.  3.6  and  3.8,  ^ /u)z ,  based  on  the  ex¬ 
perimental  results  uncorrected  for  temperature  effects,  is  es¬ 
sentially  independent  of  M  and  jj.  M,  as  predicted  by  theory. 
Correcting  for  temperature  effects  would  lessen  thi  agreement 
between  theory  and  experiment.  Since  it  is  a  natter  of  common 
experience  that  flutter  frequency  is  predicted  less  accurately 
than  the  flutter  speed,  the  correlation  with  regard  to  the 
trends  in  the  flutter  frequency  ratio  is  considered  adequate. 

The  oely  law  Mach  number  point  obtained  at  a  model  pres¬ 
sure  differential  of  2  psig  was  for  Model  4  at  a  Mach  number  of 
0.116.  No  calculations  were  carried  out  for  this  case. 

At  the  lower  Mach  numbers  and  in  the  transonic  tests  at 
CAL,  the  65-deg  models  could  not  be  flown  in  a  stable  region  at 
model  pressure  differentials  of  2  psig,  and  higher  pressure 
differentials  (up  to  10  psig)  were  required.  The  results  for 
eodel  pressure  differentials  equal  or  greatar  than  3  psig  are 
plotted  in  Figs.  3.9  through  3.12.  The  results  are  for  low 
frequency  flutter  only,  and  the  supersonic  values  are  pre¬ 
dominantly  from  Model  7.  The  high  pressure,  high  Mach  number 
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tesCs  were  made  to  determine  the  effects  of  model  pressure,  which 
will  be  discussed  later,  and  are  included  here  as  also  being  ap¬ 
propriate. 

In  general,  the  discussions  presented  in  the  preceding 
paragraphs  apply  to  the  results  of  Figs.  3.9  through  3.12,  the 
correlation  between  theory  and  experiment  being  about  the  same 
as  for  the  basic  model  configuration  with  model  pressure  differ¬ 
entials  of  2  psig.  It  should  be  noted  that  the  theoretical 
boundaries  used  in  Figs.  3.9  through  3.12  are  based  on  mode  shape 
data  obtained  at  a  pressure  differential  of  2  psig.  Again,  be¬ 
cause  of  the  marked  dependence  of  the  theoretical  results  on 
mode  shapes,  no  particular  significance  can  be  attached  to  the 
fact  that  the  experimental  results  agree  with  calculations  based 
on  mode  sets  S  and  C-l,  for  the  supersonic  region,  and  with 
mode  shape  set  A  for  the  transonic  region. 

A  sharp  rise  in  the  boundary  of  \ £  /bu^'/juu  vs  M,  at  Kach 
numbers  slightly  above  transonic,  is  evident  from  both  theory  and 
experiment.  From  the  experimental  data,  the  flutter  index, 
yf/b0o3L'fj£‘  i*  essentially  constant  for  subsonic  speeds,  with  no 
dip  in  the  transonic  region. 

The  effect  of  model  pressure  differential  on  flutter  of 
Model  7  is  shown  in  Fig.  3.13.  Increasing  model  pressure  dif¬ 
ferential  from  2  to  4  psig  causes  a  slight  decrease  in  the 
flutter  index  "Vf. / •  Further  increase  in  pressure  to 
6  psig  results  In  a  leveling  off  of  the  curve  for 

M  -  2.98,  and  a  sudden  Increase  in  the  curve  for  M  "  4.99.  .  The 
same  type  of  behavior  is  exhibited  by  the  frequency  ratio  CO^/ui%  . 

The  M  •  2.50  and  2.98  tests  were  all  made  at  constant 
stagnation  temperature  (85  to  87  °F),  so  that  correction  for 
temperature  effects,  if  any,  would  be  uniform  over  the  range  of 
pressures.  The  M  ”  4.99  tests  were  conducted  at  higher  stagna¬ 
tion  teoperatures  (121  to  175  °F),  with  the  higher  teaperatures 
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assoclaCed  with  the  higher  pressures.  In  this  case,  any  cor¬ 
rection  for  terperature  effects  would  raise  the  M  ■  4.99  re¬ 
sults  (relative  to  the  lower  Mach  number  results)  and  accentuate 
the  sudden  increase  for  the  pressure  differencials  of  5  and 
6  psig. 

As  seen  froo  Table  3.1,  the  vibration  measurements  made 
prior  to  the  M  -  2.98  tests  showed  that  the  frequencies  in¬ 
creased,  but  the  frequency  ratio  u/,Ju)^  decreased,  with  increas¬ 
ing  model  pressure  differential.  VibratiQn  measurements  made 
prior  to  the  M  »  4.99  tests  also  showed  increases  in  fre¬ 
quencies,  but  at  a  constant  frequency  ratio,  as  the  pressure 
differential  was  varied  frdm  2  to  6  psig.  Theoretically,  a 
slightly  lower  value  for  £>,/&>*  results  in  a  small  increase  in 
V / (aee  Fi8*  1.15).  Tor  the  present  case, 
these  effects  are  negligible. 

Finally,  mode-shape  effects  are  known  to  be  important, 
and  probably  account  for  the  major  portion  of  the  variation  of 
the  flutter  parameters  with  pressure  differential.  Since  it  is 
expected  that  mode-shape  effects  would  be  similar  for  super¬ 
sonic  Mach  numbers,  the  different  behavior  at  M  «  2.98  and  4.99 
cannot  be  explained  on  this  basis. 

The  variations  of  Jf?  and  with  the  root 

angle  of  attack  are  presented  in  Fig.  3.14.  Theoretical  calcu¬ 
lations  for  comparison  purposes  were  not  carried  out,  since 
many  complicating  factors  must  be  considered  in  the  analysis. 

The  most  important  are  (1)  changes  in  vibration  mode  shapes  and 
frequencies  which  are  certain  to  occur  when  the  model  deforms 
under  the  steady  aerodynamic  loads  at  angle  of  attack,  and  (2) 
the  variations  in  the  local  angles  of  attack  due  to  the  steady- 
state  deformations  and  their  effect  on  the  unsteady  aarodynaadc 
loads.  The  typically  large  deformations  experienced  by  the 
present  models  are  shown  in  Fig.  3.15  for  root  angles  of  attack 
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up  to  15  degrees,  at  H  ■  5.*  Although  models  were  flown  at 
steady-state  root  angles  of  15  degrees,  no  flutter  data  was 
sought  beyond  12  degrees,  because  of  the  excessive  deformations 
encountered  at  12  degrees  and  dynamic  pressures  required  for 
flutter. 


For  a  typical  section.  Ref.  13  demonstrates  that  the  ef¬ 
fect  of  angle  of  attack  on  flutter  is  destabilizing  according  to 
third  order  piston  theory.  For  a  flat  plate,  the  Influence  of 
angle  of  attack  may  be  estimated  by  interpreting  the  v* 

results  (for  oC  ■  0)  through  the  use  of  an  effective  jl  M 
according  to 


_  (j*  ^)iu£uuj,  <x„?o 
It  &L  MW 


where  cct  is  the  angle  of  attack  of  the  flat  plate  and  Its  assumed 
constant  over  the  plan form.  A  similar  correction  could  be  applied 
to  the  present  models,  if  the  model  possessed  a  constant  angle  of 
attack  in  its  steady  deformed  state.  This  is  not  the  case,  as  is 
evident  from  Fig.  3.15.  To  include  the  effects  of  angle  of  attack, 
when  the  local  angle  varies  over  the  planform,  one  must  start  from 
the  basic  pressure  equation  (Eq.(16),  Ref.  13).  This  approach 
would  require  a  priori  knowledge  of  the  deformed  state  which  would 
have  to  be  obtained  by  a  static  aeroelastic  analysis  or  appropriate 
experiments.  On  the  basis  of  the  preceding  discussion,  it  is  not 
evident  whether  the  effect  of  angle  of  attack  would  be  stabilising 
or  destabilizing  in  the  present  case. 

*The  dynamic  pressures,  under  which  the  photographs  shown  in 
Fig.  3.15  were  taken,  were  about  two  orders  of  magnitude  larger 
than  chose  which  would  be  experienced  during  an  actual  re-entry 
flight  (Fig.  1.3).  It  is  interesting  to  note  the  extremely 
large  deflections  which  these  models  are  capable  of  sustaining 
without  collapse  or  subsequent  change  in  structual  characteristics. 
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The  experimental  results  indicate  that  /bt re- 

nuins  constant  for  M  -  2.99  and  V,  «  4  psig,  and  that  it  in- 
<  *  * 

creases  substantially  with  angle  of  attack  for  M  «  4.99  and 

fa  ■  2  psig.  It  appears  that  the  effect  of  increasing  oc  is 
slightly  stabilizing. 

In  addition  to  the  65 -degree- sweep-delta  models,  a  number 
of  tests  were  conducted  on  five  75-degree  models.  Of  these. 

Models  15,  16,  and  17  were  observed  to  flutter  in  the  low-fre¬ 
quency  mode.  The  corresponding  experimental  results  are  presented 
in  Fig.  3.16.  The  tests  were  limited  to  Mach  number  3,  or  below. 

No  parallel  calculations  were  made  for  these  cases;  throughout 
the  low  supersonic  range  to  M  •  3,  the  leading  edges  ere  subsonic 
and  piston  theory  cannot  be  expected  to  yield  accurate  flutter 
.  predictions . 

When  compared  with  the  corresponding  curve  for  the  65-degree 
•*  models,  one  notes  that  the  rapid  rise  of  the  flutter  boundary  oc¬ 

curs  at  higher  values  of  M  for  the  75-degree  models.  Also,  for 
this  low  supersonic  range,  the  level  of  1*  slightly 

t  below  that  for  the  65-degree  models.  This  may  be  .  artially  due 

to  the  higher  frequency  ratios  obtained  with  the  lower-aspect- 
ratlo  models  (  47//*^  -  0.546  to  0.582);  the  effect  of  higher 

frequency  ratios  is  destabilizing  in  the  supersonic  range  (see 
Fig.  1.15).  In  the  transonic  region,  the  results  agree  well  with 
the  data  for  the  higher-aspect-ratio, models. 

A  second  type  of  flutter,  involving  modes  higher  than  the 
first  and  second,  which  will  be  referred  to  as  high-frequency 
flutter,  was  obtained  with  two  of  the  65-  and  75-degree  models; 
specifically.  Models  9,  10,  18,  and  20.  For  the  high-frequency 
flutter,  it  was  very  difficult  to  determine  the  flutter  boundary. 
Some  motion  persisted  throughout  all  of  the  tests,  as  could  be 
seen  from  the  response  of  the  strain  gages  mounted  on  the  models. 
Since  the  strain  gages  could  not  be  used  to  determine  the  boundary. 
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recourse  was  made  to  visually  observing  the  model  and  determining 
the  tunnel  conditions  for  which  some  motion  was  just  discernible. 
All  of  the  high-frequency  flutter  data  given  in  Table  3.1  and 
presented  in  Fig.  3.17  was  determined  in  this  way.  The  results 
for  Model  9  are  typical  of  the  behavior  of  the  models  which 
fluttered  in  the  high-frequency  mode,  and  the  following  discus¬ 
sion  will  be  mainly  restricted  to  this  particular  model. 

Referring  to  the  Model  9  flutter  data  given  in  Table  3.1 
for  ^  -  2  psig,  M  -  4.98,  and  HAfe"  1.11,  the  experimental 
value  of  the  flutter  parameter  is  7.21.  This  compares 

veil  with  the  theoretical  prediction  of  7.50,  based  on  the  re¬ 
sults  presented  in  Fig.  1.17  for  («j/“i.)*  •  1.27  (or  «1.13) 

and  at  the  experimental  value  for  jiH  of  1115.  This  value  of 
Vf  is  considerably  lower  than  that  for  low-frequency 
flutter.  It  is  a  well-recognized  fact,  that  the  closeness  in 
frequency  of  two  natural  modes  greatly  influences  the  location  of 
the  flutter  boundary,  and  the  most  critical  flutter  boundary  is 
most  often  the  one  based  on  those  modes  which  are  the  closest. 
Such  is  the  case  here  for  modes  2  and  3,  when  ^ *  2  psig.  As 
the  pressure  differential  is  increased  to  4  psig,  the  flutter 
frequency  is  very  near  that  for  the  fifth  mode,  indicating  strong 
coupling  of  the  fifth  node  with  either  a  close  fourth  or  a  close 
sixth  mode. 

Referring  to  Fig.  3.17,  it  is  seen  that  raising  the  pres¬ 
sure  from  2  to  4  psig  at  M  -  4.98  results  in  a  large  Increase  in 
the  critical  value  of  ,  mainly  as  a  consequence  of 

the  spreading  apart  of  the  frequencies  of  modes  2  and  3  (see 
Table  3.2,  Model  9,  Test  No.  3).  For  M  -  3  and  -fa  -  5  psig,  the 
value  of  the  critical  V 'fi'  somewhat  below  the  flutter 
results  based  on  lew- frequency  flutter.  A  rise  in  fa  to  8.10 
psig  is  accompanied  by  a  significant  rise  in  %/btUi2  JjZ  i  In 
fact,  the  value  of  is  in  excellent  agreement  with 
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the  low-frequency  flutter  results  for  the  other  65-degree  models. 
For  these  two  tests  (and  also  the  Model  9  test  at  M  -  2.48),  the 
recorder  monitoring  the  strain-gage  outputs  was  run  at  too  low  a 
paper  speed,  and  the  frequencies  could  not  be  counted  from  the 
record.  What  is  certain,  however,  is  that  the  flutter  frequencies 
are  above  200  cps,  and  that  the  flutter  frequency  for  fam  8.10 
psig  is  higher  than  that  for  fa  •  5.13  psig.  This  indicates  that 
modes  substantially  higher  than  the  fifth  mcd '  are  involved.  For 
the  M  •  2.48,  fa  »  8.42  psig  case,  the  str* in-gage  records  were 
not  at  all  clear.  Bursts  of  low-frequency  oscillations  are  evi¬ 
dent,  but  are  not  clean  enough  for  the  frequency  to  be  counted. 
Apparently,  at  this  pressure  the  model  may  be  close  to  flutter  in 
the  low-frequency  mode,  and  further  increases  in  ^  may  h«.ve  re¬ 
sulted  in  a  low-frequency-type  flutter. 

The  amplitude  of  the  high-frequency  flutter  was  too  low 
for  the  mode  shape  to  be  determined  from  the  high-speed  movies. 
From  visual  observations,  the  flutter  mode  had  two  node  points 
on  the  leading  edge.  The  first  was  approximately  midway  on  the 
leading  edge,  and  the  second  was  about  1  1/2  to  2  inches  up  from 
the  trailing  edge.  The  high-frequency  flutter  invariably  com¬ 
menced  with  a  pure  tip  motion,  the  node  line  extending  back  from 
the  second  node  point  on  the  leading  edge  at  en  angle  of  45  degrees 
with  respect  to  the  streamwise  direction.  It  was  the  motion  of 
these  few  square  inches  of  model  tip  which  was  observed  while 
visually  determining  the  onset  or  termination  of  flutter. 

The  results  for  Models  17  and  18  are  interesting  in  that 
the  models  were  built  to  be  identical,  yet  Model  17  fluttered  in 
the  lew  frequency  mode  while  Model  18  fluttered  in  the  high  fre¬ 
quency  mode.  Table  2.3  reveals  significant  differences  in  the 
vibration  results  for  the  higher  modes.  It  would  have  been  quite 
informative  if  Model  18  could  have  been  tested  at  the  high  pres¬ 
sures  that  Model  17  was  tested  at,  to  determine  if  the  flutter 
mode  would  indeed  have  changed  from  a  high-frequency  to  a  low- 
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frequency  one,  as  postulated  earlier.  Unfortunately,  Model  18 
was  destroyed  while  testing  at  the  lower  'a. 

Aa  a  final  presentation,  all  of  the  low-frequency  flutter 
data  have  been  combined  and  the  results  for  various  classifications 
are  presented,  in  terms  of  "^/b^z'/f^  vs  M,  in  Fig.  3.18.  In  lieu 
of  plotting  the  individual  points,  the  data  are  presented  as  bands 
encompassing  the  experimental  points  pertinent  to  one  of  the  par¬ 
ticular  classifications  shown  in  the  figure.  It  is  noted  here  once 
again  that  no  corrections  have  been  introduced  in  Fig.  3.18  to  ac¬ 
count  for  apparent -mass  and  thermal  effects. 

From  the  figure,  it  la  evident  that  the  experimental  re¬ 
sults  are  quite  consistent,  considering  the  variations  encountered 
between  some  of  the  models.  All  of  the  supersonic  results  (for 
M  ^  2.5)  lie  in  a  band  of  4  8  per  cent  about  a  median  experimental 
line.  The  region  from  M  ■  1.2  to  2.5  is  not  too  well  covered,  but 
the  trend  for  the  65-degree  models  is  reasonably  well  defined.  It 
is  also  clear  that  the  75-degree  boundary  lies  below  that  for  the 
65-degree  models  over  the  Mach  number  range  from,  say,  1.2  to  3.0. 
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Figure  3.2  Typical  Strain-Gage  Responses. 
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Figure  3.4  Photographs  of  Typical  Flutter  Models  after 
Destruction  IXiring  Flutter, 
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Figure  3.3  Experimental  Fitter  Index  '(fg/b^yS)  **  M.  and 
Comparison  with  Theory,  for  65-deg  Delta-Wing 
Models  at  a  Model  Pressure  Differential  of  pd  • 

2  pslg  and  o  •!  deft* 
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Figure  3.7  Experimental  Flutter  Parameter  (V^/l^o^)  vs  ***** 
and  Cooparison  with  Theory,  for  65-deg  Delta- 
Wing  Models  at  a  Modal  Pressure  Differential  of 
pd  -  2  psig  and  a  •  0  deg. 
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Figure  3.8  Experimental  Flutter  Frequency  Parameter  (o^/cdj) 
ve  uM,  and  Comparison  with  Theory,  for  65-deg 
Delta-Wing  Models  at  a  Model  Pressure  Differential 
of  pj  "  2  psig  and  o  »  0  deg. 
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Figure  3.9  Experimental  Flatter  Index  V1?*)  v*  **»  •“<* 

Cowparison  with  Theory,  for  65-deg  Delta-Wing 
Xodels  at  Model  Treasure  Differentials  of  ptf  2. 

3  paig  and  o  •  0  deg. 
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Figure  3.10  Experimental  Flutter  Frequency  Parameter  (m^/coj) 
vs  M,  and  Comparison  with  Theory,  for  65-deg 
'  Delta -Wing  Models  at  Model  Pressure  Differentials 
of  pd  3  p*ig  and  a  -  0  deg. 
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.  Figure  3.12  Experimental  Flutter  Frequency  Parameter  (u^/tuj) 
vi  |iH,  and  Comparison  with  Theory,  for  65-deg 
Delta-Wing  Models  at  Model  Pressure  Differentials 
of  Pj  1  3  P«U  and  o  •  0  deg. 
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Figure  3.15  Photographs  Showing  Steady-State  Deflections 
of  Model  8  at  Several  Angles  of  Attack  for 
M  «  5.  Portions  of  Root  and  Trailing  Edges 
are  indicated  by  Dotted  Lines. 
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Figure  3.18  liptrlatBtil  Boundaries  vs  M 

and  Coaparlson  with  Theory  for  all  Low* 
Frequency  Flutter.  All  Theory  Based  on 
63-d e*  Model  Data  at  pd  •  2  pslf . 
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SECTION  4 

CONCLUSIONS  AND  RECOMMENDATIONS 


The  main  conclusions  and  recommendations  that  may  be  drawn 
from  the  experimental  and  theoretical  results  of  this  study  are 
summarized  below. 

The  model  construction  technique  presented  here  is  simple 
and  yields  relatively  inexpensive  models.  It  is  particularly 
applicable  when  only  a  few  models  ere  desired  and/or  when  differ¬ 
ent  surface  materials  are  to  be  used  in  the  individual  models. 

In  general,  it  is  difficult  to  build  duplicate  models. 
Factors  which  should  lje  given  very  careful  attention  are:  (1)  the 
orientation  of  the  warp  and  fill  directions  of  the  surface  ma¬ 
terial,  for  both  single  and  multiple  layer  models,  (2)  the  type  of 
seams  used  along  the  leading  and  trailing  eiges,  and  (3)  the 
number  and  thickness  of  the  elastomer  coatings. 

The  latex-coated  nylon  and  dacron  cloth  used  as  the  cover 
material  have  been  found  to  be  sensitive  to  both  temperature  and 
humidity.  Their  properties  also  vary  slightly  with  age.  Since 
the  thermal  state  of  a  flutter  model  cannot  be  predicted  ac¬ 
curately  and  thus  compensated  for,  the  variations  of  model  struc¬ 
tural  behavior  with  temperature  are  undesirable  from  the  stand¬ 
point  of  flutter  testing  to  evaluate  analytical  techniques. 

Less  temperature-sensitive  cover  smterials  should  be  considered 
for  future  flutter  models  of  this  type. 

the  procedures  used  for  static  and  dynamic  testing  of  in¬ 
flatable  models  are  similar  to  those  used  in  programs  employing 
more  conventional-type  models.  With  the  exception  of  the  materials 
tests  and  the  mirror  mode  shape  tests,  no  new  or  special  techniques 
appear  to  be  required. 

In  vibrating  the  model  with  the  audio  system,  care  must  be 

119 


CONFIDENTIAL 


FDL  TDR  64-147 


CONFIDENTIAL 


exercised  to  not  overexcite  the  models  and  to  properly  position 
the  speaker  so  as  to  detect  all  the  weak  modes. 

Due  to  the  lightness  of  the  present  models,  the  apparent 
mass  effect  in  free-air  vibrations  is  noticeable,  amounting  to  a 
few  per  cent  for  the  first  mode  and  decreasing  in  importance  for 
the  higher  modes.  The  effect  may  be  easily  estimated  and  ac¬ 
counted  for  whenever  desirable. 

With  more  effort,  the  mirror  tests  can  be  made  to  yield 
mode  shapes  which  will  be  sufficiently  accurate  for  use  in  flutter 
calculations.  The  particular  advantage  of  mirror  mode  shapes  is 
that  the  streaovise  slopes,  which  enter  into  the  calculation  of 
the  generalized  aerodynamic  forces,  are  obtained  directly  from 
the  measurements,  while  the  deflections,  which  are  also  needed, 
are  obtained  by  integration. 

The  measurement  of  influence  coefficients  poses  no  special 
problem.  From  influence  coefficients,  one  may  calculate  the  mode 
shape  amplitudes  at  specified  points,  and  then  determine  the 
streamwlse  slopes  by  differentiation.  For  cases  where  the  theo¬ 
retical  result.",  appear  particularly  sensitive  to  the  choice  of 
assumed  modes,  a  more  laborious,  but  a  more  accurate  approach  is  ■ 
recommended.  This  approach  would  utilize  the  streamwlse  slopes 
from  the  mirror  mode  shape  tests  and  the  modal  amplitudes  from 
the  calculations  based  on  the  measured  influence  coefficients  to 
determine  the  generalized  mass  and  aerodynamic  forces. 

Any  theoretical  calculations  used  to  predict  the  static  or 
dynamic  behavior  of  Inflatable  delta-wing  models  of  the  type  pre¬ 
sented  here  should  allow  for  both  shear-  and  bending-type  de- 
.  formations.  The  effect  of  any  heavy  or  stiff  seams  along  the 
leading  and  trailing  edges  should  also  be  accounted  for. 

The  stiffness  properties  of  the  model  surface  material 
oust  be  accurately  determined  before  the  effects  of  bending-type 
deformations  can  be  included  in  vibration  or  other  aeroelastlc 
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calculations.  Dynamic  testing  methods  applied  to  cylindrical 
test  specimens  appear  to  offer  an  alternate  and  better  approach 
than  the  static  methods  now  being  used  to  determine  the  exten¬ 
sions  1  and  shear  moduli. 

The  predicted  flutter  speed  boundaries  are  in  fair  agree¬ 
ment  with  the  experimental  results.  Various  calculated  boundaries, 
based  on  different  sets  of  assumed  modes,  give  better  correlation 
over  different  Mach-number  ranges;  furthermore,  they  exhibit  suf¬ 
ficiently  different  levels,  so  as  to  suggest  significant  influences 
from  the  mode  shapes.  A  quantitative  statement  as  to  how  well  the 
analysis  predicts  flutter  of  inflatable  delta  wings  is  precluded 
by  this  dependence  of  the  theoretical  calculations  on  the  assumed 
mode  shapes.  The  main  problem  is  one  of  correctly  accounting  for 
the  camber  effects,  which  are  known  to  be  important,  and  which 
are  especially  dominant  in  the  case  of  low-aspect -ratio  wings 
(Ref.  10). 

Two  types  of  flutter  were  encountered;  the  first  involved 
the  two  lowest  vibration  modes  (low-frequency  flutter),  while  the 
second  involved  modes  higher  than  the  second  (high-frequency 
flutter) .  The  high-frequency  flutter  occurred  for  models  having 
two  very  close  higher  frequencies. 

The  influence  of  model  pressure  differential  appears  to  be 
very  slight  for  the  case  of  low-frequency  flutter.  In  contrast, 
the  effect  Is  much  more  significant  for  high-frequency  flutter. 

With  Increasing  pressure  differential,  higher  flutter  frequencies 
are  noted,  indicating  a  shift  to  higher  predominating  modes.  This 
shift  results  in  a  rise  of  the  flutter  index  ,  until 

such  time  as  the  level  reaches  that  of  the  low-frequency  flutter. 

From  the  experimental  tests  run  at  H  ■  2.99  and  4.99,  it 
appears  that  the  effect  of  increasing  the  angle  of  attack  Is 
slightly  stabilising  for  angles  of  attack  up  to  12  degrees.  It 
would  he  extremely  difficult  to  carry  out  a  theoretical  assess- 
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aent  of  thd  effect  of  angle  of  attack,  because  the  typlca/j^ 
large  deformations  experienced  by  the  aodela,  at  even  Low  *•<{!•• 
of  attack,  introduce  (1)  changes  in  vibration  mode  shape*  end 
frequencies,  and  (2)  variations  in  the  local  angles  of  attest. 

In  the  transonic  range,  the  experimental  result*  for  she 
65-  and  75-degree  mot  els  are  essentially  the  same.  For  th*  lam* 
supersonic  range,  the  level  of  the  ^/i,u curve  i*  lower 
for  the  75*degree  models.  Also  the  rapid  rise  in  the  curve, 
which  is  evident  for  all  models,  occurs  at  a  higher  Mach  mother 
for  the  75-degree  models.  These  results  may  be  partially  ds*  to 
the  higher  frequency  ratios  (  CJ,  /<, Jz  )  obtained  with  the  75- 
degree  models.  Eecause  of  the  substantially  higher  frequences 
of  the  7S-degree  models,  it  was  possible  to  obtain  f bitter  points 
in  the  Kach-number  range  1.2  to  2.5.  For  the  65>degree  models, 
flutter  points  could  not  be  obtained  in  this  Mach  number  region,, 
the  models  being  unstable  even  at  the  lowest  tunnel  dynamic  pres¬ 
sures. 

The  experimental  results  are  quite  consistent,  considering 
the  variations  encountered  between  some  of  the  models.  All  of 
the  supersonic  results  (for  M  ^  2.5)  lie  in  a  band  of  ±8  per  cimtr 
about  «  median  experimental  line  (see  Fig.  3. IS),  The  region  from* 
M  -  1.2  to  2.5  is  not  too  well  covered,  but  the  trend  for  t*m 
65-degree  models  is  reasonably  well  defined.  It  is  also  clsm X 
that  the  75-degree  boundary  lies  below  that  for  the  65-degzmm 
models  over  the  Mach-number  range  from,  say,  1.2  to  3.0. 
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APPENDIX  I 

ANALYTICAL  PROCEDURES 

A.  AN  APPROXIMATE  METHOD  OF  DEFLECTION  AND 
VIBRATION  ANALYSIS 

A  simplified  "strain-energy  approach"  is  derived 
here  to  compute  the  deflection  and  frequency  characteristics 
of  inflatable  plates.  It  is  specialized  to  the  case  of  a 
plate  of  uniform  thickness,  undergoing  small  deformations, 
i.e.,  the  linearized  formulation.  It  is  considerably  simpler 
than  the  more  general  linearized  approach  of  Ref.  1.1.  The 
agreement  between  the  two  methods  is  discussed  at  various 
stages  of  the  derivation.  For  this  purpose,  the  notation  of 
Ref.  1.1  is  adopted  insofar  as  possible. 


-ton <r 


Delta 

Planform 


Ml* 


,  Drop  Chord 
^  Small  Thickness 


Figure  X.l.  Coordinate  Systems  and  Deflection 
Rotations  for  Sisplified  Analysis* 
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Consider  Fig.  1.1,  which  depicts  elements  of  the 
plate  as  seen  in  the  x-z  and  y-z  planes.  Within  the  frame¬ 
work  of  small  deflection  theory  with  shear  deformations 
present,  and  assuming  small  uniform  thickness  ,  the  dis¬ 
placements  na7  he  expressed  as 

ic  m  Zee  V  *  Zj3  and  >/*  * 

The  strain-displacement  relations  follow; 


C  =  -  Zoe 

"  dx  *  * 


6 

n 


and 


,  dtr 
9y  d* 


s  <kc 
'■*  ~d* 


2x 


dtr  •  2*aT 
9a  Of 


The  corresponding  stress-strain  relations  for  an  isotropic 
plate  are  given  by 


-// 


*« 


/3 


J—  (T 

<S  *3 


a  6f  *3 


a- 3a) 


or  inversely. 


*  Subscripts  X  ,  H  on  oc ,  fi  ,  indicate  partial  dif¬ 
ferentiation  with  respect  to  pc»  H  • 
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<5  -7 rp-  'S  “  i fbr&*V‘ V 


0“  .*  /CiC  js 


Z('+/^ 


&  £T*  a  ,  A  fid  0“  »(^ 

C/Z  4*  >3  23  ^3  (I.3b) 


Certain  modifications  must  be  made  to  the  stress-strain  re¬ 
lations  (Eq.  (1.3))  before  applying  them  to  the  inflated, 
plate-type  wing.  Firstly,  the  shear  modulus  G  in  the 
relations 


^7a 


and  0J3 


must  be  replaced  by  ^  .  This  can  be  easily  demonstrated. 
For  exanple,  from  the  sketch  below,  one  has  for  the  shear 
stress  (T^  (for  small £•>„+«)), 

<S  flr  t1-4*) 


Similarly, 


(I. 4b) 


[W,+a 


■P<h  (*,«.<»} 


•  DIRECTION  Or  OROP  CHOROS 


.  •  .A 
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On  this  basis,  and  assuming  the  covers  to  behave  iso¬ 
tropically,  the  strain  energy  is  given  by 

U  *4-/7/  J<r£+<T6  e+r<eia+ ere  IdxchJz.  (1.5) 

s  2. JJJ  1  H  U  Zi  ZZ  tz  '*  4 i  /3  ij  Z3J  a 

Secondly,  let  the  stress-strain  relations  be  expressed  in 
the  more  general  fora 


*5*4v*A*«  Ajzz  (i,6> 

which  may  be  applied  to  orthotropic  as  well  as  isotropic 
situations. 

Substitution  of  Eqs.  (I.4a,b)  and  (1.6)  in  Eq.  (1.5) 

yields 

{a»£'+  44 

+  i;jfj ' (i.7) 

The  shear  stresses  and  &2s  are  implicitly  assumed  to 

be  uniform  over  the  depth  A. ,  so  that  the  z-integration  f6r 
the  second  of  the  two  triple  integrals  in  (1.7)  yields 

zj <i  a> 

.  *• 

On  the' other  hand,  stresses  <TJj  ,  0^  ,  and  0^x  are 

carried  by  the  covers,  so  that  the  appropriate  z-integration 
is  carried  out  as  follows: 

+JM  ur-fr]  (i*9> 

-5./  i 
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Utilizing  strain-dispincccent  relations  in  (1.7)  and 
performing  the  z-integrations  as  indicated  by  Eqs.  (1.8) 
and  (1.9),  one  obtains 


U  *  £Jj *  4 //  4>„  $  *  4 

+ ijjh1  {(*+ *  ft*  (i. io) 

3 

where  the  double  integrals  are  taken  over  the  planform  area 

5. 

The  latex-impregnated  cloth  covers  do  not  behave 
isotropically;  if  the  covers  are  treated  as  orthotropic, 
and  the  x-  and  y-axes  are  aligned  with  the  principal  di¬ 
rections,  it  follows  that 

* 

a  Eft  n  _ /?  _  A;  5*  '/**.&*£ 

u  '-Wu  *.  "  /vWfa  * 


/-/*«/** 


a*J 


<x,u> 


The  structural  parameters  £t£r  and  JEtt  are  the  ex- 
tensicmal  stiffnesses  in  the  x-  and  y-directions,  re¬ 
spectively.  is  the  Poisson's  ratio  associated 

with  a  contraction  in  the  y-  (or  x-)  direction  caused  by 
a  tensile  stress  in  the  x-  (or  y-)  direction.  The  values 
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of  these  struccural  constants  depend  on  the  stress  levels  of 
the  median  state,  the  state  about  which  small  additional 
deflections  are  assumed  to  take  place.  Let  denote  the 
work  done  b j  the  external  forces,  which  are  assumed  to  act 
normal  to  the  x-y  plane;  then 

W»  JJ  f 'h*' (1*12) 

being  the  associated  load  intensity.  According  to  the 
principle  of  minimum  potential  energy,  one  oust  have  zero 
variation  of  ,  i.e., 

SIT  =-  S(  U5-’W)^0  (1.13) 

With  appropriate  integrations  by  parts,  and  from 
variations  on  *7,  oc  ,  ,  Eq.  (1.13)  yields  the  differential 

equations 

PA  (.*+ P\  — / 

r/4  *=*4  &  ’■4  («„%)] 

t[A*  tl-lA) 

These  are  in  agreement  with  Eqs.  (42),  Ref.  1.1,  for  the 
case  of  time-independent  deflections.  When  motion  is  in¬ 
cluded,  one  Mist  add  to  the  right-hand  sides  of  Eqs.  (I.1A), 
the  inertia  terms  77700^  -5 n  tnd.  '0A!‘A  ,  respectively. 
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Here,  ?n  represents  the  mass  of  the  plate  for  unit  area 
(top  and  bottom  covers  included)  .♦  From  the  variational 
principle,  the  natural  boundary  conditions  associated  with 
Eqs.  (1.14)  are,  for  a  rectangular  plate, 

(1)  along  an  edge  where  X  "  constant 

either  W'  »o  or  oC^O 

either  oC  •  O  or  At/  Aify*0 

either  ft  m  O  or  m  & 


(2)  along  an  edge  where  ^  »  constant 
either  mO  or  %K+/i 


either  ■» J  •■O 
either  c(  »0 
either  /3~0 


or  **  O 

ot  ex„+fit~o 
or  Au^*A,^0 


These  boundary  conditions  assume  that  no  edge  forces 
are  being  applied;  they  are  slightly  different  from  those 
given  in  Eqs.  (39)  of  Ref.X.l.  This  apparent  discrepancy 
between  the  two  analyses  arises  from  a  term 


aU 


which  is  Included  in  laf.  1.1,  but  which  is  missing  from 
the  present  derivation.'  It  is  worth  mentioning  that  this 
AU  -term  does  not  change  the  system  of  Eqs.  (1.14). 
Furthermore,  this  difference  can  be  shown  to  lead  to 
practically  negligible  terms  in  the  natural  boundary  con¬ 
dition. 

For  the  present  application  to  a  delta  plan  form, 
the  method  of  solution  described  for  a  rectangular  plate 
is  inconvenient.  Instead  of  solving  the  differential 


in  corresponds  to  the  quantity  in 'Ref.  1.1. 


131 


FDL  TDR  64-147 


equations  subject  to  certain  boundary  conditions,  one  can 
return  to  the  original  energy  equations  and  use  a  Rayleigh- 
Ritz  approach.  A  particularly  simple  case  is  one  where  the 
bending  stiffnesses  are  assumed  so  large  that  the  deforma¬ 
tions  are  purely  of  the  shear  type.  For  such  a  situation, 
one  has 


cc  0  £5  =  Jj  &  fax*  rfjdxe/y  (1 . 15) 


To  facilitate  the  ensuing  analysis,  the  dimensionless 
coordinates  X,  £  are  introduced,  where 

and  dT  is  the  root  chord.  Also,  let 

O 

*(*,?,*)- jh  t.VKft?)  >  <I17> 

l»*f 

where  are  trial  mode  shapes  (or  modal  functions), 

suitable  in  the  sense  that  each  satisfies  the  geometric 
boundary  conditions.  From  Eqs.  (1.15)  -  (1.17),  on*  has 

•9‘toj <*■»> 

For  simplicity,  one  may  consider  "T^'s  which  are 
single  terms  in  powers  of  %  and  .  For  the  exaople 
presented  here,  the  following  power  foms  sre  taken: 
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m. 


ire 


The  double  integrals  of  7rf  •  and  ^  *  may  be  easily 

•  /|  */• 


evaluated,  using  the  relation 

•<rj 


where  0~  is  the  tangent  of  the  semi-apex  angle.  Equations 
(1.20)  form  a  . set  of  homogeneous  algebraic  equations  wich 
the  eigenvalues  '  Jl()  •  •  •  ,  ^  ,  and  may  be  solved  by 
any  one  of  several  well-known  techniques,  by  iteration  for. 
instance.  To  each  eigenvalue  Xj  is  associated  an  eigen¬ 
vector  gC i )  with  components 

eigenvector  defines  the  mode  shape.  Since 

and  £ .  are  functions  of  r  only,  X:'*  are  fixed  once  the 
semi-apex  angle  is  specified.  For  the  65-degree-sweep 
deltas*  considered  here,  the  first  three  eigenvalues  turn 
out  to  be: 


-  / 9.Z2V 


Xz- 77.5*9  ,  -  /f/32- 


with  the  associated  eigenvectors 


More  precisely  these  results  apply  to  model  1,  which  bn 
a  leading  edge  very  slightly  different  from  65  degrees* 

All  "65-deg  models"  have  sweeps  sufficiently  near  that  Cor 
Model  1,  so  that  these  results  apply  to  all  of  them. 
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-0.011,  713 

0.034,  604 

-0.037,  22 

‘0.106,  14 

-0.312,  80 

-0.532,  72 

-0.197,  51 

-0.167,  02 

0.641,  80 

0.052,  988 

U) 

0.303,  64 

0.962,  41 

1.000,  00 

$  H 

1.000,  00 

t  ^ 

1.000,  00 

-0.873,  19 

-0.253,  30 

-2.942,  5 

-0.906,  36 

-0.964,  02 

-2.713,  9 

0.772,  15 

0.746,  00 

0.807,  94 

0.031,  149 

-0.566,  77 

4.851,  1 

L  J 

. 

The  natural  frequencies,  for  a  particular  model,  can 
be  determined  by  inserting  the  above  eigenvalues  and  the 
model  parameters  into  Eq.  (Z.21).  The  locations  of  the  node 
lines,  which  are  the  same  for  all  models,  can  be  calculated 
by  inserting  the  assumed  modal  functions  and  the  associated 
eigenvectors,  given  above,  into  Eq.  (1.17),  and  then  setting 
the  left-hand  side  equal  to  zero  and  solving  for  jx*  as  a 
function  of  ^  (  or  y  as  a  function  of  %). 

The  deflection  of  the  wing  models  under  a  uniform, 
static  load,  again  based  on  the  assumption  of  pure  shear- 
type  deformations  only,  can  be  calculated  from  the  non- 
homogeneous  Lagrange's  equation,  which  takes  the  form 


jw 


&  2  tA  ■  $  -  two 


dot*  -/-V.T 


where  represents  the  1-th  generalized  external  force 
and  y  is  the  load  intensity.  The  system  of  non  homogeneous 
equations  may  be  solved  to  yield  the  generalized  displace¬ 
ments  J  's,  and  from  these,  the  deflections.  For  the  65- 
degree-sweep  delta  models,  the  deflection  &  turns  out  te 
he 
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ttfXtf)  *  ~~  j~  0. 00^3?  +  0. 6  fo  Ox'y-0.  0. 3/3 f  X*y 

+  Q.3J7J  *f*-  O.  A*4{  Zjf*-  0. 33^3  Z*fX 

+  A3«>f yS-  o.a/7,3?fJ  (1.23) 

B.  FLUTTER  ANALYSES 

1.  Equations  of  Motion 

The  analytical  flutter  boundaries  presented  in  this 
report  are  based  on  the  expression 

Hn  i+  2  V  i{  *  tfO  Hr*  Cl-»> 

derivable  from  Lagrange's  equation  (cf.  Ref.  1.2  t  page  407)  i 
Here, 

»  generalized  coordinate  of  the  nth  orthogonal 
mode. 


nth  generalized 


ate  the  out-of-phase  and  in-phase  aerodynamic 
f  ^  coefficients  which  will  be  discussed  below,  and 

(j)  is  the  ncJl  orthogonal  mode  shape. 

Parallel  static  divergence  boundaries  are  also  given;  they 
are  based  on  the  reduced  fora  of  Eq.  (1.24) 
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i 


» 


t 


\ 


* 


* 


H 

Mn  A  *2  A;  t  m  0  (I  • 25) 

n  *  \  y«/  7  7 


The  structural  damping  nay  be  Included  in  the  flutter 
analyses  by  applying  the  factor  (  L-fig  )  to  the  term 

•  where  ^  is  the  structural  damping  coefficient 
of  the  ntj!  mode.  The  coordinate  system  is  -that  given  in 
Fig.  1.2.  In  using  Eq.  (1.24),  the  following  are  implicitly 
assumed: 

(1)  The  flutter  or  divergence  mode  is  representable 
by  a  linear  combination  of  a  finite  number  of 
normal  (orthogonal)  modes.  In  the  applications, 
combinatiorfs  of  two  and  three  degrees  of  free¬ 
dom  are  considered. 

(2)  For  the  supersonic  flutter  calculations,  the 
generalized  aerodynamic  forces  are  quasi¬ 
steady;  they  may  be  expressed  in  terms  of  two 
components,  one  in  phase  with  the  displacement, 
the  other  in  phase  with  the  velocity.  For  a 
flat  plate,  and  according  to  first-  or  second- 
order  piston  theory. 


J. 

7 


(1.26) 


% 


did? 


where  f  ,  V  ,  and  M.  are  the  free-streaa  density, 
velocity  and  Mach  number,  respectively.  For  the  transonic 
calculations,  the  combination  of  t« 


137 


y* 


FBL  TDR  64-147 


N 

21 

i-t 


9 


aT 

-2  A,.  5  • 


trust  be  replaced  by  an  appropriate  aerodynamic  generalized 
force  term.  This  point  is  further  discussed  below. 


2.  Mode  Shapes  Used  in  the  Calculations 


In  all  of  the  analyses  presented  here,  one  of  the  following 
sets  of  mode  shapes  is  utilized.  All  are  based  on  experimental 
vibration  results  for  a  model  pressure  differential  of  2  psig. 


Set  A 

These  are  analytical  modes  derived  from  the  vibra¬ 
tion  data  of  Model  1.  The  orthogonal  modes  are  approximated 
by  polynomials  of  the  form: 

#  (1.27)* 

jg  •  (i+ax  +  bx^jf'+cy* 


The  constants  a,StC  in  the  assumed  mode  shapes 
are  determined  by  imposing  the  orthogonality  condition 

JJ  dxdj  *0  (1.28) 

3 

and  by  satisfying  the  node  line  locations  obtained  from  the 
Model  1  vibration  data.  If  the  mass  distribution  fH  is 
assumed  uniform,  Eq.  (1.28)  reduces  to  ' 

II  ?!  ^ *  0  <I*28a> 

3  . 

Forcing  the  node  line  of  the  second  mode  to  pass 

Superscripts  like  A  in  (p^  designate  the  sets  of  mode 
shapes.  ' 
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through  the  leading-edge  point  (  X  »  0.83,  %  m  0.38)  and 
to  approach  the  root  at  the  point  (  *  »  0.9 lS,  $  m  0)  re¬ 
sults  in  the  following  node  shapes 


(l- l‘lS7.IX  +  &f*‘(S‘XZJjfi l-0./8l,9f*  (1.29) 

The  choice  of  the  above  two  nodal  points  is  arbitrary. 
The  question  arises  as  to  how  well  Eqs.  (1.29)  represent  the 
actual  vibration  modes.  This  can  be  answered  only  by 
measuring  the  complete  vibration  mode  shapes.  For  two  of 
the  models,  the  mirror  reflection  technique  was  used  to 
obtain  experimental  deflection  slopes;  from  the  latter,  the 
mode  shapes  can  be  calculated,  see  Sets  C-l,  C-2  below.  A 
qualitative  confirmation  of  the  assumed  modes  was  attempted 
by  visually  observing  constant  deflection  lines,  using  the 
NASA  acceleration  technique  (cf.  Ref.  1.3)  and  comparing  these 
contour  lines  with  those  obtained  from  the  assumed  modes. 

The  theoretical  contour  lines  for  the  first  mode,  Ty* 
are  given  in  Fig.  1.3;  these  agree  with  the  visual  obser¬ 
vations  of  "l-g"- lines  as  the  amplitude  of  motion  is  varied. > 
Confirmation  of  the  second  and  particularly  higher  modes  is 
more  difficult,  and  other  modes  can  be  assumed  which  would 
agree  reasonably  well  with  experimental  l-g- lines.  Contour 
lines  for  the  second  mode  are  given  in  Fig.  1.4. 

Set  B 

This  set  is  also  comprised  of  assumed  analytical 
mode  shapes.  The  first  mode  is  identical  with  that  of 
Set  A.  The  second  mode  is  assumed  to  be  of  the  fora 

fz  ~  (1.30) 
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which  involves  no  chordwise  bending.  Retaining  the  first 
mode,  and  imposing  exactly  the  same  conditions  as  in  Set  A, 
one  arrives  at  the  second  mode  shape. 

^B>-  O-i.  ojz, ;  ?)/*+(-  *.  o*,9+  4-  059/ xjf*  (1 . 31) 

The  contour  lines  of  are  given  in  Fig.  1.5. 

Comparison  of  Figs.  1.4  and  1.5  shows  that  the  contour 
lines  for  cf>^  and  are  similar  in  nature, 

although  exhibits  much  larger  relative  deflections 

at  the  tip  and  at  the  middle  of  the  leading  edge.  For  this 
set,  an  expression  for  the  third  mode  in  the  polynomial 
form 

(/+4Z  +  j>ZZ)fZ  +(c+c/x)tfA  (1.32) 

can  also  be  determined.  The  constants  a,  *  *  •  d.  are  de-  • 
termined  by  imposing  two  orthogonality  relations 

s  s 

(a) 

and  by  making  the  node  lines  of  pass  through  the 

points  (  %  ■  1,  «  0.33)  and  (  X  -  0.5,  ^  m  Q).  Once 

more,  these  points  are  taken  from  the  experimental  node 
line  for  Model  1  and  their  choice  is  arbitrary.  This  yields, 

/o\ 

<f3  =*  (/- 2  30S3x  +  /. 3/06  X (3-&9.7~ZA7’l/)L)/i  (1.33) 

The  contour  lines  for  the  third  mode  are  given  in 
Fig.  1.6. 
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Sets  C-l  and  C-2 

These  sets  include  three  mode  shapes  derived  from 
the  mirror  tests  on  Model  8.  Subsection  D  of  Section  2 
describes  the  experimental  techniques  and  the  data  reduc¬ 
tion  procedures  for  obtaining  deflections  and  streaowise 

slopes.  cP^  and  are  obtained  for  some  55 

7  n  "3T« 

selected  points  on  the  planform.*  For  Sets  3,  where 

the  mode  shapes  are  expressed  by  ar  lytical  t  is, 

the  mode  shape  integrals 

IJtfjW*  ■  ffr.% 

may  be  obtained  in  closed  form.  For  Sets  C,  inasns 
the  data  are  furnished  at  a  large  number  of  discrete  pox. 
numerical  integration  techniques  are  more  convenient.  In 
all  cases,  the  rectangular  rule  is  applied.  The  generalized 
masses  are 


"  c*Jf *>(*,*)  fluff'd*-*? 


(1.34) 


where  the  index  p  denotes  the  points  1,  tl,  •  •  •  ,  37. 

The  weighting  factor  Tv^  is  the  mass  of  that  portion  of 
the  wing  assigned  to  point  p  ,  and 

The  mass  distribution  is  assumed  uniform,  i.e. ,  jjj  ■  con¬ 
stant,  for  all  areas  nonadjacent  to  the  leading  and  trailing 
edges.  The  overlap  between  top  and  bottom  covers  for  areas 
adjacent  to  the  leading  and  trailing  edges  is  accounted  for 

These  include  the  mirror  points,  augmented  by  some  inter¬ 
mediate  points;  they  do  not  include  the  points  1R,  9R,  .... 

37R  on  Fig.  2.11.  ’ 
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in  the  mass  distribution. 

Sets  C-l  and  C-2  are  identical  except  for  the 
generalized  aerodynamic  forces.  For  the  aerodynamic  in¬ 
tegrals,  a  weighting  factor  Ap  »  which  is  the  effective 
area  of  the  wing  assigned  to  the  point  p  ,  may  be  intro¬ 
duced.  For  exaaple. 


(1.35) 


In  general,  the  Aj,  is  simply  the  midplane  area  of 
the  segment  of  the  wing  assigned  to  the  point  p  .  For 
Set  C-l,  the  rounded  portions  of  the  leading-  and  trailing- 
edge  regions  associated  with  the  point  p  are  included  in 
the  calculation  of  Ap  .  For  Set  C-2,  these  rounded  por¬ 
tions  are  not  included  (in  the  calculation  of  Ap  .  That  is, 
the  Ap  is  the  total  midplane  area  minus  the  midplane  area 
of  the  rounded  edges  (see  Fig.  1.2).  The  justification 
for  considering  this  modified  Set  C-2  is  discussed  below. 


3.  Aerodynamic  Forces 


For  the  Mach-number  range  from  2.5  to  6.0,  piston 
theory  with  and  without  modifications  is  used  to  obtain 
the  generalized  aerodynamic  forces.  For  M  ■  1,  Landahl's 
.first-order  theory  for  low  aspect  ratio  delta  wings  is 
applied. 

Unless  stated  otherwise,  the  models  are  assumed  to 
be  set  at  zero  angle  of  attack,  and  are  thus  undistorted. 
Since  the  model  is  of  uniform  thickness  over  most  of  the 
planform,  it  is  reasonable  to  expect  that  calculations 
based  on  a  flat-plate  assuc^>tion  should  give  adequate  esti¬ 
mates  for  the  aerodynamic  forces.  The  errors  committed  by 
such  an  assumption  (at  least  for  these  slender  configurations 


'k 
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at;  moderate  Mach  numbers)  are  confined  to  the  leading-  and 
trailing-edge  regions  (see  Fig.  1.2). 

As  a  simplification,  one  may  arbitrarily  neglect  the 
contributions  of  these  blunt  regions.  Although  this  may 
appear  quite  empirical,  one  can  offer  a  qualitative  physical 
justification  for  the  two-dimensional  case  with  semicircular 
edges.  Consider  the  section  sketched  below.  The  contributions 


of  the  leading-  and  trailing-edge  shaded  regions  to  the  air¬ 
loads  associated  with  motion  come  only  from  the  "translation" 
part.  Any  rotation  about  (or  )  would  give  zero  total 
lateral  load  an  the  shaded  region  at  jPf  (or  J?z  )  for  very 
small  deflections.  Consequently,  the  contributions  from 
the  shaded  regions  are  of  order  -k.  ("  /V  »  the  re¬ 
duced  frequency).  This  means  neglect  of  these  regions  re¬ 
sults  in  an  error  in  the  "  jik. -terms"  of  the  flutter  de¬ 
terminant  elements.  At  the  high  Mach  numbers,  the  flutter 
is  very  small,  (at  least  for  flutter  involving*  the  first 
two  or  three  modes  predominantly),  and  the  resultant  ca¬ 
tenas  should  mot  influence  the  flutter  results  appreciably  * 

An  alternate  procedure,  which  is  also  somewhat  em¬ 
pirical,  is  based  on  the  use  of  the  simple  Newtonian 
theory  for  that  portion  of  the  airfoil  surface  near  the 
leading  edge  where  the  sine  of  the  flow  inclination  9 
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trices  M  is  greater  than  unity,  i.c..  Main 8  > /  •  For 
the  rest  of  the  airfoil  surface,  the  piston  theory  *s  ap¬ 
plied.  This  cct-hod  is  suggested  by  Morgan,  Runyan  and 
Ihickcl  (Ref. I. 4).  The  slight  improvement^ afforded  by  the 
latter  method  over  the  forcer  which  neglects  entirely  these 
regions,  is  not  significant  here,  for  there  are  other  phases 
in  the  dynamic  analyses  which  arc  equally  approximate. 

The  above  arguments  apply  to  the  two-diaonsional 
case  and  thus  for  zero  weep.  Since  one  would  expect  a 
dcci-easing  influence  from  bluntncss  with  increasing  leading- 
edge  sweep,  the  suggested  corrections  are  excessive.  'In 
the  analyses  presented  here  for  large  sweep  delta  models, 
one  is  perhaps  justified  in  neglecting  bluntncss  effects 
altogether,  and  in  treating  the  wings  as  flat  plate3  of  the 
same  planforms.  A  suggested  method  whereby  bluntncs3  and 
initial  angle  of  attack  could  be  included  in  the  analyses 
is  presented  below. 

For  deriving  the  unsteady  airloads  at  the  sonic 
speed,  Landahl's  first-order  theory  for  low  aspect  ratio  delta 
wings  is  utilized.  Davies  (Ref.  1.5)  presents,  in  detail^ 
derived  results  for  the  airloads  associated  with  some 
elementary  mode  shapes  of  oscillation.  However,  the  mode 
shapes  considered  there  arc  insufficient  for  the  present 
purpose.  Generalized  aerodynamic  forces  are  derived  for 
additional  modes.  In  particular,  the  required  mode  shapes 
which  constitute  the  elements  of  the  first  and  second 
(assumed)  mode  shapes  of  Model  1  are  (sec  Eq.  (1.27)) 

Zr  </  -  4,5,6,7.8) 

They  are  designated  as  elementary  modes  4,  5,  6,  7,  8, 
respectively.  No  attempt  is  made  here  to  reproduce  tha 
lengthy  derivations,  which  essentially  parallel  those  for 
the  simpler  modes  in  Ref. 1. 5.  However,  the  dimensionless 
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(D  * 


velocity  potentials  £P j  ,  their  streacwise  derivatives 

t ,  and  the  associated  generalized  forces  are  sum¬ 
marized  below.  The  notations  are  those  of  Fig.  1.2  and  is 
is  the  reduced  frequency  • 


Mode  4: 


Zii  *■  (£) 

a^5s4,  \  ‘2.  J 


Mode  S: 


¥  0,(*,  -('+  Zt*  *)•/&$ rLfr  J 


*  -jSw  (4)  -  *  zW^*y^ 


VaijT2  J  L 


Mode  6: 


•  f4 '&***  -  -  ij,  *_  r 


Mode  7: 


-  -  *■<  **¥+* **j 

g  ■■  "  ■“ 

Superscript  (1)  indicates  first-order  solution,  which  is 
essentially  the  slender  body  result. 


US 
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Mode  8: 


?7v>j  —  c-^i2jzr'[Y^Y*Js‘,l 


-  -  zik  ~L - ./Is  7^  rf*5**- s VZy‘/s’3l/ 

tJZJ*  3  ‘  * 


The  generalized  aerodynamic  forces  are  defined  as 


Q'r$?  -  7^  ffz.wfa/*-yP*jr 


(1.36) 


where  Af>j(x,y)  is  the  pressure  difference  between 
bottom  and  top  surfaces  of  the  wing  undergoing  an  oscil¬ 
lation  in  mode  j  .* 

<c.  *>[-*+#}  eC-g[-,-t»+& 

<£•  it  *  ry  €- &  t  i‘j 


*The  pressure  difference  is  related  to  "f  in  the 
manner  indicated  by  Eq.  (1.36),  see  Ref. 1.5. 
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# 


; 


C- 

Q‘2  -  127/1 2-  i*Sl 

^  ¥8  L  *  9  S  J 

4 ?-§£[- 1**#] 

Q2-gH-**+  u-f] 

Q2-gH- *r+  Sr\ 

Q8-&L-¥**PJ 

n<»_  3ir**r  it  zoik.  ixi 
»*  -  zsl  L~S  //  4  TJ 

*G-*£b-9»+*trl 

*£1 

*£-&{■*- v*+t 7 

• 

• 
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in 

Appropriate  linear  combinations  of  the  's 

yield  tiie  generalized  aerodynamic  forces  associated  with 
the  assumed  modes  in  the  transonic  flutter  analyses.  For 
example,  if  the  indices  1  and  2  denote  the  first  and 
second  orthogonal  vibration  modes  of  Set  A  (see  Eq.  (1.29)), 

,  c0> 

'  5S 


Ql 


=  Qsi  ~  '■  ?SV  +  0.9995 -0.W9  Qsl 
Cf*  =  0%  -  >■  9sy  Qsl  +  0.99%  sQ^'J  -  0.  /s/,9  Q™ 

€  =  #VY  +(-,  ttSvfQss  i-(-°-,8)rfQn 

-/  95%,  (*2+QZ)t*9its(Q£+Q2)  -o.m.*(fyoZ) 

+(->9S7j)(o.  W.s)[Ql‘l+Q%]+(ri-  9sI/X-o*V)(Qs7+Q(v) 


+  (o.  999,  s)(-  0.  tn,  9)  [o™  +  q£J 


It  can  be  easily  proved  that,  when  using  the  tran¬ 
sonic  theory,  it  is  appropriate  to  replace  the  and 

~omj  -  terms  in  the  equations  of  motion  (1.24)  according 
to  the  rule 

» 

(\irhJ,vt/)  b  fa-'-*) 

4.  Summary  of  Calculations 

Table  1.1  presents  a  summary  and  index  of  the  tvo- 
and  three-degree-of-freedom  flutter  and  divergence  analyses 
that  have  been  carried  out.  The  results  are  plotted  in 
Figs.  1.7  -  1.24  and  consist  primarily  of  the  stability 
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boundaries  ^f/bu )7  (or  Vs/bm  *7*  for  the  cases  of  static  diver¬ 
gence)  vs  (or  }L  for  the  transonic  calculations,  M  ■  l)./c 
is  the  familiar  mass  ratio  parameter,  defined  conveniently  for 
the  present  analysis  in  terms  of  the  total  mass  of  the  wing,  a, 
as 

- - iSL. - - 

'  2.  Kfjt/tr 

The  choice  of  these  parameters  as  ordinate  and  abscissa,  re¬ 
spectively,  is  natural,  in  that  it  is  in  this  form  that  the  theo¬ 
retical  results  are  obtained  when  one  uses  piston  theory  or 
s lender-body  theory.  When  comparing  theory  and  experiment,  re¬ 
sults  are  reinterpreted  in  the  more  common  form  vs  H. 

In  most  instances,  dimensionless  flutter  frequencies,^/^  are 
also  plotted  vs 

A  portion  of  these  parametric  studies  were  performed  early 
in  the  program;  specifically,  the  curves  utilizing  mode  shape 
sets  A  and  B  were  obtained  once  the  vibration  and  inertia  proper¬ 
ties  of  Model  1  became  available.  These  results  were  used  in  the 
design  phases  of  the  other  models.  The  remainder  of  these  cal¬ 
culations  were  carried  out  at  later  times,  for  purposes  of  cor¬ 
relation  between  theoretical  and  experimental  flutter  results. 

For  use  in  the  interpretation  of  the  experimental  data  (Sec¬ 
tion  3),  the  important  points  brought  out  by  these  calculations 
are  now  briefly  discussed.  4 

Figures  1.7  and  1.8  establish  flutter  and  static  divergence 
boundaries  for  the  typical  65-degree  configuration.  Two  flutter 
branches  are  shown  in  each  figure,  corresponding  to  two  possible 
solutions  from  the  flutter  determinant.  Only  the  lower  branch  is 
of  interest  here.*  The  higher  speed  boundary  for  static  diver¬ 
gence  means  that  the  model  will  flutter  before  the  static  diver¬ 
gence  condition  is  reached.  This  is  true  in  all  cases  carried 

*In  the  subsequent  Figs.  1.13  to  1.24,  only  this  lower  branch  is 
shown  whenever  two  stability  boundaries  are  obtained. 
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out,  with  some  predicting  no  divergence  at  any  speed.  These 
plots  also  indicate  that  inclusion  of  equal  structural  danping 
is  slightly  stabilizing,  and  results  in  higher  flutter  fre¬ 
quencies.  Similar  conclusions  are  reached  when  mode  shape  sec  B 
is  substituted  in  the  above  calculations,  as  evidenced  by  the 
curves  of  Figs.  1. 11  and  1.12.  Figures  1.9  and  1. 10  are  included 
to  extend  the  results  with  mode  shape  set  A  to  the  sonic  Kach 
number.  Figures  1.13  and  1.14  present  the  results  according  to 
a  three-degree-of-freedom  analysis.  The  interesting  points  about 
these  curves  are  (1)  the  lack  of  a  flutter  boundary  near  that 
associated  with  modes  1  and  2,  and  (2)  the  weak  dependence  of  the 
higher  boundary  on  mode  1,  i.e.,  modes  1,  2,  3  boundary  is  nearly 
the  same  as  that  for  modes  2,  3  alone.  A  close  examination, 
using  the  usual  "V-g"  type  flutter  analysis,  revealed  that  the 
"disappearing  branch,"  with  flutter  predominantly  made  up  of  the 
first  two  modes,  is  very  lightly  damped.  This  means  that  in  the 
analysis,  the  coupling  between  modes  2  and  3  are  seriously  over¬ 
estimated,  a  strong  possibility  considering  the  extent  of  the 
approximations  involved  in  the  mode  shapes. 

For  Figs.  1.15  to  1.24,  the  calculations  are  based  on  the 
mode  shapes  C-l  and  C-2.  Figure  1.15  shows  the  influence  of 
lowering  the  frequency  ratio,  U),/u>z  ,  to  be  a  stabilizing  one, 
but  rather  small  for  the  range  0.4-^  A  0.6.  Figure  1.16 

presents  the  influence.on  the  flutter  boundary  of  the  bluntness 
of  the  leading  and. trailing  edges,  following  the  approximate 
technique  of  Subsecticn  3  of  the  present  section.  This  effect  is 
also  slightly  stabilizing.  The  influence  of  the  frequency  ratio, 
/u>z  ,  on  the  higher  mode  flutter  (modes  2  and  3)  is  brought 
out  in  Fig.  1.17.  It  is  interesting  to  note  that  higher  *i/<*>a's 
are  stabilizing.  A  rapid  lowering  of  the  flutter  boundary  is 
noted  as  (u)z  approaches  unity.  At  the  value  of 
this  boundary  is  seen  to  be  below  that  for  flutter  involving  the 
first  two  modes,  when  one  compares  it  with  Fig.  1.7.  The  influ- 
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ence  of  the  edge  bluntness  is  seen  (Fig.  1.18)  to  be  stabilizing 
on  the  high-frequency  flutter,  as  in  the  previous  ease. 

Figures  1.19  and  1.20  show  that,  with  mode  shape  set  C-l, 
small  equal  structural  damping  coefficients  for  the  participating 
modes  have  slight  destabilizing  effects.  This  destabilizing  in¬ 
fluence  is  considerably  more  with  unequal  damping  coefficients  as 
may  be  noted  from  Figs.  1.21  and  1.22.  Finally,  Figs.  1.23  and 
1.24  present  the  results  of  three-degree-of-freedoa  analyses 
utilizing  the  mode  shapes  derived  from  the  mirror  tests.  These 
results  lead  to  the  same  conclusions  on  the  effects  of  edge  blunt¬ 
ness  and  of  structural  damping  as  those  from  the  two-degree-of- 
freedom  cases  of  Figs.  1.16  and  1.18  to  1.20. 

C.  FLUTTER  OF  BLUNT  WINGS  AT  LARGE  INITIAL  ANGLES  OF 
ATTACK 

1.  General 

It  is  desirable  to  assess  the  influence  of  the 
initial  angle  of  attack  on  the  flutter  index,  as  in¬ 
flatable  vehicles  are  expected  to  fly  at  large  cx^'s 
during  reentry. 

*  In  any  reliable  analysis,  two  major  modifications 
oust  he  introduced: 

(1)  The  changes  in  the  vibration  characteristics 
from  those  in  still  air.  It  is  expected  that 
the  large  deformations  due  to  static  lateral 
loading  alter  significantly  the  frequencies 
and  mode  shapes. 

(2)  Piston  theory  becomes  increasingly  inaccurate 
ss  Afc/Viot  appi  caches  the  values  of 
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1.*  For  M  »  5  or  6,  this  limit  is  reached  with  at. 
10  deg,  even  without  elastic  deformations. 

Beyond  this  point,  piston  theory  must  be  re¬ 
placed  with  another  aerodynamic  theory,  such 
as  the  simple  Newtonian.  The  latter  improves 
in  accuracy  with  increasing  M.sina(  . 

In  the  initial  phases  of  the  program,  some  analytical 
work  was  carried  out  on  the  simultaneous  effects  of  oC9  and 
bluntness  on  the  aerodynamic  forces  according  to  Newtonian 
theory.  However,  no  attempt  was  made  to  extend  and  apply 
the  results  of  this  study  to  particular  model  configurations, 
because  (1)  the  experimental  flutter  tests  were  limited  to 
small  angles  of  attack  and  (2)  there  exist  no  reliable  and 
fully-developed  techniques  for  obtaining  the  vibration 
characteristics  of  tjic  deformed  models.  Nevertheless, 
this  preliminary  work  is  presented  below  for  possible 
future  use;  calculations  for  some  simplified,  two-dimensional 
cases  are  also  included,  to  obtain  a  qualitative  assessment 
o'f  the  effects  of  bluntness  and  initial  angle  of  attack. 

2.  Aerodynamic  Forces  on  a  Blunt  Deformable 

Surface  at  Large  Angles  of  Attack,  According 
to  Newtonian  Theory 

For  lifting  surfaces  at  high  angles  of  attack  and  at 
high  speeds  (J^ \>y  i.  ),  the  simplest  way  of  determining 
steady  or  unsteady  aerodynamic  loads  is  that  following  the 
raw  Newtonian  theory  or  the  impact  theory  (Ref.  1.6)* 

Without  added  coup  li  cat  ion,  one  may  introduce  Lee's  Modi¬ 
fication  in  the  above,  wherein  all  pressures  are  Modified 
by  a  constant  factor.  The  justification  for  the  latter  is 
based  on  experimental  results.  This  factor  adjusts  the 

*More  precisely,  if  ftsioO  — d.  ,  where  Q  is  the 
local  surface  inclination.  Thus  piston  theory  is  invalid 
over  certain  partr  of  the  blunt  leading  edge. 
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stagnation  point  pressure  coefficient  to  coincide  with 

the  observed  C p,  .  The  pressure  coefficient  according  to 
Newtonian  theory  is 

Cp  *  Cj>  [(*-  -y")'71]  (I*37) 

f,  is  the  velocity  of  the  body  point  where  the  pressure 
is  to  be  determined,  and  71  is  the  unit  normal  vector, 
directed  in.  The  pressure  coefficient  at  the  stagnation 
point  is 

CK  "  ~  «•“> 

The  observed  varies  between.  1.8  and  2.0  depending 

on  the  configuration.  Equation  (1.37)  is,  in  a  rigorous 
treatment,  applicable  only  to  those  regions  of  the  body 
where  the  instantaneous  local  inclination  0  of  the  surface 
relative  to  the  stream  is  large,  in  the  sense  that 

Msind»l  (1.39) 

However,  for  lack  of  a  more  rational  approach,  Eq.  (1.37) 
may  be  utilized,  provided  that  Eq.  (1.39)  is  satisfied 
over  most  of  the  body.  For  the  so-called  shadow  region, 
where  (9  <0  »  is  assumed  to  be  zero.  At  high  Mach 

numbers,  it  is  consistent  with  the  above  approximations  to 
sec 

?“/:[(*-$-)•*]  a.*o> 

where  fo  SL  ,  since  the  dynamic  pressure  will 

be  much  larger  than  £  . 
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The  main  task  la  the  detatmlnatian  of  J  end  . 

The  complexity  of  this  proceed  will  depend  an  the  particular 
problem  at  hand.  Xt  may  be  possible  to  treat  the  steady 

and  unsteady  phaaea  simultaneously,  or  it  may  here  to  be  , 

carried  out  in  the  following  two  phaaea: 

(1)  Calculating  the  sharp*  under  a  static  loading, 
corresponding  to  strady  flight,  and 

(2)  Proceeding  to  the  unsteady  part,  using  the 
steady  deformed  fc!--'te  as  the  starting  point, 
assuming  small  (in^uiteaimel)  time-dependent 
deformations. 

Phase  (1)  may  or  may  not  be  affected  by  linear  aeroelastie 

analysis,  for  even  though  the  structure  may  behave  linearly, 

the  airloads  associated  with  the  static  deflections  any  be  * 

large.  To  convey  these  thoughts,  consider  the  simple  model 

of  a  cantilever  rectangular  wing,  in  its  uadsformad  state. 

(See  sketch  below).  X,  is  coincident  with  the  "center  * 

line"  of  the  crose  section,  is  normal  to  the1  midp  lane, 
i.e.,  (x,  —  plane.  Under  load,  the  wing  is  assumed  to 
behave  as  a  beam-rod.  Typically,  the  deformed  atete  at  any 
sparvise  station  Xt  may  be  cbought  of  aa  that  due  to  the 
sum  of  (a)  a  translation  4m((xAl*)  perpendicular  to  the  mid- 
plane,  and  (b)  e  rotation  a c»  about  am  axis 

%t  »  x*.  •  The  shape  of  the  cross  section  is  sssi— d  to 
remain  unaltered  under  load,  or  alternately  stated,  the 
deformations  are  deecribeble  entirely  by  the  deflection  of 
the  center  line.  Let  X,  -  /  (x, ,  )  denote  the  agnation 

of  the  surface  in  its  unloaded  position.  After  deformation, 
assume  each  section  ct  X*  remains  at  Xg  •  A*  equation 
of  the  surface  in  terms  of  body  fixed  axes  t  xj  la: 

*1 <x‘4l> 
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Three  secs  of  axes  are  now  considered. 

(1)  X ,  <4  , 2  ,  with  X  aligned  in  the  free- 

0 

stream  direction,  and  having  unit  vectors 

t.f.t 

(2)  ,  Xz  ,  Xj  with  unit  vectors  Tft  and 

(3)  x/»  >  TCj  ,  the  body-fixed  axes  with  . 

t;»  r*' .  • 

The  following  relations  hold  between  the  coordinate 
systems: 

X/  *  fx,- XA) 


and 


=  **)  s,n  *  +■  (  xi  +  co3  ** 

4*/  =  "teas  (oC'i-oi)-  Aai'n  (oe*  +  oc) 


(1.42) 


*  ~J 

~tj  =  i'sm  £(,-f  0t)  -t  A  Cos  (ote+cc) 


(1.43) 


The  above  configuration  may  be  thought  of  as  a  "sheared” 
model,  in  that  during  deformation  each  section  stays  in 
its  Xj-  p lane  and  keeps  its  shape.  Then  the  vector  ^ 
which  represents  the  absolute  velocity  (relative  to  x,  y, 
z- system)  of  a  point  on  the  body,  does  not  possess  a  com¬ 
ponent  in  the  J -direction.  Since  any  body  point  .has  zero 
velocity  relative  to  the  system  f  }  “X^  )# it 

follows  that 


At 


dW 

At 


(1.44) 
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t 


t 


Differentiated  forms  of  Eqs.  (1.42)  Chen  yield  in  conjunction 
with  Eqs.  (1.44) 

t  -  4s/»«]  +  (1.45) 

Here  the  relations 

T.  =  ~t!  cascx-hT/sw* 

J  (1.46) 

L 3  a  -  /.JstfK*  +  <3  cosoc 


have  also  been*  utilized.  The  unit  normal  vector  tt,  may  be 
obtained  from 


7L 


ivsi 


(1.47) 


where 


&  .*•  B (x{  ,xtJx^)=  xj-f(x,' jXj!)=*o 


(1.48) 


represents  the  equation  of  the  body  shape.  Also  note  that 
in  Eq.  (1.45)  one  may  set  X^m/  since  refers  to 

the  velocity  of  a  body  point.  After  some  algebra,  there 
results 


-f,e*fa+ec)+SHl  (<X.+* ) 

-/i+VT-t-F*1 


(1.49) 


V/+  (Tf+F*  ' 


(1.50) 


4 
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(1.51) 


where 

rr  _  df 

Jx/  ' 

When  Eqs.  (1.49)  and  (1.50)  are  substituted  in  £q.  (1.40), 
one  obtains 


/-/A'*')  -  /'•  if 


L+y-'f+F* 


52) 


for  those  regions  of  the  body  with  positive  inclination  9  . 

The  latter  expression  may  be  linearized  with  respect 
to: 

(1)  steady  and  unsteady  deformations  simultane¬ 
ously, 

(2)  unsteady  deformations  alone. 

Let  oC  and  A.  be  split  up  into  two  parts,  one  steady 
(denoted  by  subscript  s )  and  one  time-dependent  (denoted 
by  subscript  «.). 


Case  1:  e^Sfeetk/^s,Au  are  all  small. 

COS  (&%+x)  =  CCS  oc#-  (or,  ret^smeu 

$in(oCt+(X)Z T  (cXg-nxJcast  Smc g 

•  *  ;  ; 

&.**<*+  f  *■  m  ,  3">OC~  Xfi  CffS&oti 
Utilizing  these  relations,  one  arrives  at 
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b 


A 


jfc"?#,-/ case^  f  Z(stMa-  fctsat.). 


[/'&**+)*"><*,  t(*s  +*J**y.  +fc]J  (I*  53) 


Case  2:  o£u,Ju. 


are  small,  so  that 
cosfa +<x)  =  iiofaHCs); 

et-  fX^,  A  -  ;  s'n  <x-sir>o(jctia>scx3  ■  cos** caso^-#* z*cf3 

The  corresponding  pressure  is 

k  _  A  _ . -  I &"/'*/«,)- +<*>)- b*>k+«*)l* 

r  H-Vf+r*  [L  u 


•j«m  (f»*K+«s)+ 


54) 


Appropriate  integrations  of  Eq.  (1.53)  or  Eq.  (1.54)  lead 
to  the  total  aerodynamic  forces.  The  choice  between 
Eq.  (1.53)  and  Eq.  (1.54)  will  depend  on  the  particular 
problem  at  hand.  If  CC3  is  not  small,  i.e. ,  large  initial 
steady  deformations  are  present,  method  (2)  may  have  to  be 
used. 

The  main  disadvantage  of  the  above  derivation  Is 
that  it  assumes  simple  deformations  of  the  beam-rod  type. 

The  much  more  complicated  case  arises  when  deflec¬ 
tions  are  not  easily  describable,  as  in  the  plate-type, 
inflatable  delta  structures.  To  make  computations  tractable, 
it  is  necessary  to  resort  to  other  simplifications.  One  such 
case  is  when  one  assumes  that  surface  points  have  only  de¬ 
formations  normal  to  the  midplane  of  the  planform,  i.e., 
a  point  at  (x,  ,  xt ,  X3)  goes  to  the  deflected  position 
(under  load)  at  ( X, ,  Xii  x^-J)  »  where  &  is  the  deflection, 
positive  down.  This  approach  would  be  inaccurate  in  the 
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blunt  regions.  Following  the  previous  notation,  one  may 
write 


i  -  -7>* 

3  ~  -  /A,.  *2.)  -t  S (xt,  xit±) 


V/  f  (S'-P'fi-fr-i? 


where 

f'=  i£ 
1  ax, 


rr  dS 


After  some  manipulation,  one  finds  that 


[cosou(iLf')  +  Sm0(o  +  ~r] 

[l  +  (?-{'?+ 


(1.55) 


Since  the  mid-surface  of  the  structure  is  required  to  be 
plane,  &  must  be  taken  to  be  the  sum  of  the  steady  and 
time-dependent  deflections.  One  added  source  of  inaccuracy 
arises  here.  Newtonian  theory  requires  high  A jsinQ  , 

which  means  that  the  initial  angles  of  attack  oust  be  large. 
For  the  present  models,  large  angles  of  attack  result  in 
excessive  model  deformations,  so  that  one  can  no-  longer 
assume  this  type  of  normal  deformation.  The  preceding 
Eq.  (1.55)  is  no  longer  adequate  and  other  methods  of 
analyses  oust  be  devised.  One  possible  approach  would  be 
to  first  consider  the  deformation  of  the  wing  under  steady 
load,  and  then  treat  the  unsteady  motion  as  a  perturbation 
about  the  steady  deformed  position. 


i 

f 

i 


i 
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3.  Applications  to  Two-Dimensional  Sections 

So ne  calculations  have  been  carried  out  on  the 
typical  section  (corresponding  to  the  spanwise  sections 
of  the  bean-rod-type  structure)  to  obtain  qualitative 
assessment  of  the  effects  of  bluntness  and  initial  angle  of ( 
attack.  Specifically,  the  cases  treated  are  those  given  in' 
the  sketch  below. 


(a)  Circular  Nose  (b)  Pointed  Nose 

As  shown  in  the  sketch,  the  sections  are  situated  in  tueir 
deformed  state  under  steady  aerodynamic  loads  only.  Ex¬ 
pressions  were  obtained  for  the  unsteady  airloads  as¬ 
sociated  with  a  translation  (perpendicular  to  %f- 

axis) ,  and  a  rotation  Sr^e1**  (about  x »  X*.  ) .  The  aero¬ 
dynamic  loads  of  interest  are: 

L  Force'  perpendicular  to  x, ,  positive  in  the 

-  Xj -direction, 

Moment  about  any  axis  X,«  X*,  positive  nose  vip. 
It  is  noted  that  in  the  derivations  of  L  and  M-h(  for  %„/•  X«), 
the  axial  force  A  ,  in  the  direction  of  x,  is  also  re¬ 
quired,  although  it  does  not  explicitly  appear  in  the  final 
equations.  These  expressions  are  not  reproduced  here,  for 
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they  -ire  rather  lengthy  and  rlo  not  bear  directly  on  the 
discussions  to  follow. 

Some  observations  can  be  made  as  to  the  behavior  of 
the  "center  of  pressure"  movement  due  to  initial  angle  of 
attack  and  leading  and  trailing  edge  bluntness.  In  the 
limit  as  the  oscillation  frequency  to—yo  ,  the  -motion 
effects  vanish  and  the  center  of  rotation  X.  becomes  im- 
material. 

The  point  JL ^  Xm  about  which  the  moment  M  vanishes 
is  denoted  by  xc  .  the  center  of  pressure.  For  slow  oscil¬ 
lations,  the  ccnter-of -pressure  location  relative  to  may 
be  important  in  flutter  calculations.  gives  a 

measure  of  the  real  part  of  the  complex  Af*.  due  to  the 

rotation  £  - 

u 

The  variation  in  center-of -pressure  location, 
as  calculated  for  the  circular  nose  and  pointed  nose  sec¬ 
tions  at  various  angles  of  attack,is  presented  in  Fig. I. 25a. 
A  measure  of  bluntness  is  the  ratio  £-»/?  >  If  the  chord 
is  taken  to  be  unity) .  For  R  -yo  *  -he  section  becomes 
a  flat  plate  and  -  0.c^t  regardless  of  oC0  ,  as  one 

would  expect.  R  — y  1/2  corresponds  to  circular  and  square 
cylinders  for  the  circular  and  pointed  noses,  respectively. 
For  the  case  of  oC0m  45  deg,  X,.  remains  at  0.50  regardless 
of  the  bluntnoss  R  .  For  lower  &0's,  Xe  is  forward  of 
the  50-per-ccnt .chord,  while  the  reverse  is  true  for 
larger  of4's.  Also  the  Xc -movements  are  less  for  the 

circular  than  for  the  pointed  nose. 

The  above  calculations  are  essentially  those  of  a 
small  steady  angle  of  attack  superimposed  on  the  median 
angle  (  oC0  +<Xj  )•  For  this  case,  another  interesting 
comparison  can  be  made  between  the  force  L  for  vari>  ua 
R's  and  (oC0  +  et^ys.  The  plots  given  in  Fig.  1. 25b 
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correspond  to  the  expressions: 

Jl.  at  sin  Zfa+CXs)  for  flat  plate  (1-56) 

,  k0^ 

3=  (l-2fC)sin2fat+c(i)  for  circular  nose  (1-57) 

st  ^ - 3 R^s /V> 2 (tXg+cxJ  for  pointed  nose  with  (oc#fo%)?  45 

s  ( !-&£)$;» 2fa,+ oO  ] 

+  £(»-2cas2fo+0^)  for  pointed  nose  with  («.+<%) 4 4&Jy 

The  delta  wings,  such  as  those  considered  in  the 
experimental  program,  have  spanwise  cross  sections  which 
f  are  similar  to  the  circular  nose  typical  section.  -S- 

hewever  is  variable,  increasing  as  one  takes  sections 
further  away  from  the  root.  The  above  results  (Eqs.  (1.56)  - 
‘  (1.57))  suggest  that  perhaps  one  can  account  approximately 

for  the  leading-  and  trai ling-edge  bluntnesses  by  discarding  f 
from  aerodynamic  consideration  the  rounded  portions  of  the 
leading  and  trailing  edges.  A  similar  conclusion  was 
reached  in  Section  B  of  this  appendix,  based  on  physical 
considerations . 
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Figure  Z.9.  Theoretical  Flutter  Parameter  (V^/h^oj)  and 
Static  Divergence  Pa  reset  er  ve 

uM  for  Stated  Condition*. 
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1.2 


Degrees  of  Freedom:  1—  8  2—  Modes 
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Figure  1. 10.  Theoretical  Flutter  Frequency  Parameter 
(ac  j/ojj)  va  uM  for  Stated  Conditlona. 


173 


'T>*r 


FDL  TDR  64-147 


Figure  1.11.  Theoretical  Flutter  Parameter  and 

Static  Divergence  Parameter  (V^/ bQn)j)  vi 
hM  for  Stated  Conditiona. 
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Oegrees  of  Freedom:  Ilia  2~  Modes 
Mode  Shapes:  B 

Aerodynamics:  Piston  Theory,  Flat  Plate 


(^•)2=0220.9 

g  =  0,  0.05 


V, 


Correspond  lo  Lower  Branches  of  -—i— .  Ve.  uM 

A  bo"2 


.L_ 


g=0 

‘g=0.i 


■  r 

Correspond  to  Upper  Branches  of  ^  vs.  pH 


g=0 


_ I _ 

500 

Figure  1.12. 
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va  pH  for  Stated  Conditions. 
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Figure  1.14.  Theoretical  Flutter  Frequency  Parameter 
(dBf/'oj)  ^  uH  for  Stated  Conditions. 
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Oegrees  of  Freedom:  l“8  ^Modti 
Modes  Shapes:  C-l 

Aerodynamics:  Piston  Tbeory.Flct  Plot# 
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Figure  1.15.  Theoretical  Flutter  Parameter  (V^/b^) 
and  Flutter  Frequency  Parameter  (m^ «j) 
vs  UM  for  Stated  Condition#. 
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Figure  1.16.  Theoretical  Flutter  Parameter 

and  Flutter  Frequency  Parameter  (od^/cd^) 
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Degrees  of  Freedom:  2.—  a  3“  Modes 
Mode  Shapes:  C-i 

Aerodynamics:  Piston  Theory,  Flat  Plate 
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Figure  X.17.  Theoretical  Flutter  Parameter  (V^c/t^ouj) 
and  Flutter  Frequency  Parameter  ( 
vs  uM  for  Stated  Conditions. 
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Degrees  of  Freedom:  2— 0k  3—  Modes 
Mode  Shapes:  C-1  ,C-2 
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Figure  1.18.  Theoretical  Flutter  Parameter  (V^/f^avg) 
and  Flutter  Frequency  Parameter  (cd if/a^) 
vs  pM  for  Stated  Conditions. 


Figure  1.19.  Theoretical  Flutter  Per— *  ter  (Vj/b^) 
and  Flutter  Frequency  Parameter  (og/aj) 


ve  uM  for  Stated  Conditlooa. 
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pM. 

Figure  1.20.  Theoretic*!  Flutter  Parameter  (Vj/b^) 
end  Flutter  Frequency  Parameter  (cof/cog) 
vi  for  Stated  Condltloni. 
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Figure  1.22.  Theoretical  Flutter  Parameter  (Vj/b^) 
and  Flutter  Frequency  Parameter 
vs  uM  for  Stated  Conditions. 
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Figure  1.23.  Theoretical  Flutter  Paraaeter  (Vj/b^a^) 
and  Flutter  Frequency  Paraaeter  (Of/ap 
va  uM  for  Stated  Conditions. 
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25.  Variations  of  Center>of •Pressure  Location  and  Force  L  with  Bluntness 
and  Initial  Angle  of  Attack,  for  Various  TVro-Dimensional.  Sections. 
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APPENDIX  IX 

MEMBRANE-TYPE  FLUTTER  OF  AN  INFLATABLE  PANEL 
by  Marten  T.  Landahl 

A.  SIMPLIFIED  MATHEMATICAL  MODEL 

Following  Miles  (Ref.  II.  1)  we  will  consider  the 
stability  of  two-dimensional  waves  travelling  on  a 
stationary  panel  of  infinitive  length  with  a  deformation 
of  the  fora 

n*  s*  Pee.  j 'hS'e  >  (ii. l) 

where  is  the  instantaneous  deformation  (assumed  in¬ 
finitesimal)  ,  -4  »  (ztt/x)  is  the  wave  number,  C  is  the 
wave  velocity,  and  X  the  wave  length.  It  will  be 
assumed  that  the  panel  is  exposed  to  flow  of  a  constant 
Mach  number  M—  (y/a.)  and  density  P  on  both  sides 
(Fig.  II. I).' 


This  model  will  be  useful  for  an  approximate  analy 
sis  of  waves  of  short  wave  length  so  that  edge  effects 
may  be  neglected.  Also,  it  can  be  used  as  an  approxima¬ 
tion  for  nonunifora  outside  flow  (as  would  be  created  by 
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Che  disturbances  due  Co  Che  vehicle  itself),  provided  the 
wave-length  X  Is  so  small  thaC  flow  parameters  vary 
liccle  over  Che  distance  of  one  wave  length.  A  more 
general  analysis  would  require  different  flows  on  Che  two 
sides  of  Che  panel,  buC  this  would  unduly  complicaCe  the 
analysis.  For  Che  case  of  vacuum,  or  very  low  velocity 
on  one  side  the  present  analysis  will  hold  If  the  aero¬ 
dynamic  forcing  term  in  the  equations  below  is  reduced  by 
one  half  its  value. 

B.  AERODYNAMIC  CONSIDERATIONS 

The  disturbance  pressure  created  by  a  small  ampli¬ 
tude  travelling  wave  of  the  form  (II.  1)  may  easily  be 
calculated,  since  it  must  be  the  same  as  that  on  a  steady 
wavy  wall  in  a  stream  of  velocity  (V-fi).  One  thus  easily 
finds  that  the  pressure  induced  on  the  upper  surface  of  the 
panel  ia 


\ll2y 

Vir«Y'-fT 

(The  denominator  is  simply  the  Prandtl-Glauert  correction). 
Similarly,  for  a  wave  that  is  supersonic  with  respect  to 
the  stream  [(V-c)>a.  J  one  obtains 


(IX. 3) 


V 


« 
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Since  there  is  an  equal  pressure  disturbance,  but  witb 
opposite  sign,  on  the  lever  surface,  the  total  aero¬ 
dynamic  loading  on  the  panel  is 


Lk(x-  ct)  Z$k  (Pcf 


w-e. 


ik(x-ct) 


(11.4) 


for  a  subsonic  wave,  and 


£  St 


«I.5> 


for  a  supersonic  wave. 

These  loadings  are  conveniently  expressed  in  fora 
of  an  “aerodynamic  impedance" 


'  ztt(V-c.f 
.  -h-rfit-Sjf- 

2.  fit  (V-cf 


for  subsonic  wave  (II. 6) 


for  supersonic  wave 
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C.  STRUCTURAL  CONSIDERATIONS 

The  two-dimensional  inflatable  panel  possesses 
shear  deformation  and  rotary  inertia,  and  as  such  can  be 
represented  as  a  Timoshenko  beam  (see  for  example.  Ref. 

II. 2,  Chap.  3,  Eqs.  3-4  to  3-6).  The  governing  panel 
equation  would  then  be 

-  (r<-  jrY'hx  *  w *  *  fr * 

4  J  CII.7) 

'  t  + .%[ t  ~  jprS** 

where  D  *»  ES it  /z  is  the  overall  panel  bending 
stiffness,  X  =  is  the  rotary  inertia,  **»  Is 

the  mass  per  unit  length,  p^K  is  the  shear  stiffness 
due  to  internal  pressure  and  £  is  loading  per  unit  length. 
The  contribution  of  drop  chord  mass  to  X  and  nt  may  be 
neglected. 

To  give  the  panel  the  travelling  wave  deflection 
1**-  of  Eq.  (II.  1)  would  require  a  loading  ^  of  the 
type 


ik  (t-cty 


<H.8> 


Upon  introduction  of  this  and  Eq.  (II.  1)  into  Eq.  (II. 7) , 
one  obtains 
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As  with  the  aerodynamic  forcing  term  we  write  this  result 
in  fora  of  a  ’’structural  impedance1'  ZiS  »  given  by 


luJL 

S 


t* 

:  m. 


(f  <■&-*) 


It  is  instructive  to  express  this  as  follows: 

z  _ 

s  (c.x-  S) 


(XX. 11) 


where 


c‘  mz{-ir*  r*wN*tT  x  r" 


/•*  <fhti .  a  .  bi\.  (ii. u> 

4  ’£(-sr*r*27V*'ir  J  x**'  r* 


3>  ■  6*  (11.14) 

X  *  X4* 

The  ing>edance  of  the  structure  is  thus  zero  when  excited 
at  the  wave  velocities  C&£Ct*  and  C  *  •  Hence, 

these  are  the  phase  velocities  for  free  surface  waves 
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(in  vacuum).  At  C=»Ca  the  structure  has  infinite  rigid¬ 
ity,  i.e.,  it  acts  as  a  "dynamic  vibration  damper.”  It 
turns  out  that 


C,  <  4  <  q 


(11.15) 


so  that,  in  principle,  the  structural  impedance  varies 
with  wave  velocity  as  shewn  in  Fig.  11.2. 


Figure  II. 2 

% 

To  simplify  the  discussion  we  will  consider  the  following 
two  extreme  cases: 

a)  Pure  Shear 

Pure  shear  deformation  is  obtained  for  D-m*  with  }>k 
fin'.te.  Then,  (11.12)  -  (11.14)  give 
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t 

r*~+>  A£ 

i  tn 

(11.16) 

» 

x  a  J> 

^  f/0  ^  «■—'#» 

(II. 16a) 

Hence 

(11.17) 

% 


% 


4 


b)  Pure  Bending 

In  this  case  we  let  keeping  D  finite.  This 

gives 


'  r^v-Tn. 

(11.18) 

(11.19) 

r  *■  ,  hi. 

• 

(11.20) 

which  in  turn  leads  to 

• 

% 

(■A  w 

(11.21) 

J 

Practical  configurations  tend  to  lie  between  these  two 
extremes,  perhaps  closer  to  the  shear  case. 

In  both  cases,  the  structural  ispedance  is  of  the 


I 

I 
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fora 

where,  for  pure  shear. 


y/l  *  772 

C» 


(IX. 22) 


(11. 23) 

(11.24) 


whereas,  for  pure  bending, 

=.  m  +  Z&* 

c/.  *'» 

7n+  Li 


(11.25) 

(11.26) 


The  expression  (11.22)  will  be  used  in  Che  following  dis¬ 
cussion.  By  inserting  typical  parameter  values  it  turns 
out  that  £e  is  of  the  order  100  -  200  ft/sec.,  i.e., 
rather  snail  coaparcd  to  flight  velocities  of  interest. 
This  is  of  crucial  importance  for  the  instability  problem 
as  will  he  seen  below. 


D.  STRUCTURAL  DAMPING 

To  avoid  complications,  structural  damping  will  be 
introduced,  without  consideration  of  how  it  originates 
physically,  slnply  by  replacing  by 

Ce  (f-  i£j  (11.27) 
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i 


f 


where  £,&  is  an  equivalent  ’’damping  ratio”  such  that 
the  free  surface  waves  are  reduced  in  amplitude  by  the 
factor 


-  ZTtS- 

& 

(11.28) 

over  each  cycle  of  an  oscillation.  In  principle,  £& 
could  be  obtained  from  experiments.  It  will  be  assumed 
that  the  damping  is  small  so  that  terms  that  are  quadratic 
in  £  may  be  neglected. 

With  danping,  the  impedance  (11.22)  that  before 
was  purely  real  now  becomes  complex,  so  that 

3,  - 

+  *  % 'si 

(11.29) 

where 

4- 

(11.30) 

(11.31) 

B.  EIGENVALUE  PROBLEM 

Self-sustained  motion  can  occur  when  the  aero¬ 
dynamic  and  structural  impedances  match,-  i.e.,  when 

ZA  (c)  =  Zs  (c.)  (11.32) 

where  ZA  is  given  by  (IX.  6),  and  (11.22)  will  be  used 
for  Thg  •  This  equation  can  then  be  solved  for  the 
(generally  coop  lex)  eigenvalue  (a)  C  of  the  problem  for 
a  given  set  of  system  parameters  and  A  .  If  the  imagi¬ 
nary  part  of  the  eigenvalue  C  is  positive,  the  waves  will 
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grow  exponentially  with  time,  i.e.,  flutter  sets  in.  In 
view  of  the  fact  that  the  aerodynamic  impedance  changes 
from  purely  real  (for  real  c  )  for  a  subsonic  wave  to 
purely  imaginary  for  a  supersonic  wave,  the  two  cases  are 
best  discussed  separately. 


F.  STABILITY  OF  SUBSONIC  WAVES 


The  simplest  approach  for  a  wave  that  is  subsonic 
relative  to  the  stream  is  to  first  consider  a  panel  with 
zero  structural  damping  for  which  the  eigenvalue  relation 
(11.32)  takes  the  form 


zfA(V-o)‘- 


(11.33) 


The  arch  type  of  a  subsonic  wave  is  the  incompressible 
one,  M**0.  The  incompressible  cases  have  been  thoroughly 
discussed  by  Landahl  (Ref.  II. 3)  and  by  Brooke  Benjamin 
(Ref.  II. 4). 

For  real C  *  (11.33)  is  purely  real  and  expresses  a 
balance  between  the  aerodynamic  stiffness  (on  the  left* 
hand  side  of  the  equation)  and  the  mechanical  stiffness  • 
inertial  reaction  (right-hand  side).  It  was  shown  in 
Ref.  II. 3  that,  for  MmO ,  the  two  roots  of  (11.33)  are 
real  for  sufficiently  low  V*  and  thus  correspond  to  two 
neutrally  stable  waves  labelled  Class  A  and  Class  B. 
respectively.  The  Class  A  wave  speed  is  essentially  the 
root  O for  fmO  *  modified  by  the  presence  of  the  air 
flow  and  thus  increased.  The  Class  B  wave  speed  corresponds 
to  the  free-surface  wave  speed  C  “  ,  modified  (lowered) 

by  the  flew.  For  Increasing  forward  speeds  or  densities, 
there  may  eventually  be  no  wave  speed  for  which  the  two 
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sides  of  (11.33)  balance,  and  the  resultant  mechanical 
stiffness  will  consequently  always  be  insufficient  to  bal¬ 
ance  the  aerodynamic  forces.  In  such  a  case,  the  two  eigen- 
values  for  C  become  complex  conjugates  indicating  the  on¬ 
set  of  a  violent  Class  C  or  Kelvin-Eelmholtz-type  insta¬ 
bility.  This  was  the  type  of  flutter  considered  by  Miles 
(Ref.  II. 1). 

The  effect  of  (small)  structural  damping  on  sta¬ 
bility  was  considered  in  detail  by  Landaul  (Ref.  II. 3). 

He  demonstrated  that  dasping  may  make  the  Class  A  wave 
(mildly)  unstable,  but  stabilizes  the  Class  B  wave.  In 
contrast  it  has  very  little  effect  an  the  Class  C  wave.  It 
was  also  shown  by  Brooke  Benjamin  (Ref.  II. 4)  that  for  high 
,  density  ratios  ,  the  speed  margin  between  the  on¬ 

set  of  Class  A  and  Class  C  instabilities  is  inversely  pro¬ 
portional  to  yt l.  It  would  therefore  seem  that  Class  C  in¬ 
stability  is  the  type  of  primary  technical  interest. 

It  also  follows  from  (11.33)  that  the  wave  velocity 
for  either  Class  A  or  Class  C  is  always  less  than  C,  . 

Since  for  cases  of  practical  interest  for  the  present 
problem  is  likely  to  be  small  (  <  200  ft/sec,  say)  a  sub¬ 
sonic-type  wave  cannot  be  maintained  for  free  stream  speeds 
ouch  beyond  .sonic,  say  beyond  .  From  a  practical 

standpoint,  the  stability  cf  supersonic  waves  is  therefore 
more  of  interest. 


G.  STABILITY  OF  SUPERSONIC  HAVES  . 

For  a  wave  that  is  supersonic  with  respect  to  the 
free  stream  the  eigenvalue  relation  (11.32)  becomes 


-zlfifV-cf  _ 
vV/jt-AP-i'  " 


(11.34) 
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> 


Here  we  have  included  structural  damping  because  without  it 
no  neutrally  stable  wave  can  exist.  By  separating  the  real 
and  imaginary  parts  of  (II. 34)  we  find  that  flutter  sets  in 
for 


(11.35) 

(11.36) 


Here  we  have  neglected  terms  quadratic  in  .  For  Mach 
numbers  well  into  the  supersonic  region  (transonic  Mach 
numbers  will  be  considered  later)  Q  may  be  neglected 
compared  to  and  we  thus  obtain  from  (11.36) 

v-  <u-37> 

where 

A' 

It  follows  from  (11.37)  that  since  the  flutter  velocity  is 
proportional  to  the  yfjfc  ,  the  long  waves  are  more  unstable 
than  short  ones.  However,  the  analysis  was  based  oh  the 
assunftion  that  the  wave  length  is  much  smaller  than  the 
chord.  This  is  due  to  the  fact  that  the  structural  damping 
decreases  faster  with  decreasing  wave  number  (as  k  )  than 
does  the  aerodynamic  out-of-phase  forcing  pressure.  It  Is 
therefore  unlikely  that  membrane  type  flutter  will  occur 
before  ordinary-type  flutter. 
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V 


I 


» 


H.  EFFECT  OF  BOUNDARY  LAYER 

Equation  (11.37)  indicates  that  for  M~**L  the 
flutter  velocity  goes  to  zero.  It  wa3  shown  by  McClure 
(Ref.  II. 5)  that  this  conservative  result  is  due  to  the 
neglect  of  the  boundary  layer.  McClure's  calculations  were 
based  on  shear-flow  stability  concepts  and  are  too  compli¬ 
cated  to  be  reproduced  here.  However,  the  basic  features 
may  be  studied  by  using  a  very  simple  flow  model  proposed 
by  Fung  (Ref.  II. 6).  In  this,  the  boundary  layer  is  repre¬ 
sented  by  a  layer  of  constant  velocity  less  than  the  free- 
stream  velocity  (see  Figure  11.3). 


JL 

Mi,  Pd ,  a*  <* 

1 

i 

Figure  II. 3 


Assuming  that  the  wave  length  is  much  longer  than  the 
thickness  of  the  boundary  layer,  wie  may  treat  the  flow  in 
the  boundary  layer  as  one-dimensional.  One  then  finds  by 
matching  pressures  and  deformations  across  the  slip  line 
thee 


-z:4f(v-cf 

J 


(11.38) 


This  result  reduces  to  the  previous  one,  e.g.,  (11.6),  for 

eh 
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Again,  it  has  been  assumed  that  there  is  an  identi¬ 
cal  flow  over  the  lower  panel  surface.  Since  the  impedance 
(11.38)  now  has  both  a  real  and  imaginary  part,  the  eigen¬ 
value  problem  becomes  much  more  complicated.  However,  con¬ 
siderable  insight  might  be  gained  by  a  simplified  treatment 
along  the  following  lines': 

Since  the  negative  of  the  imaginary  part  of  ZtA 
represents  the  force  component  in  phase  with  the  transverse 
plate  velocity,  it  gives  a  measure  of  the  energy  input  per 
unit  area  which,  for  a  neutrally  stable  wave,  oust  be  bal¬ 
anced  by  the  structural  energy  dissipation  in  the  panel 
which  is  proportional,  to  £  .  Whereas  according  to  (II. 6) 

this  energy  input  is  infinite  when  the  wave  travels  at 
sonic  speed,  (11.38)  shows  that,  in  fact,  for  this  speed 
the  boundary  layer  will  make  it  zero.  By  separating  the 
real  and  imaginary  parts  of  ZA  ,  one  finds  that,  for  given 
?,  V,  C,  ,  Vj  ,  &  and  .the  imaginary  part,  £AT 

is  minimum  for  a  Mach  number  given  by 


■&d 


tOHtlr 


(11.39) 


The  minimum  value  is  then  found  to  be 


(11.40) 


with  the  associated  real  part  being  equal  but  with  opposlta 
sign,  thus 
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i 


r 


\ 


fc  ' 


V/,.  - 

at  m*t. 


k  &  -A 


The  boundary- layer  model  would  require  that 


V<  =  *V 


(11.41) 


(11.42) 


where  is  same  constant  less  than  unity  but  probably 
greater  than  0.5.  Since  the  wave  velocity  is  mainly  con¬ 
trolled  by  Cg  ,  the  propagation  velocity  of  free  surface 
waves  in  vacuum,  we  may  therefore,  using  the  same  arguments 
as  before,  assume  that  C  may  be  neglected  compared  to  V"  , 
and  hence  also  compared  to  "]£  .  Then  (11.40)  and  (11.41) 
simplify  to 


(11.43) 


4 


where 


Now  substituting  these  back  into  the  eigenvalue  relation 
(11.32)  we  find  that,  at  the  maxiimim  energy  input  condition. 


(11.44) 
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and 


l£.  ^  =  m  I'fcJ'-c)  (11.45) 

The  second  of  these,  (11.45),  merely  shews  that,  for  high 
mass  ratios,  the  wave  velocity  will  be  slightly  below  Cg  . 

The  first  relation  is  the  interesting  one  and  represents  an 
energy  balance  at  the  condition  of  maximum  aerodynamic 
energy  input.  Solving  it  for  the  flutter  velocity  we 
obtain 

•  x  ~  ('-*<)  •  <«•«> 

which,  again,  shows  that  the  low  wave  numbers,  l.e.,  large  • 

wave  lengths  are  the  most  unstable  one.  It  Is  also  in* 

teresting  to  note  that  thickening  of  the  boundary  layer  can 

stabilize  membrane-type  flutter  at  transonic  speeds.  This 

should  also  hold  for  panel  flutter  in  general,  the  physical 

explanation  being  that  the  low-speed  boundary-layer  flow 

acts  as  a  "cushion"  to  lower  the  aerodynamic  forces  exerted 

by  the  free  stream  due  to  panel  deflections. 

1.  CONCLUSIONS 

A  simplified  analysis  of  the  stability  of  travelling 
waves  over  a  flat  ..wo -dimensional  inflatable  panel  indicates 
tha~  short-wave  i.~it ability  is  not  likely  to  occur.  The 
reason  is  that  the  stiffness  and  resistance  of  the  structure 
decreases  faster  (as  )  than  does  the  aerodynamic  forcing 
term  (as  it).  Therefore,  the  longest  waves  would  be  s»st 
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unstable  (i.e.,  have  lowest  flutter  velocity)  which  points 
to  the  ordinary  type  (bending-torsion)  flutter  as  being  the 
most  critical.  The  situation  would  be  different  for  a 
curved  inflatable  panel,  or  membrane,  because  then  the  hoop 
stresses  generated  by  the  travelling  waves  would  act  as  a 
distributed  elastic  support  and  the  structural  stiffness 
would  go  to  a  constant  as  ^  -*0  ,  thus  making  possible 

for  instability  to  occur  at  some  intermediate  wave  number. 

The  situation  would  then  be  similar  to  that  realized  by 
Dugundji,  Dowell,  and  Perkin  (Ref.  II. 7)  In  their  experiments 
with  a  panel  on  continuous  elastic  foundations,  in  which 
travelling-wave- type  flutter  indeed  was  found.  * 

The  very  simple  and  crude  analysis  of  the  effect  of 
a  boundary  layer  shows,  as  previously  demonstrated  by 
McClure  (Ref.  II. S)  that  the  boundary  layer  removes  the 
singularity  in  pressure  predicted  by  linearized  inviscid 
theory  for  a  wave  travelling  at  sonic  velocity  relative  to 
the  air.  Hence,  thickening  of  the  boundary  layer  has  a 
stabilizing  effect  on  waves  travelling  at  transonic  speeds. 


A 


* 
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UST  OF  SYMBOLS  FOR  APPENDIX  II 


t 


\ 


t 


CL 

C 

A 

l 

A 

Ttl 

/> 

h 

t 

u7- 

X 

2) 

B 

M 

V 

z 


Speed  of  sound 
Wave  velocity 

Functions  defined  by  Eqs.  (11.12)  -  (11.14) 
Depth  of  panel,  t.ee  Fig.  11.1 

(-vCT  > 

Wave  number  (  *  27c/ X.  ) 

Mass  per  unit  length 
Pressure 

Differential  pressure  is  shear  stiffness 

due  to  internal  pressure) 

Time 

Time-dependent  lateral  deflection  of  panel 

Streamwise  coordinate 

Panel  bending  stiffness  (  -  EM  / 2-  ) 

Young's  modulus 
Mach  number  (  **  V/<a.  ) 

Free-stream  velocity 
Ixppedance 


§ 

C 

X 

h 

? 

* 


Related  to  M  according  to  {^4^1 ) 
Thickness  of  cover 
Draping  ratio 
Wave  length 

Relative  mass  parameter,  see  Eq.  (IZ.37) 
Free-stream  density 
Constant  in  Gq.  (11.42) 


Subscripts,  etc. 

tL  Quantities  at  base  of  boundary  layer,  except  for 

the  case  of  ,  see  Fig.  II. 3 
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£ 

A 

F 

X 

R  «r  'ffce. 

s 

oo 

D 

< 

( '  > 


"Effective"  values,  see  Eqs.  (11.22)  -  (11.27) 

Aerodynamic 

At  flutter 

Imaginary  part 

Real  part 

Structural 

Free- stream  conditions 
Denotes  amplitude 
Differentiation  with  respect  to  x 
Differentiation  with  respect  to  time 
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APPENDIX  III 

STRUCTURAL  ANALYSIS  OF  AN  INFLATED,  AIEMAT-TYPE  STRUCTURE 

by  James  W.  Mar 


A.  INTRODUCTION 


The  equations  governing  the  behavior  of  a  sandwich- 
type  shell  structure  consisting  of  two  membranes  held  apart 
by  internal  pressure  and  shear  chords  will  be  developed. 
Details  of  the  construction  are  shown  in  Fig.  IIX.l. 


AIR  UNOER 
PRESSURE 


Figure  IIX.l 

The  development  presented  herein  forms  the  basis  for 
the  analysis  of  inflated  structures.  Specialization  and 
.sioplifications  will  be  presented  in  a  subsequent  report. 

B.  THE  GEOMETRY  OF  THE  DEFORMED  MIDDLE  SURFACE 

The  initial  undeformed  middle  surface  of  the  shell 
is  described  by  the  equation 

(iii.d 

where  £  is  the  position  vector  to  points  on  the  middle  sur¬ 
face  and  g* f  are  coordinates  on  the  middle  surface 


211 


FDL  TDR  64-147 


(soe  Fig.  III. 2) 


UNDEFORMED 
MIDDLE  SURFACE 


Figure  III. 2 

Ocher  parameters  of  the  middle  surface  are  listed 


as  follows: 

covariant  base  vectors: 

£  -  — 

* 

(111-2) 

unit  normal: 

.  • 

*  -  a,x?t 

■rzr 

(1X1-3) 

.  ♦ 

discriminant: 

*  - a „**-(*, >f 

(1X1.4) 

fint  fundamental  tensor: 

(XIX .3) 

second  fundamental  tensor: 

l  - 

(III. 6) 
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In  the  above  formulae  and  in  the  subsequent  deriva¬ 
tions  the  indicial  and  summation  conventions  of  tensor 
analysis  will  be  used.  Thus,  Greek  letter  indices  imply 
the  numbers  one  or  tw>  and  a  repeated  index  is  to  be  sutxned 
from  one  to  two.  Also,  the  vector  product  of  two  vectors 
(as  in  Eq.  (111.3))  will  be  symbolized  by  an  "x*  and  the 
scalar  product  of  two  vectors  (as  in  Eq.  (Ill. 5)  will  be 
symbolized  by  a  dat). 

The  deformed  middle  surface  is  located  by  the  position 
vector  as  depicted  in  Fig.  III. 3. 


Figure  I II. 3 


Point  o  on  the  unde  formed  middle  surface  is  displaced  to 
the  point  labeled  0  after  deformation.  This  is  expressed 

as 

*0(1',  iz)  ~  (1x1.7) 
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Other  parameters  of  the  deformed  middle  surface  which 
will  be  of  interest  are  listed  as  follows: 


covariant  base  vectors: 

dg*  * 

(III. 8) 

unit  normal: 

V-  A>  *  A* 

(III. 9) 

discriminant: 

Ast  Att  A22  "  (AJ 

(III. 10) 

first  fundamental  tensor: 

second  fundamental  tensor: 

A.  ->£•>£ 

*  * 

3  -  ft.  id* 

(111. 11) 

(111.12) 

contravariant  base  vectors: 

N  *  A 

(111.13) 

contravariant  metric  tensor: 

’  % 

(III.1A) 

Christoff  el  symbols  of  the  second  \ 

•“  U-r-& 

(Ill. 15) 

Weingarten  formulae: 

if  f 

(IH.  16) 

Gauss  relations: 

(III. 17) 

The  covariant  surface  permutation  tensor,  , 

introduced  in  Eq.  (III.  13),  and  its  contravariant  counter¬ 
part  are  defined  to  have  the  following  values: 
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t 


t„  -£a-o 


(III. 18) 


5,t  --t„-  sr 


(III. 19) 


The  symbol 

values 


is  a  Kronecker  delta  with  the 


(111. 20) 

(111.21) 


C.  THE  GEOMETRY  OF  THE  DEFORMED  SHELL:  THE  STRAIN  TENSOR 

X  -  The  geometries  of  the  uudeforacd  and  deformed  shell 

are  taken  with  reference  to  imaginary  middle  surfaces  which 
are  defined  as  the  locus  of  all  points  midway  between  the 
upper  and  lower  membranes. 
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An  examination  of  Fig.  ZI1.1  discloses  the  following 
geometrical  parameters: 


k 


s  *  OP 


vector  along  shear  chord  from  poinC  o  on 
the  middle  surface  to  a  point  p  on  the 
upper  membrane  of  undeformed  shell 

length  of  shear  chord  in  undeformed  shell, 
i.e.,  spacing  of  membranes  in  undeformed 
shell 

vector  along  shear  chord  in  deformed  shell, 
i.e.,  S  is  the  orientation  of  the  drop 
chord  after  deformation 


2-S  Z]s{  length  of  shear  chord  in  deformed  shell 
(note  that  this  is  not  spacing  of  the  m 
brane  in  the  deformed  shell) 

Point  P  location  of  material  point  f>  after 
deformation 


Point  0 

km 

u.(i‘A2) 


location  of  point  o  after  deformation 
spacing  of  membranes  in  deformed  shell 
displacement  vector  of  middle  surface 


bt 

Gii£, 


base  vectors  on  undeformed  upper  mrntirsne 

base  vectors  on  deformed  upper  instirsne 
associated  with  & 


216 


FDL  TDR  64-147 


ds^  element  of  arc  length  on  undeformed  upper 

membrane 

element  of  arc  length  on  deformed  upper 
membrane 

i 

The  base  vectors  on  the  undeformed  and  deformed 
upper  surface  are  expressed  as 


'd'C  /  dn.  -r  1  L<r-* 

4  —  — —  -h  2b - as  CL  —  —  rt, 

4  *  df1  *  Z  *  r 

/•  _»  3</* 

%  “  &  *  J? 


(in.  22) 


(III. 23) 


tfe  define  the  strain  tensor  for  the  upper  surface  in 
the  usual  manner 


(III. 24) 


Kf,  *  (III'25) 

It  will  be  convenient  to  express  5  in  terms  of  Ax 


and  A/ 


S*  -  +  3*  ft  (III. 26) 

Then  'the  following  relations  can  be  obtained 


f*  -  (s'l~  Ks’)%r+(*''e„+s\)*  im-tn 
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where  che  vertical  bar  before  the  cC  signifies  covariant 
differentiation  with  respect  to  f*  in  the  metric  of  the 
middle  surface  of  the  deformed  shell. 

Relativa  elongations,  or  cxtensional  strain  in  the 
usual  engineering  sense,  are  defined  along  the  two  coordinate 
curves  as 

£  =  >  LzUs^SlL  (no  sun  on  *  )  (III.29) 

*  lt*f  I 

This  leads  to  the  result 


(no  sum  on  oc  )  (111.30) 


In  what  follows,  the  relative  elongations  will  be 
assumed  to  be  small  relative  to  unity.  Thus 


(no  na  on  oc  )  (111.31) 


The  shearing  deformation  is  defined  as  the  change  in 
angle  formed  by  the  base  vectors  ^  and  6*  •  bet 

be  the  angle  between  ^_^and  ^  ,  and  let  b«  tha 

angle  between  and  •  ***** 


(111.32) 


This  leads  to  the  result 


%,2  i-2/iz 


0+e,X'+£*)JMn 


(111.33) 
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f 


l 


* 


ft 


I 


4 


If  the  base  vectors  $  are  orthogonal  (n^vd  generally 
the  are  chosen  such  that  this  is  true) ,  then  -  -|r 
and  »  o  •  The  ^le  <^/2>  will  also  be  assumed  to  be 

snail  such  chat  sin^S.  cf>/z  •  Thus  for  orthogonn  t.  coordinates 
on  the  undefortaed  membrane  the  shearing  angle  is  approxi¬ 

mated  for  the  case  of  small  strains  by 


<111.34) 


The  transverse  shearing  deformations  of  the  shell  is 
defined  by  considering  the  shear  chords  (see  Fig.  HI. 5) . 


The  angle  between  the  base  vectors  and  *5  will 
be  denoted  by  (p^  because  the  original  angles  awe 
90  deg  j since  the  normal  vector  7v  is  perpendicular  to  . 
Proceeding  as  before 

cosfS-6  \ -  d^jt-  (no  suaon  ct  >  (III.35) 

W  izllfj  ■ 


By  definition  (see  Eq.  (III. 33))  the  shearling  com¬ 
ponents  of  the  strain  tensor  are  expressed  as 


(III. 36) 
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Thus  Chere  is  obtained 

tr  i(^)  cm.3» 

Since  Che  shear  rigidity  of  Che  pressurized  sandwich 

is  relatively  low,  the  transverse  shearing  angles  will 

not  be  assumed  to  be  small  at  this  stage  in  the  development. 

• 

It  is  interesting  to  note  that  Eq.  (III. 37)  repre¬ 
sents  the  usual  first  approximation  to  transverse  shear  de¬ 
formations  in  shell  theory.  However,  in  the  case  of  the 
shell  construction  under  consideration,  Eq.  (III. 37)  is  the 
exact  expression  because  the  shear  chords  remain  straight. 

D.  DEFINITION  OF  FORCE  AND  MOMENT  RESULTANTS 

It  will  be  found  convenient  to  define  force  and 
moment  resultant  tensors.  These  are  obtained  by  summing 
the  stresses  across  the  thickness  of  the  shell  to  obtain  a 
set  of  equipollent  forces  and  moments  which  are  assumed  to 
act  at  the  mid-surface  of  the  shell.  An  isolated  element 
of  volume  of  the  deformed  shell  will  be  acted  upon  by  the 
following  kinds  of  forces: 

(a)  Membrane  stresses  in  the  upper  and  lower  membrane 

(b)  The  pressure  of  the  gas 

(c)  The  stresses  in  the  shear  chords 

Each  of  these  will  be  separately  discussed. 

The  isolated  element  of  volume  ta  carved  out  of  the 
shell  by  the  vectors  a^and  as  depicted  in 

Fig.  III. 6.  On  the  upper  and  lower  membranes,  the  edges  of 
the  element  of  volume  are  defined  by  the  vectors  * 

h£ ,  H,  ** Vg 1  ,  and  where  the  enclosed  sub- 
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Figure  III. 6 


scripts  (*.)  and  (/  )  are  tags  identifying  the  ’upper*  and 
’lower1  medbranes, respectively.  These  vectors,  which  are  in 
reality  another  set  of  base  vectors  for  the  upper  and  lower 
membranes ,  are  defined  by  the  equations 

+  X  Jvf  (III. 38) 

“ r  Z  Z  if" 


X  -  ~  --L  (III. 39) 

n  r  r  Z  df  Z  d*r 

It  should  be  borne  in  mind  that  these  base  vectors  are  dif¬ 
ferent  froa  those  used  in  subsection  C  of  this  Appendix. 

The  taeobrane  stress  vectors  associate**  with  the  upper 
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and  lower  ncr.brancs  arc  shown  acting  on  the  volume  decent 

of  Fig.  III. 6  in  Fig.  III. 7.  These  are  depicted  as  ~<rm‘ 

«•  W 

and  acting  on  the  upper  a r.d  lower  membranes,  respec¬ 

tively 


Figure  III. 7 


Since  it  is  assumed  that  there  exists  a  membrane 
state  of  8 tress  in  the  membranes,  the  stress  vectors  on 
the  upper  membrane  are  written  in  component  fora  as 


(no  sum  on  ot  )  (HI. 40) 


Vith  the  introduction  of  the  Weingarten  relations  the 
base  vectors  become 


(III .41) 
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and  hence 


(no  sun 
on  oC) 


(HI. 42) 


Similarly  for  the  lower  membrane 


(no  sum  (111.43) 
on  <*) 


The  equipollent  force  resultant  vector  for  the  mem¬ 
brane  stresses  acting  on  the  length  dl£^  is  defined  by 
the  relation 


where  ^(u.)  and  are  the  thickness  of  the  upper  and 

lower  membranes,  respectively.  Substituting  Eqa .  (111.42) 
and  (HI. 43)  into  Eq.  (III. 44)  yields 

T'JQ  -  (T»-v"£effi.+V~ijgiT  (III.45) 

% 

where  there  has  been  introduced 

.  r  ^  4>  Tm  +  JKl  %l)  C  '  «“■«> 

y“'r»  4;  C~  Ju  *£T  »“•«> 

and  it  is  known  that 
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H  - 

Hu 

H " 

(III. 43) 

The  components  of  the  force  resultant  tensor  are 

defined  by  writing 

f'JF'- 

rrAr +i’# 

(III. 49) 

T4aF  - 

(III. 50) 

By  cocqparing  Eq.  (III. 49)  with  Eq.  (III. 45), 

there 

is  obtained 

i 

[t"-v”£b;J 

(III. 51) 

-L.  j 

vr  1 

f  T*Vflfj 

(III. 52) 

'  y'x  iq 

2»?»; 

(III. 53) 

In  a  similar  manner  the  remaining  components  can  bm 

determined. 

j  Ta-virfs;j 

(III. 54) 

'r*'-v‘*£sjJ 

(111.55) 
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*♦ 


« 


» 


I 

\ 


(111.56) 


The  equipollent  moment  resultant  vector  is  defined  as 


which  becooes  with  the  introduction  of  Eqs.  (III. 42),  (III. 43), 
(III. 46)  and  (111.47), 


This  leads  to  the  following  definitions  for  the 
moment  resultant  tensor 


m’tf' -Ja  mnA*-jA  maA'  (in. 59) 

=  -/T  *U$‘  ~  4a  ma A,  (m. 60) 

A  comparison  of  Eqs.  (III.59)  and  (III. 58)  leads  to 
the  following  components 


*'«  JJfvvr'f  *] 

r"gs;] 

T‘'££'l 


(111.61) 

(111. 62) 

(111. 63) 


223 


FDL  TDR  64-147 


*  vr  r  2 


(111.56) 


The  equipollent  moment  resultant  vector  is  defined  as 

" '  f'  "  *  ft., %  <ln ' 57) 

which  becones  with  the  introduction  of  Eqs.  (III. 42),  (111.43), 
(III. 46)  and  (III .47), 

-  p*{(V‘'-7% "ipv^jjpj  cm. 58) 

This  leads  to  the  following  definitions  for  the 
moment  resultant  tensor 

*  ^ m"  A*- ^/a1  m “a1  (HI. 59) 

"  *Ja  muX,  (in. 60) 


A  comparison  of  Eqs.  (III. 59)  and  (III. SB)  leads  to 
the  following  components 


w'.-Hrv1'-  r’i-’j,'] 

ijri  *  ,J 


(111. 61) 

(111. 62) 

(111. 63) 
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n “  -  T  2rfs;j  (in.M) 

The  pressure  of  Che  gas  on  the  element  of  volume  Is 
dcpiccod  in  Fig.  III. 8  where,  for  purposes  of  clarity,  the 
pressure  is  shown  acting  on  only  one  face. 


It  will  be  necessary  to  define  vectors 

««•«> 

in  order  to  delineate  the  areas  on  which  the  pressures  are 
acting .  ’  Since  the  pressure  is  a  force  which  acts  perpen¬ 
dicular  to  a  surface,  the  equipollent  force  resultant  acting 
on  the  right  face  of  the  volume  element  shown  in  Fig.  III. 8 
is  defined  from  the  equation 
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r 

2  .  *  -ft  * 


P  K s  -J  P ^  IP  (iii-66> 

where  It  should  be  recognized  thee  the  vector  is 

the  unit  vector  which  is  perpendicular  to  the  shaded  area 
dtj  .  Th*  vector  H  *  is  obtained  from  the  vector  pro¬ 
duct  of  7t  and  7T,  : 

rt'M-r&r  <«.«) 

The  substitution  of  Eq.  (111.67)  into  Eq.  (111.66) 
leads  to 

PS.-  -J  Oa.« 


In  a  manner  analogous  to  the  treatment  of  the 
brane  stresses  we  write 


ys-  spA'+spA' 

*Ja"  »  *Ia 


(111.69) 


(III. 70) 


By  cockering  Eqs.  (111.68)  and  (III. 69),  we  arrive 


4? 

p  -  f*  i 

r>  J-4! 
z 


(111.71) 
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t-£ 


(111.72) 


Similarly 


t'-ii  &(-<*)<■  & 


M. 

r 


/* 


if  i  is‘v-° 


(111.73) 


(111.74) 


The  shear  chords  are  in  a  state  of  uniaxial  tension 
caused  by  the  pressurized  gas  and  their  orientation  in 
space  is  along  the  vector  $  .  This  means,  since  the 
element  of  volume  (see  Fig.  111.6)  is  carved  out  with  edges 
orthogonal  to  the  middle  surface,  that  many  of  the  shear 
chords  are  cut  by  the  creation  of  the  element  of  volume. 

We  can  visualize  the  several  shear  chords  as  a  thicket  of 
wires  with  ends  on  the  surface  (see  Fig.  III. 9).  There 
oust  be  placed  on  these  ends  the  forces  representing  the 
stresses  in  the  shear  chords.  It  will  be  assumed  'that  there 
are  a  sufficient  density  of  shear  chords  such  that  their 
effect  can  be  replaced  by  an  effective  stress.  <T  ,  with  a 
line  of  action  along  3  . 

Let  us  isolate  a  tetrahedron  from  the  thicket  of 
shear  chords  (see  Fig.  III.  10).  The  tetrahedron  has  as  a 
base,  the  vectors  ^  dij  and  and  an  altitude  along 

the  vector 
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.0^1 


•IMAGINARY 

MIOOLE 

SURFACE 


•  REPRESENTATION 
of  Thicket* 


Figure  III. 9 


*,d «* 


Figure  III. 10 

The  fourth  face  (canted  face)  is  taken  perpendicular 
to  J  ,  i.e.,  it  is  orthogonal  to  the  shear  chords.  On  this 
fourth  face  the  stress  vector  is  that  caused  by  ehe  tension 
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in  the  shear  chords  and  is  shown  as  rf  where  is  the 
unit  vector  taken  along  the  shear  chords.  The  tetrahedron 
is  in  equilibrium  and  this  is  represented  by  the  equation 


tdA;  -  dt'dAt+  d\jAt  +  7‘JdAJ  (III. 75) 

-f'l  -f, 

where  a  ,d.  are  stress  vectors  acting  on  the  faces  of  the 
tetrahedron  which  are  perpendicular  to  /f  ^//^and  7?  re¬ 
spectively  (see  Fig.  III.  10).  Note  also  that  * t'A*  and 
are  the  forces  (per  unit  length  of  the  aid-surface)  which 
arc  equipollent  to  the  stresses  in  the  severed  shear  chords 
between  the  upper  and  lower  membranes  of  the  element  of 
volume  (see  Fig.  III. 9).  The  areas  of  the  faces  on  which 
these  vectors  are  acting  are  represented  as  d.All(lAZldAi 
with  oLk^  as  the  area  of  the  fourth  face.  Geometrically 


(111.76) 


We  can  express  5  (see  Eq.  (III. 26))  in  terms  of  the  contra- 
variant  base  vectors: 

^  ft*  *  ^  N  (111.77) 

and  hence  it  is  seen  that 

(111.78) 

(111.79) 

(111.80) 
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Then  the  equilibrium  equation  becomes 
<r*s  -  +•  ~2*Ja?s£+  d\  (III. 81) 

This  last  expression  gives  us  license  to  write 
dt'JJF  =*  d'^A*  +  <^3)l7^  (III. 82) 

c/*^  *  <JZHAoi’t  d(lHN  (III. 83) 

Next  we  isolate  an  element  of  volume  from  the  thicket 
with  one  edge  along  and  two  other  edges  perpendicular  to 
£  (see  Fig.  III. 11). 


V 


I 


Figure  1II.U 

The  vectors  t and  ^  are  unit  vectors  which  together 
with  5  form  a  emtually  orthogonal  set.  The  fourth  or 
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< 
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In  a  similar  fashion 

2Z~ 


(III. 59) 


(III. 90) 


By  cocabining  Eqs.  (Ill  89)  and  (III.90)  with  Eqs. 
(III. 82)  and  (III. 83)  we  obtain 


<rV-s  -  +  d(3)/*N 


(III. 91) 


Thus,  by  using  Eq.  (111.26),  we  obtain 

jflct  »  B  oi 
<L  —  O'  SrS 


(III. 92) 


<r*s^ - 


(111.93) 


E.  EQUATIONS  OF  EQUILIBRIUM 

There  has  been  defined  in  the  previous  section  force, 
resultants  and  moment  resultants  which  are  equipollent  to 
the  stresses  in  the  two  membranes,  the  internal  pressure, 
and  the  shear  chords.  This  equipollent  system 'is  defined 
for  the  middle  surface  of  the  shell  and,  therefore,  the 
equilibrium  equations  are  written  for  an  element  of  area 
on  the  imaginary  middle  surface.  It  is  convenient  to  de¬ 
fine  quasi- force  and  moment  vectors  as  follows: 

Lr  **  -fA  [I  Am+  (III. 94) 
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(III. 95) 

P-4T  g  A* Af 

(III. 96) 

(III. 97) 

The  conversion  of  Che  preceding  vectors  to  forces  is 
accomplished  sinp ly  by  multiplying  the  appropriate  d^ e.g., 
the  force  acting  on  the  mid-surface  due  to  4*  is  Lt ' d-g* • 

We  also  define  a  loading  force  vector  F  as 

r.jx-frK^vj  (hi. 9a> 


All  of  these  are  shown  in  Figs.  III.  12  and  III.  13  (two 
figures  arc  used  for  purposes  of  clarity) . 


The  vector  equation  of  force  equilibrium  is  obtained 
by  summing  all  the  force  vectors  in  Figs.  III. 12  and  III. 13. 
There  results 


+ 


fJa-o 


(111.99) 


The  vector  equation  of  moment  equilibrium  is  obtained 
by  summing  the  moment  vectors  (represented  by  double-headed 
arrows  in  Fig.  III.  13)  and  taking  moments  of  all  the  force 
vectors  about  point  o  .  The  result  is 

fpr.+  Arx  £r)~0  (III.100) 
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We  make  use  of  Che  following  auxiliary  relations: 


w  nr-  n)i  } 


+  Va  +  dt*rtf) 

df  1 


>fr » S%  U‘r$Jt^k’iJr} 


V 


2_*f‘ 


(III. 102) 


(111.103) 


(III. 104) 


(HI.  105) 


(III. 106) 


(III. 107) 


These  latter  relations  lead,  when  substituted  into  the 
vector  Eqs.  (111.99)  and  (III. 100),  to  the  following  tensor 
equations  of  force  and  moment  equilibrium: 


r\r-{frq\-rB^r.$£A~+F--o_  mi.™. 


FDL  TDR  64-147 


F.  STRESS- STRAIN  RELATIONS  FOR  SMALL  STRAIN 

The  upper  and  lower  membranes  of  an  inflated  shell 
will  most  likely  be  constructed  of  a  woven  fabric.  Although 
the  microscopic  behavior  of  such  a  material  is  very  difficult 
>  to  specify,  the  gross  macroscopic  behavior  can  be  satis¬ 

factorily  approximated  by  treating  the  membrane  as  a  material 
with  two  axes  of  symmetry.  The  3tress-strain  relation  can 
^  then  be  written  in  the  general  form 

E“^evrez  (in.  in) 

for  general  curvilinear  coordinates.  We  will  call  the 
E'KP&X  the  elasticity  tensor  and  it  obeys  the  symmetry 
relations 

E^a £^Sr=.  E^rB-E  <9r"1  (UI.U3) 

*  • 

If  we  choose  a  set  of  local  rectangular  cartesian  coordinates 
{  to  coincide  with  the  axes  of  orthotropy,  the  stress-strain 

relations  become 

c  £5C  (III.  114) 
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<  (mP*er+r\*.'cL*'r)=o 

acr  \  P  ' 


(m.iio) 

(III. Ill) 


* 
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F.  STRESS- STRAIN  RELATICUS  FOR  SMALL  STRAIN 


The  upper  and  lower  membranes  of  an  inflated  shell 
will  most  likely  be  constructed  of  a  woven  fabric.  Although 
the  microscopic  behavior  of  such  a  material  is  very  difficult 
to  specify,  the  gross  macroscopic  behavior  can  be  satis¬ 
factorily  approximated  by  treating  the  membrane  as  a  material 
with  two  axes  of  symmetry.  The  stress-strain  relation  can 
then  be  written  in  the  general  fora 


(III. 112) 


for  general  curvilinear  coordinates.  We  will  call  the 
£*^*r  the  elasticity  tensor  and  it  obeys  the  symmetry 
relations 


gfiotev ^  ^  ^  BTatfl 


(III. 113) 


If  *•  choose  a  set  of  local  rectangular  cartesian  coordinates 
to  coincide  with  the  axes  of  ort  hot  ropy,  the  stress-strain 
relations  beeoae 


G~otP 


(III. 114) 
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where  <3*°^  and  £9T  are  physical  coaponents.  Since  the 
material  is  assumed  orthotropic,  there  are  only  five  inde¬ 
pendent  components  of  the  elasticity  tensor. 


g£Mil 


-a^OV 


t-12.it  * 


r 


(III. 115) 


If  we  indicate  the  transformation  from  rectangular 
cartesian  coordinates  ^  ,  to  general  curvilinear  coordinates 
5*.  as 


(III. 116) 


then  the  elasticity  tensor  in  the  £  coordinates  is  trans¬ 
formed  to  the  coordinates  as 

EmP*v =s  13*  (III.  117) 

For  an  isotropic  material  the  elasticity  tensor  takes 
on  the  form 


where 
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I 


< 
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/*  ** 

with  £  and  *>  being  the  usual  Young's  modulus  and  Poisson's 
ratio,  respectively.  Tlio  are  the  octric  tensors  of  the 
membrane  surface. 


G.  INCREMENTAL  THEORY  LEADING  TO  URGE  DISPLACEMENTS 

The  equations  which  have  been  developed  for  the 
large  deformations  of  an  inflated  shell  are  much  too  dif¬ 
ficult  to  solve  explicitly.  Wie  will  propose  a  method  of 
solution  which  arrives  at  the  final  result  by  a  scries  of 
linearized  equations. 

In  Fig.  III. 14  the  imaginary  middle  surfaces  at 
three  different  stages  of  deformation  arc  shown.  The  first 
stage  is  the  undofotmed,  the  second  will  be  called  the  de¬ 
formed,  and  the  third  will  be  called  the  further-deformed. 

It  will  be  assumed  that  the  additional  displacements  leading 
to  the  further-deformed  shell  are  sufficiently  small  so  that 
the  squares  and  products  of  changes  in  all  quantities  can 
be  neglected. 

Let  us  first  define  and  calculate  the  geometric 
parameters  of  the  imaginary  middle  surface  of  the  further- 
deformed  shell.  The  displacement  vector  connecting  the  de¬ 
formed  and  further-deformed  middle  surfaces  will  be  denoted 
by  where  the  A  signifies  a  small  change.  All 

vectors  and  tensors  will  be  denoted  as  the  tensor  of  the 
deformed  shell  plus  a  A  change,  e.g.,  the  base  vectors  are 
denoted  by  4/?^  ,  the  fundamental  metric  tensor  by 

etc.  We  assume  that  squares  of  A  changes 

are  negligible. 

Our  primary  interest  lies  in  the  A  changes.  For 
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Figure  111.14 
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exanp le 


and  hence 


5  it  ?dui 
2$*  + 


(111.119) 


(III. 120) 


We  will  write  in  terms  of  components  with  respect 
to  vector  directions  of  the  deformed  middle  surface 

&U.  -  Ap  +  AuJ-N  (Ill. 121) 

This  leads  to  the  following  relation  for  A  A, ( 

M,-  [ ^'j- [A“-'en*  %?]« 

where  the  vertical  bar  signifies  covariant  differentiation 
with  respect  to  the  metric  of  the  deformed  middle  surface. 

•  We  recall  the  basic  definition  for  the  base  vectors 

A^+aA^  -  (A,+aX)-(y  (iu.123) 

As  a  consequence  of  the  assumption  regarding  the 
magnitude  of  the  A  changes,  we  neglect  the  term 
and  approximate  AA^  as  follows: 

A**.’***  (III.  124) 


♦ 
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In  a  similar  fashion,  the  changes  in  the  parameters 
of  the  strained  middle  surface  can  be  found. 


A  A‘'m  -  /t'V1*  AA?f 

(III. 125) 

A  A  -  AA~aA„ 

(III. 126) 

AU  -gfi'Z  +  Xn*] 

(III. 127) 

AA*  =.  A**aAa  + 

(III. 128) 

AB  «.  AN-  A/* 

(III. 129) 

A  a;* 

(III. 130) 

a/*1  «  AA'.  +  A“- 

IfiJ  w . 

(1X1.131) 

The  increments  in  the  components  of  strain  on  the 
upper  membrane  are  defined  in  the  traditional  manner: 


(III. 132) 


where  the  index  U  denotes  the  upper  membrane.  A  similar 
expression  can  be  written  for  the  lower  meaibrane. 


The  relationships  of  the  various  vectors  connecting 
on  the  upper  membrane  of  the  deformed  shell  to  the  saaw 
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point  on  £'  of  the  further-deforced  shell  are  depicted  in 
Fig.  111.15. 


t  The  vector  joining  O'  on  the  imaginary  middle  surface  to 

£'  is  seen  to  be  denoted  by  where 

At  -  AsTAr  +  (in.  133) 

We  observe  that  the  position  vector  to  £  is  given 
hy 

K+  &  (111.134) 

and  hence 

d(8u+  +  cL(& %+At)  (III .  135) 
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This,when  substituted  into  Eq.  (Ill.  132),  leads  to  the 
following  expression  for  the  increments  in  the  components 
of  strain  on  the  upper  membrane 


wherein  we  have  dropped  terms  involving  products  of  A 
changes. 

In  a  similar  manner  the  components  on  the  lower 
membrane  can  be  obtained 


t* 

(in.  137) 

It  will  be  convenient  to  define  an  average  strain 

Ay* 

and  «  curvature  strain  • 

r  *  %  / 

(III. 138) 

Ai’  - 

<S»  *  l  H  f' 

(III. 139) 
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The  substitution  of  Eqs.  (III.  136)  and  (111.137)  into  Eqa. 
(111.138)  and  (III. 139)  yields 

A/*  .  —  Hi}  (111.140) 

Ak0 

z\  «  iff  T  fi  ds« 


-a.2tf  I 


(III. 141) 


wherein  the  nonlinear  terms  have  been,  omitted  from  Eq. (HI.  140) 
The  derivatives  of  s  ,  Ait  ,  end  At'  are  as  follows: 


ds 


[s%-By]t +[’%*§]*  <ut'i42) 


- (m  u3> 


g-M-Zty+yX+’SP  «*•«*> 


There  results  for  the  components  of  strain 
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-f  /  4a  /  +4tf /-££  4«r 
*  *1  *lfi  ?'*  1® 


d*y-i 

*  [s°l  -Bp][  **]*„■ 

+[s%  -<<M  ~BJ**']4rr 

4r*^]hX+'#] 
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The  first  three  terms  for  in  Eq.  (111.145)  are 

recognizable  as  of  the  same  form  as  those  appearing  in  the 
infinitesimal  theory  of  thin  shells.  The  remainder  of  the 
terms  in  Eq.  (111.145)  can  be  shown  to  be  proportional  to 
the  square  of  the  spacing  between  the  membranes.  Again, 
these  terms  are  usually  omitted  from  thin  shell  theory.  It 
is  not  clear  at  this  stage  that  omission  from  the  theory  of 
inflated  shells  is  justified  since  some  of  the  proposed  con¬ 
figurations  are  not  "thin.”  Similarly,  only  the  first  two 
terms  for  in  Eq.  (III. 146)  are  retained  in  thin  shell 

theory. 

Tfcs  increment  ip  the  transverse  shearing  strain 
(see  Eq.  (III. 37))  is  defined  by  the  equation 

(HI. 147) 


v  from  which  there  is  obtained 

•  i"  (%*'  +  (III. 148) 

In  terms  of  the  increments  in  5  and  tL  ,  Eq.  (III.  148) 
becomes 

% 

4f 

(  The  Inflated  sandwich  shell  remains  in  equilibrium 

in  its  further-deformed  state.  We  express  the  force  and 
moment  resultants  acting  on  the  middle  surface  in  the 
J.  .  further -deformed  state  as  the  force  and  moment,  respectively. 
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in  the  deformed  state  plus  a  small  increment.  Thus,  the 
two  equilibrium  equations  given  by  Eq.  (111.108)  become 


nK+O+r*  ^)nr 


wherein  the  double  vertical  bars  preceding  a  subscript 
signifies  covariant  differentiation  with  respect;  to  the 
metric  of  the  further-deformed  middle  surface.  For  example, 
the  covariant  derivative  of  the  force  resultant,  /"W™ 
with  respect  to  the  metric  of  the  further -deform*-  .ate 
is  given  by  * 


If  Eq.  (III. 151)  is  substituted  Into  Eq.  (111. ISO) 
together  with  an  analogous  expression  for  ( d'rm+  A**")  11^. 
there  will  result,  since  Eq.  (III.  108)  can  be  factored  out. 
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Similarly,  the  equation  of  force  equilibrium  in  the 
normal  direction  (sec  Eq.  (III. 109))  and  .the  moment  equi¬ 
librium  equations  (see  Eq.  (III. 110))  become 


*flr- 

-*•  ^  ^  HI .  1S3) 

in. *“l  (III.1J4) 

The  increments  in  tensor  quantities  «  and  CL 
which  are  related  to  the  stress  in  the  shear  chords  need  to 
be  explicitly  displayed.  With  Eq.  (111.92)  as  the  model 
we  write 


249 


k 


) 


]  >■>  ! 


and  hence 


if*. 


*  /  «  y 

<r  /«*•*.*« 


►  .*  j 


(ni.ijij) 


H.  FINITE  DEFLECTION  OF  AN  INITIALLY  FLAT,  INFLATED, 

UNIFORM  SANDWICH  PLATE 

In  this  section  the  equations  derived  thus  far  will 
be  applied  to  an  Initially  flat  inflated  sandwich  plate. 

We  will  assume  the  deflections  are  finite  which  aeons  the 
rotations  arc*  of  a  magnitude  such  that  the  else  of  Che 
angle  of  rotation  can  be  replaced  by  the  angle.  fVirthensora, 
only  those  rotations  occurring  about  the  base  veceors  2^ 
are  finite;  rotations  abwut  the  normal  are  Infinitesimal  in 
magnitude.  The  plate  is  rectangular  in  shape  and  will  ha 
described  by  coordinates  u*  as  shown  in  Fig.  XIX. Id. 
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CoW  available  to  WIC  does  not 
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Figure*  III .  io 
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Copf  available  to  DTIC  does  not 
permit  full*  legible  iepraductioa 


(111.162) 

u  —  Urir+  itrtj 

(III. 163) 

t  m  +TS 

(111.164) 

(III. 165) 

V  X  %«) 

(III. 166) 

ta*  i  {s.  +K.) 

(III. 167) 

The  following  consents  are  epproprlata: 

a*  The  5J  coopooent  of  ^  (see  Eq.  (111.26)  has 
bean  sec  equal  to  unity  la  order  to  be  consistent  with  the 
assumption  of  small  rotations. 

*  Only  the  nonlinear  tans  Involving  the  finite 
■*d  by  WJ*  are  retained. 
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r--*[*£*rS]  (III-l77) 
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(in.  i79) 

»i*  ■*  rAv  (in.  180) 

/"=  -tcf'^N  (in. 181) 

=  £***  (111.182) 

fcs?S*  (III.  183) 

^  •  fa**  (III.  184) 

&*f,mkL  an.usj 
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The  equilibrium  equations  assume  the  fens  (ace  Eqs. 
(III. 108)  -  (III. 110), 
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vtwra  <ht  la  t ha  laltlal  |twtrm  wa 
pmant  <£  .  feta  that  alaaa  /J*  la  a 
forca  aoalllhriaa  aaaatlaaa  haaaKf 


(in.  if  3) 


Aa  atraw  atrala  eolation  far  da  « 

wrlctao  la  tha  fan  (aaa  If.  (m.112)), 

4r^.«^r5r 


(ns.  if?) 


abaca  la  tha  alaatlaity  faaaar  apt* 
axtatlac  atata  of  atcaaa.  far  tha  flat  ani 
tha  aan  rhlahnaaa  aaffaaa  tcaatai  la  thla 
following  clap  la  rotation  aaa  ha  nai  la  fli 

(111.197). 


flaaaaf  ft* 


fair  /* 

it 


(in.tM) 


9 


kl  .  jq,  .  rm  *  ****** 

total  oi  18  oaten—.  For  tte  oolatloa  m  kart  tie  mpm 
i  «tma-di#kHMBe  nlatiat  $Um  ky Ip.  (XlJ.ttS)- 

(III.  117),  tte  21m  ofulUWl—  afoattaao  gfc—  ky  ft)*. 
(III.  153),  (ZZZ.191)  ml  (m.lftk),  oaf  tte  rti  mm 
« train  relation*  giraa  kjr  te«  (XZX.XM)  -  (HZ.1M)>  . 

It  ateuld  te  >f —of  tte*  ml  of  tte  etplUkrt— 
•qutloM  cm  te  paseltkte  tte o  tkroo  groove  of  Um. 
Om  group  ooMloto  of  oteagoo  te  tte  fteoo  (or  aomt) 
gooaltooto  coupled  rngotter  rite  tte  gaoontry  of  tte  otell 
la  lto  fofoiaod  iuu,  o.g .,  fr%,Y8'  •  . 

group  oonoioto  oc  too  soroo  (or  ONONC/nouton  oc  too 
dofo—1  otteo  oouplof  taootter  with  tte  oteagoo  te  tte 
got— try,  o.g.,  «r  .  Ite  tkiri  group  am 

tte  tee romteo  te  tte  applied  loefe,  If*  end  If4  .  Stem 
^  tte  teflteotf  ooadvleh  otell  ooa  defleet  uoro  seedily  te--  tea 

nnrual  dlraetioa  tten  te  tte  toogwMol  direction,  our  te- 
tultleu  loodo  no  to  o^oet  Ctet  tte  tone  loeohriog  oteagoo 
la  tte  >00001  ftaaluntol  toaoor  am  aoco  log— oof  tkaa 
tteoo  lavolrteg  nhoagoo  1a  tte  ftnritoffol  oyatels.  Ite 
It  a ay  te  Mgooff  tte  aoay  of  tea  tamo  aay  te  aagleeted 
ao  they  ora  la  tte  tteory  of  tkte  otello.  ten r,  aatil 
oa  actual  configuration  la  toofclof,  tte  tolatlro  ortfon  of 
aogoltudo  ooaaot  te  oaaerotoly  dotoraiaof. 


i 
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SYMBOLS  (Cont'd) 


-v> 

Poisson's  ratio 

%)  >  °(il 

Stress  vectors  on  upper  and  lower 

membranes 

<r* 

Effective  stress  of  shear  chords 

T**  T** 

C(*)  *  L(i) 

Stress  tensor  on  upper  and  lower 

membranes 

Associated  v/ith  Undeformed  Middle  Surface 

£ 

Covariant  base  vectors 

a  a.m* 

*0  > 

First  fundamental  tensor 

v* 

Second  fundamental  tensor 

u 

Unit  normal  vector 

t(«W) 

Position  vector 

Christoffel  symbol  of  second 

Associated  with 

Deformed  Middle  Surface 

*< 

Covariant  base  vectors 

\.A+ 

First  fundamental  tensor 

Second  fundamental  tensor 

n 

Unit  normal 

*0 

Position  vector 

Christoffel  symbol  of  second  V<nd 
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Asserted  Marks 

H coasted  Greek  Indices  to  be  summed  from  1  to  2. 


(  ), 


Signifies  covariant  differentiation 
with  respect  to  metric  of  undefoxned 
middle  surface,  e.g.. 


*  <r 


<m 


c  > i . 


Signifies  covariant  differentiatica 
with  respect  to  metric  of  deformed 
middle  surface,  e.g.. 


n.-*?' */«/'"  W' 


o-r 


Vector  product 
'Scalar  product 
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11.  Wm.UUTU1  mt  t*<> 


Ift  MHUCT 

Ibis  rsport  cows  •  pacggi  of  research  on  analytical  etndlas  sat  saqwrt- 
orotal  verification  of  ths  vibration  and  flatter  oharacterlstlea  sf  bfliOU 
atrueturee,  ss  •  first  stop  tosard  wtatoUablag  Assign  criteria  ahdob  will  lnneri 
that  lifting  surfaces  for  Tutors  lnflstsbls  vehicles  will  bo  free  fro*  djrnnlu 
Instabilities.  tbs  rososroh  effort  was  aalaly  sosflssi  to  tho  IMiT  -  typo 
•true two,  wtaoro  tho  spacing  of  tbs  ewfboes  Is  mM»w4  h?  darts  halt  1*  tension 
by  intend  pressure* 

A  asthed  Is  presented  for  noaatraottag  18KU  aodols  which  are  suitable  for 
use  as  wird  tunnel  flutter  aodols.  In  edition  ta  tbs  stented  iftiaUa  Ante, 
the  vibration  node  ■hapee  of  s  Aslts  rlsg,  as  AetscatneA  by  «*e  of  atrrors  novated 
on  tbs  nodal,  are  presented,  low  opsrtaestal  results  on  tb»  pregartle*  of  tbs 
aoAel  surface  naterld  are  also  IneludaA. 

Maty-five  flutter  points  vers  obtalaeA  for  osntllevereA,  half  open  delta, 
wing  nodule  over  the  Mftch  sober  mgf  fron  0.116  ts  6.0.  9a  Sddltlsn  be  Aster- 
■living  an  eaperlnental  boundary  frr  a  basis  Aslts  sr*Tl,  tbs  effasts  ef  landing- 
edgs  rjaep,  eagle  of  attack,  sad  ths  pressure  differential  between  the  nodal  and 
ash  lent  air  preasures  ware  lnresil^ted  at  anwanl  agninai  Hash  wabars. 

A  tbaoretlsal  laraatlgatlon  eUbeaas  type  flatter  of  an  inflatable  penal 
and  the  development  of  the  egwtlene  governing  tbo  shau  tarsi  behavior  of  aa  In¬ 
flated  AZMAT  Motion  are  presetted  U  segarate  oggaadloee; 


